Re: [ExternalEmail] [ferret_users] Problem to read descriptor file in new version of Ferret

[Russ Fiedler <russell.fiedler@xxxxxxxx>]

Hi Arnold,

I found some of your CICE files on Raijin and I can see  why Ferret is getting confused. On opening a file (not a descriptor)


[raf599@raijin1 raf599]$ ferret -nojnl
     NOAA/PMEL TMAP
     FERRET v6.842 
     Linux 2.6.32-279.19.1.el6.x86_64 64-bit - 02/01/13
     12-Mar-14 10:32    

yes? use iceh.0001-01.nc
 *** NOTE: Unspecified or unsupported ordering of axes in variable Tinz
 *** NOTE: The default ordering will be used
 *** NOTE: Unspecified or unsupported ordering of axes in variable Tsnz
 *** NOTE: The default ordering will be used
 *** NOTE: Axis TIME marked as TI
 *** NOTE: but it has different directions in grids. Axes in grids may be inconsistent.

There is a strange grid


yes? show grid GBP7
    GRID GBP7
 name       axis              # pts   start                end
 NI        X                  360 r   1                    360
 NJ        Y                  300 r   1                    300
 NKSNOW    Z                    1 r   1                    1
 NC        T                    5 r   1                    5
 TIME      ENSEMBLE             1mr   16-JAN 12:00         16-JAN 12:00
 normal    F


5 dimensions?

And these variables that you didn't show

TINZ     ice internal temperatures        1:360     1:300     1:1       1:5       1:1       ...
 TSNZ     snow internal temperatures       1:360     1:300     1:1       1:5       1:1       ...


    float Tinz(time, nc, nkice, nj, ni) ;
        Tinz:units = "C" ;
        Tinz:long_name = "ice internal temperatures" ;
        Tinz:coordinates = "TLON TLAT VGRDi NCAT time" ;
        Tinz:cell_measures = "area: tarea" ;
        Tinz:missing_value = 1.e+30f ;
        Tinz:_FillValue = 1.e+30f ;
        Tinz:cell_methods = "time: mean" ;
        Tinz:time_rep = "averaged" ;
    float Tsnz(time, nc, nksnow, nj, ni) ;
        Tsnz:units = "C" ;
        Tsnz:long_name = "snow internal temperatures" ;
        Tsnz:coordinates = "TLON TLAT VGRDs NCAT time" ;
        Tsnz:cell_measures = "area: tarea" ;
        Tsnz:missing_value = 1.e+30f ;
        Tsnz:_FillValue = 1.e+30f ;
        Tsnz:cell_methods = "time: mean" ;
        Tsnz:time_rep = "averaged" ;


So in fact ferret has to put time as the fifth axis.

I did a ncks excluding these 2 variables on a couple files, made a descriptor and Ferret worked just fine


yes? use dummy.des
yes? show data
     currently SET data sets:
    1> ./dummy.des  (default)
 name     title                             I         J         K         L         M         N
 ANGLE    angle grid makes with latitude   1:360     1:300     ...       ...       ...       ...
 ANGLET   angle grid makes with latitude   1:360     1:300     ...       ...       ...       ...
 NCAT     category maximum thickness       1:5       ...       ...       ...       ...       ...
 TLAT     T grid center latitude           1:360     1:300     ...       ...       ...       ...
 TLON     T grid center longitude          1:360     1:300     ...       ...       ...       ...
 TSFC     snow/ice surface temperature     1:360     1:300     ...       1:2       ...       ...
 ULAT     U grid center latitude           1:360     1:300     ...       ...       ...       ...
 ULON     U grid center longitude          1:360     1:300     ...       ...       ...       ...
 VGRDI    vertical ice levels              1:1       ...       ...       ...       ...       ...
 VGRDS    vertical snow levels             1:1       ...       ...       ...       ...       ...
 AICE     ice area  (aggregate)            1:360     1:300     ...       1:2       ...       ...
 CONGEL   congelation ice growth           1:360     1:300     ...       1:2       ...       ...
 DARDG2DT ridge area formation rate        1:360     1:300     ...       1:2       ...       ...
 DIVU     strain rate (divergence)         1:360     1:300     ...       1:2       ...       ...
...


So it looks like you need to move Tinz and Tsnz to another file or somehow encode the dimensions nkice and nksnow with nc or maybe eliminate them since they are just 1  long.


Cheers,
Russ






On 11/03/14 17:11, Arnold.Sullivan@xxxxxxxx wrote:

Hi there,

 

I have problem to read the same descriptor file in new version of Ferret

 

