[Met_help] [rt.rap.ucar.edu #90346] History for unable to read 2D field

[John Halley Gotway via RT]

----------------------------------------------------------------
  Initial Request
----------------------------------------------------------------

Here is the 2D field file:

/scratch3/BMC/wrf-chem/KaYee/project_smoke/standalone_upp_met/Testing_HRRR_SmokeV5/test_met
/WRFTWO*
The var we need is "COLMD" for WRFTWO

The 3D field data is WRFNAT* and the var we need is "MASSDEN"
/scratch3/BMC/wrf-chem/KaYee/project_smoke/standalone_upp_met/Test_new/WRFNAT.Grb*

Ka Yee

----------------------------------------------------------------
  Complete Ticket History
----------------------------------------------------------------

Subject: unable to read 2D field
From: John Halley Gotway
Time: Fri May 24 10:14:33 2019

Ka Yee,

I moved this question over to met-help.  Thanks for pointing me to
your data on theia.  I'm working in this directory:

/scratch4/BMC/dtc/John.H.Gotway/MET/MET_Help/wong_data_20190524

I started by loading the MET module:

module use /contrib/modulefiles
module load met/8.1

Next, I ran wgrib2 to confirm that the requested data really exists:

[John.H.Gotway at tfe04 wong_data_20190524]$ wgrib2 WRFTWO.GrbF12 | grep
COLMD
79:29065860:d=2018082000:COLMD:entire atmosphere (considered as a
single layer):12 hour fcst:
[John.H.Gotway at tfe04 wong_data_20190524]$ wgrib2 WRFNAT.GrbF12 | grep
MASSDEN
21:12688735:d=2018081612:MASSDEN:1 hybrid level:12 hour fcst:
...


I'm able to run plot_data_plane on the WRFNAT file fine:

plot_data_plane WRFNAT.GrbF12 MASSDEN.ps 'name="MASSDEN"; level="L1";'

But running the WRFTWO file reveals this error:

plot_data_plane WRFTWO.GrbF12 COLMD.ps 'name="COLMD"; level="L0";'
ERROR  : MetGrib2DataFile::read_grib2_grid() -> found unrecognized
grid definition (32769)

Next, I used wgrib2 to extract that single record #79:

wgrib2 -d 79 -grib_out WRFTWO_COLMD.GrbF12 WRFTWO.GrbF12

And run the result through plot_data_plane but still get the same
error.  Lastly, I ran wgrib2 with the -V verbose option and got:

[John.H.Gotway at tfe04 wong_data_20190524]$ wgrib2 -V
WRFTWO_COLMD.GrbF12
1:0:vt=2018082012:entire atmosphere (considered as a single layer):12
hour fcst:COLMD Column-Integrated Mass Density [kg/m^2]:
    ndata=794802:undef=0:mean=4.73127e-06:min=9.37527e-
22:max=0.00083322
    grid_template=32769:winds(grid):
        I am not an Arakawa E-grid.
        I am rotated but have no rotation angle.
        I am staggered. What am I?
        (953 x 834) units 1e-06 input WE:SN output WE:SN res 56
        lat0 -10.590603 lat-center 54.000000 dlat 121.813000
        lon0 220.914154 lon-center 254.000000 dlon 121.813000
#points=794802

I've never seen this type of output from wgrib2 before!  Suffice it to
say, this isn't a projection type that's supported in MET... or wgrib2
either.  Looks like there's a problem in how the WRFTWO.GrbF12 file is
being created.

Thanks,
John

------------------------------------------------