Thread Subject: find function not reporting all column values

Hi,

I have two datasets: one on elevation and one on tidal height and want to create a matrix that will give a value of 1 for when the tidal height exceeds elevation and a 0 for when the tidal height is less than or equal to the elevation. I had a colleague help me to write a script and it works when I run it on two of my datasets, but when I run it on my other datasets, I seem to be missing some values. For example, one of my elevation datasets has 441 entries and my tidal series is 11090 entries. When I run my script however, instead of getting something that's 441x11090, I get something that is 359X11090 and I'm not sure where the other 82 rows have disappeared to. I've pasted in the script I'm using below. Any suggestions on what might be causing the problem or suggestions for solutions would be greatly appreciated. Thanks.

load Mugu_1_Elevations.txt
load Mugu_M1_0_11012009.txt
date=Mugu_M1_0_11012009(:,1);
ht=Mugu_M1_0_11012009(:,7);
elev=Mugu_1_Elevations(:,4);
elev_Mugu=zeros(length(elev),length(date));
for
i=1:1:length(date)
ind=find(elev<ht(i));
inun_Mugu1(ind,i)=1;
end

On Nov 9, 3:01 pm, "Valerie " <valerie.a.c...@gmail.com> wrote:
> Hi,
>
> I have two datasets: one on elevation and one on tidal height and want to create a matrix that will give a value of 1 for when the tidal height exceeds elevation and a 0 for when the tidal height is less than or equal to the elevation.  I had a colleague help me to write a script and it works when I run it on two of my datasets, but when I run it on my other datasets, I seem to be missing some values.  For example, one of my elevation datasets has 441 entries and my tidal series is 11090 entries.  When I run my script however, instead of getting something that's 441x11090, I get something that is 359X11090 and I'm not sure where the other 82 rows have disappeared to.  I've pasted in the script I'm using below.  Any suggestions on what might be causing the problem or suggestions for solutions would be greatly appreciated.  Thanks.  
>
> load Mugu_1_Elevations.txt
> load Mugu_M1_0_11012009.txt
> date=Mugu_M1_0_11012009(:,1);
> ht=Mugu_M1_0_11012009(:,7);
> elev=Mugu_1_Elevations(:,4);
> elev_Mugu=zeros(length(elev),length(date));
> for
> i=1:1:length(date)
> ind=find(elev<ht(i));
> inun_Mugu1(ind,i)=1;
> end

Why did you change the name of the output matrix from elev_Mugu to
inun_Mugu?

This line here:
> elev_Mugu=zeros(length(elev),length(date));
is very important because it sets all the data to zero, then the for
loop resets those data that satisfy your criterion to 1. By changing
the name, you have completely screwed up the routine your colleague
wrote for you. You need to either change the name back or to change
the above line to:
inun_Mugu=zeros(length(elev),length(date));