Question

我在nco工作并遇到了这个难题。我必须将360天的日历重新安排到365天的正常格里高利日历中。在360天的日历中，每个月有30天。我的想法是在格里高利日历中添加一天到一个月应该包含31天，并从2月（闰年或非闰年）中删除一天或两天。我知道如何删除天数（只需将你需要的东西复制到另一个文件中）。

是否有一种简单的方法可以复制特定月份的最后一天并将其添加到同一个月？

我的文件分辨率为3小时，我希望使用名为tas（温度）的变量。下面我有最重要的部分：

dimensions:
rlon = 424 ;
rlat = 412 ;
time = UNLIMITED ; // (2880 currently)

variables:
double rlon(rlon) ;
    rlon:standard_name = "grid_longitude" ;
    rlon:long_name = "longitude in rotated pole grid" ;
    rlon:units = "degrees" ;
    rlon:axis = "X" ;
double rlat(rlat) ;
    rlat:standard_name = "grid_latitude" ;
    rlat:long_name = "latitude in rotated pole grid" ;
    rlat:units = "degrees" ;
double time(time) ;
    time:standard_name = "time" ;
    time:units = "days since 1949-12-01 00:00:00" ;
    time:calendar = "360_day" ;
    time:long_name = "time" ;
    time:axis = "T" ;
float tas(time, rlat, rlon) ;
    tas:grid_mapping = "rotated_pole" ;
    tas:_FillValue = 1.e+20f ;
    tas:missing_value = 1.e+20f ;
    tas:standard_name = "air_temperature" ;
    tas:long_name = "Near-Surface Air Temperature" ;
    tas:units = "K" ;
    tas:coordinates = "lon lat height" ;
    tas:cell_methods = "time: point" ;

我期待听到想法。

Answer 1

我用更长的方法解决了这个问题。首先，我将年度文件分成几个月。我每个月都会添加并删除必要的日期，然后连接几个月。

如果有人有更快的方法请分享

Answer 2

如果360天数据存储在12个月度文件中，那么您可以使用ncrcat记录附加feature以及负的hyperslab索引从末尾向后计数并将最后一条记录追加到自身

ncrcat --rec_apn -d time,-1 in.nc in.nc

请务必先备份数据，因为这会覆盖输入。

复制月份的最后一天，并使用nco将其添加到netcdf文件的末尾

2 个答案: