在记录维度 (nco) 上连接 ncdf 文件和块

问题描述

我正在尝试将大量（~4000）个 ncdf 文件连接到一个文件中。每个输入文件都是一个空间栅格，具有 x 和 y 维度。

我正在尝试与 ncecat 合作：

ncecat -4 -L 5 -D 2 --open_ram --cnk_csh=1000000000 \
   --cnk_dmn record,2000 --cnk_dmn x,10 --cnk_dmn y,10 \
   $input_files output.nc

这给了我这样的东西：

netcdf 测试 { 维度： 记录=无限； //（当前 6 个） y = 11250 ; x = 15000 ;变量： 浮动带1（记录，y，x）； Band1:long_name = "GDAL 频段编号 1" ; Band1:_FillValue = -3.4e+38f ; Band1:grid_mapping = "transverse_mercator" ; Band1:_Storage = "chunked" ; Band1:_ChunkSizes = 1,10,10 ; Band1:_DeflateLevel = 5 Band1:_Filter = "|1,5; Band1:_Shuffle = "true" ; Band1:_Endianness = "小";

，并且记录维度实际上没有分块。 我想我可以先运行这个命令，然后在输出文件上使用 ncks 来修复记录变暗并再次重新分块，但是，由于 ncks 需要将所有内容读入 ram，而且也是另一个耗时的操作，我正在寻找一个告诉 ncecat 它还应该将 record-dim 视为分块维度的方式。我还没有找到办法做到这一点。

解决方法

您的命令看起来格式正确，但我会提出一些意见。首先，您看到的行为可能是一个错误，因为该命令应根据请求生成大小为 2000 的记录维度块。其次，请阅读分块文档 here。这导致添加 --cnk_plc=cnk_xpl 选项可能会有所帮助。第三，我建议您使用 ncrcat 而不是 ncecat 来连接和分块文件。如here 所述，前者比后者占用内存更少。