使用 NCO 连接每日 TRMM netCDF 文件时如何添加时间维度？

问题描述

我从 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 上的文档中，我无法弄清楚如何实现这两种方式中的任何一种。或者是否有第三种更好的方法来实现这一点（添加相关时间信息的串联）？

解决方法

由于数据文件不遵循 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 的数组工具用于通用算术数组。