NCL-有没有办法从索引信息中减去相同的变量、不同的值（相同的文件）？

Question

我有一个 3D nc 变量 PP = (time, lat, lon) 一整年的累积降水量，我想计算每日值。

我想在示例中做一些链接，但找不到任何示例。

我试图做一个循环，其中减法应该在 [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 提议，欢迎他们。

Answer 1

您可以使用两个 NCO shell 命令执行此操作

下面的命令将记录上移一位

ncks -v TOTALRAIN --msa_usr_rdr -d time,1,364 -d time,364 in.nc shift.nc

现在使用 ncbo 找出差异

ncbo -v TOTALRAIN shift.nc in.nc diff.nc

现在 netcdf 文件 diff 应该包含您的每日 pp 此文件中 TOTALRAIN 的最终记录将为零

...亨利

Answer 2

在 cdo 中，您可以使用分解累积变量

cdo deltat in.nc diff.nc 

这与此问题有关，请参阅那里的答案以获取更多详细信息，包括插入初始时间步长。