如何在使用NCO串联每日TRMM netCDF文件时添加时间维度？

Question

我从 https://disc.gsfc.nasa.gov/datasets 下载了几天的每日TRMM 3B42数据。文件名的格式为 3B42_Daily.yyyymmdd.7.nc4.nc4，但文件不包含任何时间维度。因此，当我使用 ncecat 连接各种文件时，结果文件中缺少日期信息。我想知道如何在组合数据集中添加时间信息。

时间戳似乎是全局属性的一部分。以下是来自 ncdump 的相关部分：

$ ncdump -h ~/Downloads/3B42_Daily.19980730.7.nc4.nc4 
netcdf \3B42_Daily.19980730.7.nc4 {
dimensions:
    lon = 201 ;
    lat = 201 ;
variables:
    float lon(lon) ;
...trimmed...
// global attributes:
        :BeginDate = "1998-07-30" ;
        :BeginTime = "01:30:00.000Z" ;
        :EndDate = "1998-07-31" ;
        :EndTime = "01:29:59.999Z" ;
...trimmed...

当我尝试使用ncecat3B42_Daily.199808？？.7.nc4.nc4 /tmp/daily.nc4它给出了

$ ncdump -h /tmp/daily.nc
netcdf daily {
dimensions:
    record = UNLIMITED ; // (5 currently)
    lon = 201 ;
    lat = 201 ;
variables:
    float lon(lon) ;
...trimmed...
// global attributes:
        :BeginDate = "1998-08-01" ;
        :BeginTime = "01:30:00.000Z" ;
        :EndDate = "1998-08-02" ;
        :EndTime = "01:29:59.999Z" ;
...trimmed...

保留第一个文件的全局属性中的时间信息，这不是很有用。

当我试图在python中使用xarray时，我面临着同样的问题——再次，我无法找到在连接数据时如何使用全局属性中包含的时间信息。

我能想到两个可能的解决方案，但我不知道如何实现它们。

1. 首先使用一些命令“手动”向每个文件添加时间戳，然后使用ncecat
2. 以某种方式ncecat可以读取全局属性，并在连接时将其转换为维度和变量

从http://nco.sourceforge.net/nco.html的留档，我不知道如何实现这两种方式中的任何一种。或者有第三种更好的方法来实现这一点（添加相关时间信息的串联）？

Answer 1

由于数据文件不遵循 CF 约定，因此在使用 ncecat 连接文件后，您可能必须手动创建时间坐标。如果时间是规则的，则只需要几个命令，例如，

ncecat -u time in*.nc out.nc
ncap2 -O -s 'time[time]={0.0,1.0,2.0,3.0,4.0}' out.nc out.nc
ncatted -a units,time,o,c,"days since 1998-08-01" out.nc
    

或者将ncap2的数组工具用于通用算术数组。