2

我有一个全年累积降水的 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提议，欢迎他们。

4

2 回答 2

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 文件差异应该包含您的每日 pp 此文件中 TOTALRAIN 的最终记录将为零

...亨利

于 2019-03-31T10:35:46.347 回答
0

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

cdo deltat in.nc diff.nc

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

于 2021-10-21T20:39:05.520 回答