netcdf - 使用 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...

当我尝试使用 ncecat 3B42_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 的数组工具来生成通用算术数组。