NetCDF分块空间合并为全球数据：如何合并32块45°×45°SPEI输出瓦片

我目前正在使用全球降水（ppt）和潜在蒸散发（pet）数据计算SPEI。受硬件资源限制，我先将全球ppt与pet数据切分为32个分块，每个文件覆盖45°×45°范围，包含1958-2020年共756个月度数据，文件命名为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和ncks将维度顺序从time,lat,lon重排调整为lat,lon,time。
得到的SPEI输出维度为lat,lon,time，因此我又使用ncpdq将维度重排回time,lat,lon。

每个分块的SPEI输出都覆盖45°×45°范围，包含1958-2020年共756个月度SPEI数据。

我最终需要将32个输出文件合并为一个全球SPEI输出文件，尝试使用cdo mergegrid但结果不符合预期，请问cdo或nco中是否有类似处理geoTIFF时gdal_merge功能的命令可以解决这个问题？

SPEI输出示例：

问题解决更新 #

我已经按照建议使用cdo collgrid完成了所有数据的合并，合并结果如下：

内容的提问来源于stack exchange，提问作者user97103