我承认这程序写的毫无美感可言,可是ncl这种语言就不配拥有美感。
omega=precip_I(0,:,:)*0
printVarSummary(omega)
do i=0,nlat-1
do j=0,nlon-1
do k=0,ntime-2
if(.not.ismissing(precip_I(k,i,j)))then
if((precip_I(k,i,j)*precip_I(k+1,i,j)).lt.0) then
omega(i,j)=omega(i,j)+1
end if
end if
end do
end do
end do
precip_mean_cycle = ntime*2/omega
copy_VarCoords(precip_I(0,:,:),precip_mean_cycle)
printVarSummary(precip_mean_cycle)