[netcdf-java] Using GridDataset to read GSD HRRR file

I would like to read a GSD generated HRRR file in netcdf format as a GridDataset. I can open the file as a GridDataset, but when I try to ask for the grids, it says there are no grids. Using getDetailInfo(), it looks like it builds some coordinate systems for each of the variables in the file, but then when trying to get them, they don't appear to exist. What makes a file a grid?

Below is the output from getDetailInfo() and ncdump on the GSD HRRR netcdf file.

Here is the information from getDetailInfo()
2011-01-06 15:15:04,047 [main] INFO org.wdssii.ncingest.NDFDIngest - Opening /mnt/home/jbrogden/hrrr1027416000306.nc for reading

GeoReferencing Coordinate Axes
Name__________________________Units_______________Type______Description
x(west_east=1799) km GeoX synthesized GeoX coordinate from DX attribute y(south_north=1059) km GeoY synthesized GeoY coordinate from DY attribute Time(Time=1) secs since 1970-01-01 00:00:00Time synthesized time coordinate from Times(time)


----------------------------------------------------
  made CoordinateTransformVariable: Lambert
Parsing with Convention = WRF
 Coordinate Axis added = x type= GeoX
 Coordinate Variable Alias added = x for dimension= west_east
 Coordinate Axis added = y type= GeoY
 Coordinate Variable Alias added = y for dimension= south_north
Coordinate Variable added = Time for dimension Time = UNLIMITED; // (1 currently)
 Coordinate Axis added = Time type= Time
 created maximal CoordSystem 'Time y x' for var= VIL
 assigned maximal CoordSystem 'Time y x' for var= RADARVIL
 assigned maximal CoordSystem 'Time y x' for var= ECHOTOP
 assigned maximal CoordSystem 'Time y x' for var= T2M
 assigned maximal CoordSystem 'Time y x' for var= TD2M
 assigned maximal CoordSystem 'Time y x' for var= U10M
 assigned maximal CoordSystem 'Time y x' for var= V10M
 assigned maximal CoordSystem 'Time y x' for var= P_SFC
***assign (implicit coordAxisType) coordTransform Lambert to CoordSys= Time y x


