如何使用ncl、cdo或nco及时分解累积的netcdf变量？

Question

我有一个3DNC变量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时间，1,364 -d时间，364 in.nc shift.nc

现在使用ncbo来找出区别。

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

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

...Henry

Answer 2

在cdo中，您可以使用以下命令分解累积的变量

cdo deltat in.nc diff.nc 

这与this question相关，更多详细信息，包括插入初始时间步，请参考那里的答案。