附一段自己写的连续干日数的程序,希望你能用上
do i=1,ni
do j=1,nj
k=0
do yri=1,39
kk=k+dayi(yri)
do day=k+1,kk
d=day-k
if(var(i,j,day).lt.1)then
n(i,j,d,yri)=1
do day2=day+1,kk
if(var(i,j,day2).lt.1) then
n(i,j,d,yri)= n(i,j,d,yri)+1
else
exit
endif
enddo
else
n(i,j,d,yri)=0
endif
enddo
k=k+dayi(yri)
enddO
enddo
enddo
print*,'n(1,1,1,1)=',n(1,1,1,1)
do i=1,ni
do j=1,nj
k=0
do yri=1,39
kk=dayi(yri)
nn(i,j,yri)=n(i,j,k+1,yri)
do day=k+2,kk
if(nn(i,j,yri).lt.n(i,j,day,yri))then
nn(i,j,yri)=n(i,j,day,yri)
endif
enddo
enddo
enddo
enddo
print*,'nn(1,1,1)=',nn(1,1,1)
do i=1,ni
do j=1,nj
do yri=1,39
nnn(i,j)=nnn(i,j)+nn(i,j,yri)
enddo
nnn(i,j)=nnn(i,j)/39.0
enddo
enddo |