From 29a8db1eda1a5496bbb336f57d67766e6c50115f Mon Sep 17 00:00:00 2001
From: Oliver Lemke <olemke@gmail.com>
Date: Tue, 14 Jan 2025 09:43:38 +0100
Subject: [PATCH] Take median on both sides

---
 atovs_read_find.m | 8 ++++----
 1 file changed, 4 insertions(+), 4 deletions(-)

diff --git a/atovs_read_find.m b/atovs_read_find.m
index bc7b76c..76062c4 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,:));
-- 
GitLab