diff --git a/atovs_read_find.m b/atovs_read_find.m index bc7b76c2a613e1178b0838e7c7534e25bb449832..76062c4e23cb86d66c6ca35f7370a267c0bf070a 100644 --- a/atovs_read_find.m +++ b/atovs_read_find.m @@ -144,10 +144,10 @@ y2 = double(squeeze(data.dsvcounts21)); y3 = double(squeeze(data.dsvcounts31)); y4 = double(squeeze(data.dsvcounts41)); for i = intr_start:intr_end - A(1,i) = median(y1(i:i+10))-median(y1(i+60:i+70)); - A(2,i) = median(y2(i:i+10))-median(y2(i+60:i+70)); - A(3,i) = median(y3(i:i+10))-median(y3(i+60:i+70)); - A(4,i) = median(y4(i:i+10))-median(y4(i+60:i+70)); + A(1,i) = median(y1(i:i+10))-(median(y1(i-70:i-60))+median(y1(i+60:i+70)))/2; + A(2,i) = median(y2(i:i+10))-(median(y2(i-70:i-60))+median(y2(i+60:i+70)))/2; + A(3,i) = median(y3(i:i+10))-(median(y3(i-70:i-60))+median(y3(i+60:i+70)))/2; + A(4,i) = median(y4(i:i+10))-(median(y4(i-70:i-60))+median(y4(i+60:i+70)))/2; end M(1) = max(A(1,:)); M(2) = max(A(2,:));