Извлеките переменную NetCDF и создайте новую NetCDF.

Question

Мне нужна помощь в работе с файлами NetCDF. Всего у меня 10 файлов за 10 лет соответственно. Каждый год имеет несколько (одинаковых) переменных, некоторые из них также охватывают дневные значения. Здесь я покажу вам один пример структуры:

      (base) thess2ice@local:rhone_smb_modelling $ ncdump -h Rhone_AWS_1990.nc
netcdf Rhone_AWS_1990 {
dimensions:
    x = 402 ;
    y = 852 ;
    time = 1460 ;
variables:
    double x(x) ;
        x:_FillValue = NaN ;
        x:standard_name = "x" ;
        x:long_name = "longitude" ;
        x:units = "degrees_east" ;
    double y(y) ;
        y:_FillValue = NaN ;
        y:standard_name = "y" ;
        y:long_name = "latitude" ;
        y:units = "degrees_north" ;
    float HGT(y, x) ;
        HGT:_FillValue = -9999.f ;
        HGT:units = "m" ;
        HGT:long_name = "Elevation" ;
    float ASPECT(y, x) ;
        ASPECT:_FillValue = -9999.f ;
        ASPECT:units = "degrees" ;
        ASPECT:long_name = "Aspect of slope" ;
    float SLOPE(y, x) ;
        SLOPE:_FillValue = -9999.f ;
        SLOPE:units = "degrees" ;
        SLOPE:long_name = "Terrain slope" ;
    float MASK(y, x) ;
        MASK:_FillValue = -9999.f ;
        MASK:units = "boolean" ;
        MASK:long_name = "Glacier mask" ;
    int64 time(time) ;
        time:units = "hours since 1990-01-01 00:00:00" ;
        time:calendar = "proleptic_gregorian" ;
    double T2(time, y, x) ;
        T2:_FillValue = NaN ;
        T2:units = "K" ;
        T2:long_name = "Temperature at 2 m" ;
    double RRR(time, y, x) ;
        RRR:_FillValue = NaN ;
        RRR:units = "mm" ;
        RRR:long_name = "Total precipitation (liquid+solid)" ;
    double ACC(y, x) ;
        ACC:_FillValue = -9999. ;
        ACC:units = "mm yr^-1" ;
        ACC:long_name = "Accumulation from RRR_solid" ;
    double MELT_I(y, x) ;
        MELT_I:_FillValue = -9999. ;
        MELT_I:units = "mm yr^-1" ;
        MELT_I:long_name = "Melt from PDD" ;
    double MELT_S(y, x) ;
        MELT_S:_FillValue = -9999. ;
        MELT_S:units = "mm yr^-1" ;
        MELT_S:long_name = "Melt from PDD" ;
    double SMB(y, x) ;
        SMB:_FillValue = -9999. ;
        SMB:units = "mm yr^-1" ;
        SMB:long_name = "SMB from PDD" ;

Мне нужны данные, обрабатываемые в качестве входных данных для модели. Переменная, которую мне нужно извлечь из каждого из 10 файлов NetCDF, — это переменная SMB, которая является только годовым значением для каждой ячейки сетки. Итак, я хотел бы построить NetCDF в форме:

      (year, y, x) for the SMB variable

Я уже знаю команду ncks для извлечения только переменной SMB, но мне не удается применить ее сразу к нескольким файлам (скажем, ко всем файлам nc в текущем каталоге) и впоследствии свести их в один файл NetCDF, охватывающий все 10 лет.

Может ли кто-нибудь помочь мне с этим? Было бы замечательно!

Есть

Answer 1

унтер-офицер ncecatкоманда, задокументированная здесь, делает именно то, что вам кажется нужным:

      ncecat -u time -v SMB Rhone_AWS_199*.nc out.nc