零点法计算周期

【零点法计算周期】我承认这程序写的毫无美感可言,可是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)

    推荐阅读