----------------------------------------------------
netcdf /mnt/home/jbrogden/hrrr1027416000306.nc {
 dimensions:
   Time = UNLIMITED;   // (1 currently)
   DateStrLen = 19;
   west_east = 1799;
   south_north = 1059;
   num_snapshots = 13;
 variables:
   char Times(Time=1, DateStrLen=19);
   float VIL(Time=1, num_snapshots=13, south_north=1059, west_east=1799);
     :FieldType = 104; // int
     :MemoryOrder = "XYZ";
     :description = "VERTICALLY INTEGRATED LIQUID WATER";
     :units = "kg m-2";
     :stagger = "";
     :coordinates = "XLONG XLAT";
float RADARVIL(Time=1, num_snapshots=13, south_north=1059, west_east=1799);
     :FieldType = 104; // int
     :MemoryOrder = "XYZ";
     :description = "VERTICALLY INTEGRATED LIQUID WATER FROM Ze";
     :units = "kg m-2";
     :stagger = "";
     :coordinates = "XLONG XLAT";
float ECHOTOP(Time=1, num_snapshots=13, south_north=1059, west_east=1799);
     :FieldType = 104; // int
     :MemoryOrder = "XYZ";
     :description = "echo top height FROM Ze";
     :units = "m";
     :stagger = "";
     :coordinates = "XLONG XLAT";
   float T2M(Time=1, num_snapshots=13, south_north=1059, west_east=1799);
     :FieldType = 104; // int
     :MemoryOrder = "XYZ";
     :description = "2-m temperature";
     :units = "K";
     :stagger = "";
     :coordinates = "XLONG XLAT";
   float TD2M(Time=1, num_snapshots=13, south_north=1059, west_east=1799);
     :FieldType = 104; // int
     :MemoryOrder = "XYZ";
     :description = "2-m dew point";
     :units = "K";
     :stagger = "";
     :coordinates = "XLONG XLAT";
   float U10M(Time=1, num_snapshots=13, south_north=1059, west_east=1799);
     :FieldType = 104; // int
     :MemoryOrder = "XYZ";
     :description = "10-m u-component of wind";
     :units = "m/s";
     :stagger = "";
     :coordinates = "XLONG XLAT";
   float V10M(Time=1, num_snapshots=13, south_north=1059, west_east=1799);
     :FieldType = 104; // int
     :MemoryOrder = "XYZ";
     :description = "10-m v-component of wind";
     :units = "m/s";
     :stagger = "";
     :coordinates = "XLONG XLAT";
   float P_SFC(Time=1, num_snapshots=13, south_north=1059, west_east=1799);
     :FieldType = 104; // int
     :MemoryOrder = "XYZ";
     :description = "surface pressure";
     :units = "Pa";
     :stagger = "";
     :coordinates = "XLONG XLAT";
   char Lambert;
     :grid_mapping_name = "lambert_conformal_conic";
     :latitude_of_projection_origin = 38.500003814697266; // double
     :longitude_of_central_meridian = -97.5; // double
     :standard_parallel = 38.5; // double
     :_CoordinateTransformType = "Projection";
     :_CoordinateAxisTypes = "GeoX GeoY";
   double x(west_east=1799);
     :units = "km";
     :long_name = "synthesized GeoX coordinate from DX attribute";
     :_CoordinateAxisType = "GeoX";
     :_CoordinateAliasForDimension = "west_east";
   double y(south_north=1059);
     :units = "km";
     :long_name = "synthesized GeoY coordinate from DY attribute";
     :_CoordinateAxisType = "GeoY";
     :_CoordinateAliasForDimension = "south_north";
   double Time(Time=1);
     :units = "secs since 1970-01-01 00:00:00";
     :long_name = "synthesized time coordinate from Times(time)";
     :_CoordinateAxisType = "Time";


Here is the header from ncdump.
ncdump netcdf ncdump {
dimensions:
        Time = UNLIMITED ; // (1 currently)
        DateStrLen = 19 ;
        west_east = 1799 ;
        south_north = 1059 ;
        num_snapshots = 13 ;
variables:
        char Times(Time, DateStrLen) ;
        float VIL(Time, num_snapshots, south_north, west_east) ;
                VIL:FieldType = 104 ;
                VIL:MemoryOrder = "XYZ" ;
                VIL:description = "VERTICALLY INTEGRATED LIQUID WATER" ;
                VIL:units = "kg m-2" ;
                VIL:stagger = "" ;
                VIL:coordinates = "XLONG XLAT" ;
        float RADARVIL(Time, num_snapshots, south_north, west_east) ;
                RADARVIL:FieldType = 104 ;
                RADARVIL:MemoryOrder = "XYZ" ;
RADARVIL:description = "VERTICALLY INTEGRATED LIQUID WATER FROM Ze" ;
                RADARVIL:units = "kg m-2" ;
                RADARVIL:stagger = "" ;
                RADARVIL:coordinates = "XLONG XLAT" ;
        float ECHOTOP(Time, num_snapshots, south_north, west_east) ;
                ECHOTOP:FieldType = 104 ;
                ECHOTOP:MemoryOrder = "XYZ" ;
                ECHOTOP:description = "echo top height FROM Ze" ;
                ECHOTOP:units = "m" ;
                ECHOTOP:stagger = "" ;
                ECHOTOP:coordinates = "XLONG XLAT" ;
        float T2M(Time, num_snapshots, south_north, west_east) ;
                T2M:FieldType = 104 ;
                T2M:MemoryOrder = "XYZ" ;
                T2M:description = "2-m temperature" ;
                T2M:units = "K" ;
                T2M:stagger = "" ;
                T2M:coordinates = "XLONG XLAT" ;
        float TD2M(Time, num_snapshots, south_north, west_east) ;
                TD2M:FieldType = 104 ;
                TD2M:MemoryOrder = "XYZ" ;
                TD2M:description = "2-m dew point" ;
                TD2M:units = "K" ;
                TD2M:stagger = "" ;
                TD2M:coordinates = "XLONG XLAT" ;
        float U10M(Time, num_snapshots, south_north, west_east) ;
                U10M:FieldType = 104 ;
                U10M:MemoryOrder = "XYZ" ;
                U10M:description = "10-m u-component of wind" ;
                U10M:units = "m/s" ;
                U10M:stagger = "" ;
                U10M:coordinates = "XLONG XLAT" ;
        float V10M(Time, num_snapshots, south_north, west_east) ;
                V10M:FieldType = 104 ;
                V10M:MemoryOrder = "XYZ" ;
                V10M:description = "10-m v-component of wind" ;
                V10M:units = "m/s" ;
                V10M:stagger = "" ;
                V10M:coordinates = "XLONG XLAT" ;
        float P_SFC(Time, num_snapshots, south_north, west_east) ;
                P_SFC:FieldType = 104 ;
                P_SFC:MemoryOrder = "XYZ" ;
                P_SFC:description = "surface pressure" ;
                P_SFC:units = "Pa" ;
                P_SFC:stagger = "" ;
                P_SFC:coordinates = "XLONG XLAT" ;

// global attributes:
                :TITLE = " OUTPUT FROM WRF V3.2 MODEL" ;
                :START_DATE = "2010-10-01_16:00:00" ;
                :WEST-EAST_GRID_DIMENSION = 1800 ;
                :SOUTH-NORTH_GRID_DIMENSION = 1060 ;
                :BOTTOM-TOP_GRID_DIMENSION = 51 ;
                :DX = 3000.f ;
                :DY = 3000.f ;
                :GRIDTYPE = "C" ;
                :DIFF_OPT = 1 ;
                :KM_OPT = 4 ;
                :DAMP_OPT = 1 ;
                :DAMPCOEF = 0.02f ;
                :KHDIF = 0.f ;
                :KVDIF = 0.f ;
                :MP_PHYSICS = 8 ;
                :RA_LW_PHYSICS = 1 ;
                :RA_SW_PHYSICS = 1 ;
                :SF_SFCLAY_PHYSICS = 2 ;
                :SF_SURFACE_PHYSICS = 3 ;
                :BL_PBL_PHYSICS = 2 ;
                :CU_PHYSICS = 0 ;
                :SURFACE_INPUT_SOURCE = 1 ;
                :SST_UPDATE = 0 ;
                :GRID_FDDA = 0 ;
                :GFDDA_INTERVAL_M = 0 ;
                :GFDDA_END_H = 0 ;
                :GRID_SFDDA = 0 ;
                :SGFDDA_INTERVAL_M = 0 ;
                :SGFDDA_END_H = 0 ;
                :WEST-EAST_PATCH_START_UNSTAG = 1 ;
                :WEST-EAST_PATCH_END_UNSTAG = 1799 ;
                :WEST-EAST_PATCH_START_STAG = 1 ;
                :WEST-EAST_PATCH_END_STAG = 1800 ;
                :SOUTH-NORTH_PATCH_START_UNSTAG = 1 ;
                :SOUTH-NORTH_PATCH_END_UNSTAG = 1059 ;
                :SOUTH-NORTH_PATCH_START_STAG = 1 ;
                :SOUTH-NORTH_PATCH_END_STAG = 1060 ;
                :BOTTOM-TOP_PATCH_START_UNSTAG = 1 ;
                :BOTTOM-TOP_PATCH_END_UNSTAG = 50 ;
                :BOTTOM-TOP_PATCH_START_STAG = 1 ;
                :BOTTOM-TOP_PATCH_END_STAG = 51 ;
                :GRID_ID = 1 ;
                :PARENT_ID = 1 ;
                :I_PARENT_START = 0 ;
                :J_PARENT_START = 0 ;
                :PARENT_GRID_RATIO = 1 ;
                :DT = 13.44f ;
                :CEN_LAT = 38.5f ;
                :CEN_LON = -97.5f ;
                :TRUELAT1 = 38.5f ;
                :TRUELAT2 = 38.5f ;
                :MOAD_CEN_LAT = 38.5f ;
                :STAND_LON = -97.5f ;
                :POLE_LAT = 90.f ;
                :POLE_LON = 0.f ;
                :GMT = 16.f ;
                :JULYR = 2010 ;
                :JULDAY = 274 ;
                :MAP_PROJ = 1 ;
                :MMINLU = "USGS" ;
                :NUM_LAND_CAT = 24 ;
                :ISWATER = 16 ;
                :ISLAKE = -1 ;
                :ISICE = 24 ;
                :ISURBAN = 1 ;
                :ISOILWATER = 14 ;

--
=======================================================================
Jeff Brogden    CIMMS/OU
Jeff.Brogden@xxxxxxxx



  • 2011 messages navigation, sorted by:
    1. Thread
    2. Subject
    3. Author
    4. Date
    5. ↑ Table Of Contents
  • Search the netcdf-java archives: