我有一个3D nc变量PP =全年的累积降水量(时间,纬度,经度),我想计算每日值。
我想在示例中做一些链接,但找不到任何示例。
我试图做一个循环,其中减法应该在[index + 1]减去[index]的值之间进行减法(如下所示)
t = f->time(:) ;size = 365
pp = f->TOTALRAIN(:,:,:) ; time, lat, lon dimensions
do i = 0,dimsizes(t)-1
DailyPp = pp(i,:,:) - pp(i-1,:,:)
end do
但是那样会显示错误。
我需要每天的值,并且只将降水量作为输入信息,并认为这样的循环应该是最好的方法。
如何从变量中进行计算?
我正在使用ncl,但是如果您有cdo或nco命题,欢迎他们。
答案 0 :(得分:1)
您可以使用两个NCO Shell命令来完成此操作
以下命令将记录上移一个
ncks -v TOTALRAIN --msa_usr_rdr -d时间1,364 -d时间364 in.nc shift.nc
现在使用ncbo找出差异
ncbo -v TOTALRAIN shift.nc in.nc diff.nc
现在netcdf文件diff应该包含您的每日PP 该文件中TOTALRAIN的最终记录为零
...亨利