使用 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 B42_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 时，我遇到了同样的问题 - 同样，我找不到如何在连接时使用全局属性中包含的时间信息数据。

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

首先使用一些命令“手动”为每个文件添加时间戳，然后使用 ncecat
不知何故 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 的数组工具用于通用算术数组。

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

How to add time dimension when concatenating daily TRMM netCDF files using NCO?

netcdf

nco

netcdf4