Take median around i

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-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;
+	A(1,i) = median(y1(i-5:i+5))-(median(y1(i-70:i-60))+median(y1(i+60:i+70)))/2;
+	A(2,i) = median(y2(i-5:i+5))-(median(y2(i-70:i-60))+median(y2(i+60:i+70)))/2;
+	A(3,i) = median(y3(i-5:i+5))-(median(y3(i-70:i-60))+median(y3(i+60:i+70)))/2;
+	A(4,i) = median(y4(i-5:i+5))-(median(y4(i-70:i-60))+median(y4(i+60:i+70)))/2;
 end
 M(1) = max(A(1,:));
 M(2) = max(A(2,:));