NetCDF 空间合并到全局数据

Question

目前我使用全局降水量(ppt)和潜在蒸散量(pet)数据来计算SPEI。由于我的硬件资源有限，我将全局ppt和ppt数据分成了32个部分，每个文件覆盖45x45deg，包含756个数据——从1958-2020每个月(tile01.nc, tile02.nc, ... tile32.nc)

• 例如，我使用 cdo sellonlatbox,-180,-135,45,90 in.nc out.nc 或 ncks -d lat,45.,90。 -d lon,-180.,-135. in.nc -O out.nc
• 根据 SPEI 脚本的要求，我使用 ncpdq 将维度从 time,lat,lon 重新排序并固定为 lat,lon,time和 ncks。
• 从 SPEI 输出中，我得到了 lat,lon,time 的输出。所以我使用 ncpdq 重新排序维度，使其成为 time,lat,lon。
• 每个图 block SPEI 输出覆盖 45x45deg 并包含 756 个 SPEI 数据 - 从 1958 年到 2020 年每月

最后，我需要将所有输出(32 个文件)合并为一个，这样我将获得全局 SPEI 输出。我尝试使用 cdo mergegrid 但结果不是我所期望的。 cdo 或 nco 是否有任何命令来解决这个问题，如果我们正在处理 geoTIFF 格式，其功能类似于 gdal_merge ？

下面是SPEI输出的例子

更新

我按照下面 Robert 的建议设法使用 cdo collgrid 合并了所有数据。结果如下:

Answer 1

我相信您想为此使用 CDO 的 distgrid 和 collgrid 方法。

首先，运行这个:

cdo distgrid,4,8 in.nc obase

这将按照您希望的方式拆分文件。

然后对文件进行必要的后处理。

然后使用collgrid合并文件:

cdo collgrid obase* out.nc

根据您拆分文件的方式，您可以在当前的工作流程中使用 collgrid 代替 mergegrid。