FERRET v6.842

In New version I got this the time axial should be under “L” :

******************************************************************************************

 

    1> ./PIC2C-07.iceh.des  (default)

name     title                             I         J         K         L         M         N

TLON     T grid center longitude          1:360     1:300     ...       ...       ...       ...

TLAT     T grid center latitude           1:360     1:300     ...       ...       ...       ...

ULON     U grid center longitude          1:360     1:300     ...       ...       ...       ...

ULAT     U grid center latitude           1:360     1:300     ...       ...       ...       ...

NCAT     category maximum thickness       ...       ...       ...       1:468     ...       ...

VGRDI    vertical ice levels              ...       ...       1:1       ...       ...       ...

VGRDS    vertical snow levels             ...       ...       1:1       ...       ...       ...

TMASK    ocean grid mask                  1:360     1:300     ...       ...       ...       ...

TAREA    area of T grid cells             1:360     1:300     ...       ...       ...       ...

UAREA    area of U grid cells             1:360     1:300     ...       ...       ...       ...

ANGLE    angle grid makes with latitude   1:360     1:300     ...       ...       ...       ...

ANGLET   angle grid makes with latitude   1:360     1:300     ...       ...       ...       ...

HI       grid cell mean ice thickness     1:360     1:300     ...       ...       1:1       ...

HS       grid cell mean snow thickness    1:360     1:300     ...       ...       1:1       ...

TSFC     snow/ice surface temperature     1:360     1:300     ...       ...       1:1       ...

AICE     ice area  (aggregate)            1:360     1:300     ...       ...       1:1       ...

UVEL     ice velocity (x)                 1:360     1:300     ...       ...       1:1       ...

 

Old version: FERRET v6.51

******************************************************************************************

 

    1> ./PIC2C-07.iceh.des  (default)

name     title                             I         J         K         L

TLON     T grid center longitude          1:360     1:300     ...       ...

TLAT     T grid center latitude           1:360     1:300     ...       ...

ULON     U grid center longitude          1:360     1:300     ...       ...

ULAT     U grid center latitude           1:360     1:300     ...       ...

NCAT     category maximum thickness       1:5       ...       ...       ...

VGRDI    vertical ice levels              1:1       ...       ...       ...

VGRDS    vertical snow levels             1:1       ...       ...       ...

TMASK    ocean grid mask                  1:360     1:300     ...       ...

TAREA    area of T grid cells             1:360     1:300     ...       ...

UAREA    area of U grid cells             1:360     1:300     ...       ...

ANGLE    angle grid makes with latitude   1:360     1:300     ...       ...

ANGLET   angle grid makes with latitude   1:360     1:300     ...       ...

HI       grid cell mean ice thickness     1:360     1:300     ...       1:468

HS       grid cell mean snow thickness    1:360     1:300     ...       1:468

TSFC     snow/ice surface temperature     1:360     1:300     ...       1:468

AICE     ice area  (aggregate)            1:360     1:300     ...       1:468

UVEL     ice velocity (x)                 1:360     1:300     ...       1:468

 

Descriptor file looks like this

******************************************************************************************

 

&FORMAT_RECORD

   D_TYPE             = '  MC'

   D_FORMAT           = '  1A'

/

 

&BACKGROUND_RECORD

   D_TITLE            = 'sea ice model output for CICE'

   D_T0TIME           = '01-JAN-0001 00:00:00'

   D_TIME_UNIT        = 86400

   D_TIME_MODULO      = .FALSE.

   D_CALTYPE          = 'proleptic_gregorian'

/

 

&MESSAGE_RECORD

   D_MESSAGE          = ' '

   D_ALERT_ON_OPEN    = .FALSE.

   D_ALERT_ON_OUTPUT  = .FALSE.

/

 

&EXTRA_RECORD

/

 

&STEPFILE_RECORD

   S_FILENAME         = 'iceh.0001-01.nc'

   S_AUX_SET_NUM      = 0

   S_START            = 15.5

   S_END              = 15.5

   S_DELTA            = 1

   S_NUM_OF_FILES     = 1

   S_REGVARFLAG       = ' '

/

 

&STEPFILE_RECORD

   S_FILENAME         = 'iceh.0001-02.nc'

   S_AUX_SET_NUM      = 0

   S_START            = 45

   S_END              = 45

   S_DELTA            = 1

   S_NUM_OF_FILES     = 1

   S_REGVARFLAG       = ' '

/

 

******************************************************************************************

Any suggestion?

 

Regards,

 

Arnold