Q: how to create a multi-file netCDF data set

To: ferret_users@xxxxxxxxxxxxxxxxxxx
Subject: Q: how to create a multi-file netCDF data set
From: Bob Yu <yuy@xxxxxxxxxxxxxxxxxx>
Date: Thu, 08 May 2003 10:24:44 -0400
References: <5.1.0.14.0.20030506155858.00ad8e90@mail.flinders.edu.au><3EB7EE54.364EE5BA@noaa.gov> <3EB95189.24D4A3F3@nemo.gsfc.nasa.gov>
Sender: owner-ferret_users@xxxxxxxxxxxxxxxxxxxx

Dear ferreters,

As a new user of ferret, I currently need to build up a dataset from multiple netCDF files. The dataset will have exact data structure as in the individual files, except the T dimension elements. individual file has only one element in T, while the built dataset will has, say 4 elements from the 4 files. I created a descriptor file and "linuxizered" it into F90 format, as is indicated by the user's guide. I got the following error. Could you please tell me what is wrong in the BACKGROUNG_RECORD?

        NOAA/PMEL TMAP
        FERRET v5.50
        Linux 2.4.3-12smp - 01/15/03
        08-May-03 10:00
yes? set data y.des
**TMAP ERR: Error in namelist record
             Error in: BACKGROUND_RECORD
             Data set: ./y.des

(NOTE: I did have put a space before "$" and have two spaces in ' MC' ' 1A', as Ansley pointed out.)

------------ Here is the y.des linuxizered from yy.des (below) ----------

&FORMAT_RECORD
        D_TYPE = ' MC',
        D_FORMAT = ' 1A'
/
&BACKGROUND_RECORD
        D_TITLE = 'Model output HYCOM expt 01.9',
        D_T0TIME = '16-JAN-0030 00:00:00',
        D_TIME_UNIT = 3600.0,
        D_CALTYPE = 360_DAY
/
&MESSAGE_RECORD
/
&EXTRA_RECORD
/
&STEPFILE_RECORD
        s_filename     = 'archv.0030_016_00_mix.nc' ,
/
&STEPFILE_RECORD
        s_filename     = 'archv.0030_046_00_mix.nc' ,
/
&STEPFILE_RECORD
        s_filename     = 'archv.0030_076_00_mix.nc' ,
/
&STEPFILE_RECORD
        s_filename     = 'archv.0030_106_00_mix.nc' ,
/
&STEPFILE_RECORD
        s_filename     = 'archv.0030_136_00_mix.nc' ,
/
&STEPFILE_RECORD
        S_FILENAME = '**END OF STEPFILES**'

------------------ Here is the yy.des ------------------

**********************************************
* my descriptor test
**********************************************
$FORMAT_RECORD
        D_TYPE = ' MC',
        D_FORMAT = ' 1A'
$END
$BACKGROUND_RECORD
        D_TITLE = 'Model output HYCOM expt 01.9',
        D_T0TIME = '16-JAN-0030 00:00:00',
        D_TIME_UNIT = 3600.0,
        D_CALTYPE = 360_DAY
$END
$MESSAGE_RECORD
$END
$EXTRA_RECORD
$END
$STEPFILE_RECORD
        s_filename     = 'archv.0030_016_00_mix.nc' ,
$END
$STEPFILE_RECORD
        s_filename     = 'archv.0030_046_00_mix.nc' ,
$END
$STEPFILE_RECORD
        s_filename     = 'archv.0030_076_00_mix.nc' ,
$END
$STEPFILE_RECORD
        s_filename     = 'archv.0030_106_00_mix.nc' ,
$END
$STEPFILE_RECORD
        s_filename     = 'archv.0030_136_00_mix.nc' ,
$END
**************************************************
$STEPFILE_RECORD
        S_FILENAME = '**END OF STEPFILES**'
$END
**************************************************

