计算连续干天日数

大海绵1990

我的数据是17个站点分开的，每个站点里面的数据是这样的，第一列是年，第二列是月，第三列是日，第四列是降水量:
1960    1       1       0
1960    1       2       0.6
。。。。
2011    12    31      0
我想要计算每年持续时间最长的干天（日降水量<1mm）日数，所以我写了这样的程序，刚接触matlab所以只能用最简单的语法来拼凑，
for n = 1:17
    data = prec_all{n,1};
% assume that your data contains 52 years and started from 1960
TDim=52; startYear=1960;
Result=zeros([TDim,1]);
% count rainfall days by a for loop
   for i=1:TDim
    initial=datenum(startYear+i-1, 1, 1)-datenum(startYear, 1, 1)+1;
    final=datenum(startYear+i-1, 12, 31)-datenum(startYear, 1, 1)+1;
     m=final-initial+1;
     extract=zeros([m,1]);
   extract(1:m,1)=data(initial:final,4);
   for z=1:m
       a=0;
       days=[];
   if extract(z,1)<1
       a=a+1;
   else a=a;
   days=a;
   end

   end
   end
end
程序是不对的，从第三个for循环下来结果就不对了，高手们能帮我指点下吗？有更好的办法能告诉我吗？拜托拜托！

二爷名声在外

我给楼主一个思路吧，不知道准确不，仅供参考。
楼主可以做一个17*总日数的数组，对每一个降水量小于1的日期都进行计算其后的干燥日数（如果大于1就计为0），然后每一列选取最大值就好了。
计算干燥日数：
for i=1:length-1(length:总日数）
if rain(i,4)<1
for j=i+1:length
if rain(j,4)>=1
break;
days(i,n)=j-1+1;
end
end
end

大海绵1990

二爷名声在外 发表于 2014-7-7 13:07
我给楼主一个思路吧，不知道准确不，仅供参考。
楼主可以做一个17*总日数的数组，对每一个降水量小于1的日 ...

谢谢你啊，给了我启发

唐屠夫

很实用  不错   

温文尔雅

貌似改一改很有用啊~

rjys1219

新手学习中。。。。

whuliuqiushi

谢谢分享，好东西

小叶子

思路也可以Fortran