diff --git a/atovs_read_find.m b/atovs_read_find.m
index a0eb5e95a844eecf2aee48b9ba186efa67ee18ba..02da2c64a76ba001565dae0fb8a96a4add2ae1d8 100644
--- a/atovs_read_find.m
+++ b/atovs_read_find.m
@@ -12,6 +12,7 @@ ret.mhsfile = filename;
 ret.x = -1;
 ret.y = -1;
 ret.maximum = -1;
+ret.maxcounts = -1;
 ret.lon=-1;
 ret.lat=-1;
 ret.height=-1;
@@ -139,15 +140,15 @@ if intr_end < intr_start
     disp('Not enough data in file')
     return
 end
-y1 = double(squeeze(data.dsvcounts11));
-y2 = double(squeeze(data.dsvcounts21));
-y3 = double(squeeze(data.dsvcounts31));
-y4 = double(squeeze(data.dsvcounts41));
+y{1} = double(squeeze(data.dsvcounts11));
+y{2} = double(squeeze(data.dsvcounts21));
+y{3} = double(squeeze(data.dsvcounts31));
+y{4} = double(squeeze(data.dsvcounts41));
 for i = intr_start:intr_end
-	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;
+	A(1,i) = median(y{1}(i-5:i+5))-(median(y{1}(i-70:i-60))+median(y{1}(i+60:i+70)))/2;
+	A(2,i) = median(y{2}(i-5:i+5))-(median(y{2}(i-70:i-60))+median(y{2}(i+60:i+70)))/2;
+	A(3,i) = median(y{3}(i-5:i+5))-(median(y{3}(i-70:i-60))+median(y{3}(i+60:i+70)))/2;
+	A(4,i) = median(y{4}(i-5:i+5))-(median(y{4}(i-70:i-60))+median(y{4}(i+60:i+70)))/2;
 end
 M(1) = max(A(1,:));
 M(2) = max(A(2,:));
@@ -160,6 +161,7 @@ if size(ret.x, 1)>1
     ret.y = ret.y(1);
 end
 ret.maximum = repmat(maximum, size(ret.x));
+ret.maxcounts = y{ret.x}(ret.y);
 ret.lon=((double(data.lonfov45(ret.y))+double(data.lonfov46(ret.y)))/2E4)'; % longitude of nadir at maximum of light curve
 ret.lat=((double(data.latfov45(ret.y))+double(data.latfov46(ret.y)))/2E4)'; % latitude of nadir at maximum of light curve
 ret.height=(double(data.height(ret.y))/10)'; % altitude of satellite at maximum of light curve