使用CDO将累积变量转换为netcdf文件中的时间步长值

Question

我有一个netcdf文件，在一个带有一个变量的网格上有大约100个timesteps，它是在timesteps上累积的。我现在感兴趣的是计算每个时间步对变量值的贡献（即连续时间步的差值）

目前，我使用以下顺序：

1. 要将每个时间步提取到一个新文件中，我使用cdo seltimestep,$i ...
2. 使用cdo sub $i ${i-1} ...将每个差异计算到一个新文件中
3. 并最终将这些新文件与cdo mergetime ...合并到一个结果文件中

在我看来，这似乎是非常繁琐的和不理想的表现。由于时间步长太多，我无法使用cdo管道，因此需要同时创建许多文件

有没有更好的解决方案，可以用cdo（或其他类似nco/ncl的东西）将累积变量转换为时间步长值

Answer 1

如果要使用cdo，不需要所有这些循环和编写大量文件，只需使用函数deltat：

cdo deltat in.nc diff.nc 

与python解决方案一样，这将比您使用的循环快几个数量级，并且具有作为命令行单行程序的优势

或者，更不简洁的是，如果您知道长度，您可以将这两个系列区别开来（我展示了这一点，因为这种技术在其他上下文中可能很有用）：

# calculate number of steps in the file:
nstep=$(cdo -s ntime in.nc)

# do difference between steps 2:n and steps 1:(n-1)
cdo sub -seltimestep,2/$nstep in.nc -seltimestep,1/`expr $nstep - 1` in.nc diff.nc

累积字段上的Postscript 请注意，上面的解决方案和本页上发布的两个python解决方案产生的输出比输入少一个timestep，，即它们丢弃了第一个timestep。在某些情况下，例如，如果有一个模型通量场在预测中累积（似乎是这样），则不希望放弃第一个时间步（因为这是从预测开始时的零到第一步的累积）。在这种情况下，您可以提取第一步并将其插入文件的“前面”，如下所示：

cdo mergetime -seltimestep,1 in.nc diff.nc diff_with_step1.nc 

您还应该确保对python解决方案也这样做

您可以将整个过程作为一个管道（有时管道可能会导致总线错误或seg故障，通常可以使用“-L”选项来执行顺序操作）

cdo mergetime -seltimestep,1 in.nc -deltat in.nc diff_with_step1.nc

如果你有seg故障，试试这个

cdo -L mergetime -seltimestep,1 in.nc -deltat in.nc diff_with_step1.nc

这是为了防止舍入和准确性问题，如果您有压缩数据（即键入NC_SHORT）：

cdo -L -b f32 mergetime -seltimestep,1 in.nc -deltat in.nc diff_with_step1.nc