我承認(rèn)這程序?qū)懙暮翢o美感可言，可是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)
