零点法计算周期

我承认这程序写的毫无美感可言，可是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)