------------------Here is the cdl for individual file -----------------
netcdf archv.0030_016_00_mix {
dimensions:
        MT = UNLIMITED ; // (1 currently)
        Latitude = 67 ;
        Longitude = 180 ;
variables:
        double MT(MT) ;
                MT:long_name = "model time" ;
                MT:units = "days since 0001-01-16 00:00:00" ;
                MT:calendar = "360_day" ;
        double Date(MT) ;
                Date:long_name = "date" ;
                Date:units = "day as %Y%m%d.%f" ;
        float Latitude(Latitude) ;
                Latitude:units = "degrees_north" ;
        float Longitude(Longitude) ;
                Longitude:units = "degrees_east" ;
                Longitude:point_spacing = "even" ;
                Longitude:modulo = "360 degrees" ;
        float mixed_layer_u_velocity(MT, Latitude, Longitude) ;
                mixed_layer_u_velocity:coordinates = "Date" ;
                mixed_layer_u_velocity:_FillValue = 1.267651e+30f ;
                mixed_layer_u_velocity:valid_range = -104.691f, 72.17733f ;
                mixed_layer_u_velocity:units = "cm/s" ;
                mixed_layer_u_velocity:long_name = "mix.l. u-velocity [01.9H]" ;
        float mixed_layer_v_velocity(MT, Latitude, Longitude) ;
                mixed_layer_v_velocity:coordinates = "Date" ;
                mixed_layer_v_velocity:_FillValue = 1.267651e+30f ;
                mixed_layer_v_velocity:valid_range = -82.96401f, 81.37826f ;
                mixed_layer_v_velocity:units = "cm/s" ;
                mixed_layer_v_velocity:long_name = " mixl. v-velocity [01.9H]" ;
        float surface_boundary_layer_thickness(MT, Latitude, Longitude) ;
                surface_boundary_layer_thickness:coordinates = "Date" ;
                surface_boundary_layer_thickness:_FillValue = 1.267651e+30f ;
                surface_boundary_layer_thickness:valid_range = 2.891365f, 623.0751f ;
                surface_boundary_layer_thickness:units = "m" ;
                surface_boundary_layer_thickness:long_name = "bnd.layr.thickness [01.9H]" ;
        float mixed_layer_thickness(MT, Latitude, Longitude) ;
                mixed_layer_thickness:coordinates = "Date" ;
                mixed_layer_thickness:_FillValue = 1.267651e+30f ;
                mixed_layer_thickness:valid_range = 2.881365f, 1150.608f ;
                mixed_layer_thickness:units = "m" ;
                mixed_layer_thickness:long_name = "mix.layr.thickness [01.9H]" ;
        float mixed_layer_density(MT, Latitude, Longitude) ;
                mixed_layer_density:coordinates = "Date" ;
                mixed_layer_density:_FillValue = 1.267651e+30f ;
                mixed_layer_density:valid_range = 14.88034f, 28.1508f ;
                mixed_layer_density:units = "sigma" ;
                mixed_layer_density:long_name = "mix.layr.dens      [01.9H]" ;

// global attributes:
                :Conventions = "CF-1.0" ;
                :title = "HYCOM GLBa2.00" ;
                :institution = "Code 971, NASA/GSFC" ;
                :source = "HYCOM archive file" ;
                :experiment = "01.9" ;
                :history = "archv2ncdf2d" ;
data:

MT = 10440 ;

Latitude = -40.22064, -38.67636, -37.09803, -35.48612, -33.84122, -32.16404,
    -30.45541, -28.71628, -26.94775, -25.15103, -23.32746, -21.47852,
    -19.60579, -17.71101, -15.79601, -13.86273, -11.91322, -9.949614,
    -7.974132, -5.989064, -3.996755, -1.999594, -2.924075e-10, 1.999594,
    3.996755, 5.989064, 7.974132, 9.949614, 11.91322, 13.86273, 15.79601,
    17.71101, 19.60579, 21.47852, 23.32746, 25.15103, 26.94775, 28.71628,
    30.45541, 32.16404, 33.84122, 35.48612, 37.09803, 38.67636, 40.22064,
    41.73051, 43.20571, 44.64608, 46.05155, 47.42214, 48.75796, 50.05918,
    51.32603, 52.55884, 53.75795, 54.92377, 56.05677, 57.15743, 58.22628,
    59.26389, 60.27082, 61.24769, 62.19513, 63.11375, 64.00423, 64.8672,
    65.70332 ;

Longitude = 34, 36, 38, 40, 42, 44, 46, 48, 50, 52, 54, 56, 58, 60, 62, 64,
    66, 68, 70, 72, 74, 76, 78, 80, 82, 84, 86, 88, 90, 92, 94, 96, 98, 100,
    102, 104, 106, 108, 110, 112, 114, 116, 118, 120, 122, 124, 126, 128,
    130, 132, 134, 136, 138, 140, 142, 144, 146, 148, 150, 152, 154, 156,
    158, 160, 162, 164, 166, 168, 170, 172, 174, 176, 178, 180, 182, 184,
    186, 188, 190, 192, 194, 196, 198, 200, 202, 204, 206, 208, 210, 212,
    214, 216, 218, 220, 222, 224, 226, 228, 230, 232, 234, 236, 238, 240,
    242, 244, 246, 248, 250, 252, 254, 256, 258, 260, 262, 264, 266, 268,
    270, 272, 274, 276, 278, 280, 282, 284, 286, 288, 290, 292, 294, 296,
    298, 300, 302, 304, 306, 308, 310, 312, 314, 316, 318, 320, 322, 324,
    326, 328, 330, 332, 334, 336, 338, 340, 342, 344, 346, 348, 350, 352,
    354, 356, 358, 360, 362, 364, 366, 368, 370, 372, 374, 376, 378, 380,
    382, 384, 386, 388, 390, 392 ;
}

I need you guys help, please!!

-- 
(Bob) Yunyue YU
Code 971, Goddard Space Flight Center
Greenbelt, MD 20771
(301)614-6850  yuy@nemo.gsfc.nasa.gov

References:
- Concatting cdf files
  - From: Hemerson Tonin
- Re: Concatting cdf files
  - From: Ansley Manke
- creating a multi-file netCDF data set
  - From: Bob Yu

Previous by thread: Re: creating a multi-file netCDF data set
Next by thread: Test -- please ignore

[Thread Prev][Thread Next][Index]

Dept of Commerce / NOAA / OAR / PMEL / TMAP