[Met_help] [rt.rap.ucar.edu #83618] History for Multiple records of AOD in a grib file
Julie Prestopnik via RT
met_help at ucar.edu
Mon Jan 22 10:44:04 MST 2018
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
Hi
I have an input GRIB1 file with multiple records of AOD,before I used
copygb to interpolate that into (-90:90;-180:180} grid to use for
Grid-Stat. I am wondering how to read multiple AOD records in that file, is
there any Grib table specific number or template I can use ?
MET able to read the interpolated file record when I use this command :
plot_data_plane out.grb aod.ps 'name="AOD"; level="A0";' and generated the
attached figure.
Attached are original and interpolated grib1 files.
Thanks.
Partha
----------------------------------------------------------------
Complete Ticket History
----------------------------------------------------------------
Subject: Multiple records of AOD in a grib file
From: John Halley Gotway
Time: Wed Jan 17 11:33:22 2018
Partha,
I'm confused, when I run wgrib on the data you sent, I don't see any
mention of AOD. But I do get a warning message about an undefined
parameter table:
[johnhg at number5]% wgrib out.grb
Undefined parameter table (center 7-0 table 141), using NCEP-opn
1:0:d=15062500:MSLMA:kpds5=129:kpds6=200:kpds7=0:TR=10:P1=0:P2=0:TimeU=1:atmos
col:anl:NAve=0
2:129684:d=15062500:KX:kpds5=133:kpds6=200:kpds7=0:TR=10:P1=0:P2=0:TimeU=1:atmos
col:anl:NAve=0
3:259368:d=15062500:SX:kpds5=134:kpds6=200:kpds7=0:TR=10:P1=0:P2=0:TimeU=1:atmos
col:anl:NAve=0
4:389052:d=15062500:MCONV:kpds5=135:kpds6=200:kpds7=0:TR=10:P1=0:P2=0:TimeU=1:atmos
col:anl:NAve=0
5:518736:d=15062500:CWAT:kpds5=76:kpds6=200:kpds7=0:TR=10:P1=0:P2=0:TimeU=1:atmos
col:anl:NAve=0
6:648420:d=15062500:TSLSA:kpds5=137:kpds6=200:kpds7=0:TR=10:P1=0:P2=0:TimeU=1:atmos
col:anl:NAve=0
All of these 6 records appear to have different GRIB codes. Based on
the
question you asked, I was expecting to find multiple records with the
*same* GRIB code but were different in some other way (i.e. different
level
type or different statistical processing type).
You say these are all AOD, but the plot_data_plane command you sent
does
*not* work for me:
plot_data_plane out.grb aod.ps 'name="AOD"; level="A0";'
ERROR : plot_data_plane -> trouble getting field "name="AOD";
level="A0";"
from file "out.grb"
Please check to see if you have the $MET_GRIB_TABLES environment
variable
set. If so, please send me the table(s) to which it points.
Thanks,
John
On Wed, Jan 17, 2018 at 10:35 AM, Partha Bhattacharjee - NOAA
Affiliate via
RT <met_help at ucar.edu> wrote:
>
> Wed Jan 17 10:35:05 2018: Request 83618 was acted upon.
> Transaction: Ticket created by partha.bhattacharjee at noaa.gov
> Queue: met_help
> Subject: Multiple records of AOD in a grib file
> Owner: Nobody
> Requestors: partha.bhattacharjee at noaa.gov
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=83618 >
>
>
> Hi
>
> I have an input GRIB1 file with multiple records of AOD,before I
used
> copygb to interpolate that into (-90:90;-180:180} grid to use for
> Grid-Stat. I am wondering how to read multiple AOD records in that
file, is
> there any Grib table specific number or template I can use ?
>
> MET able to read the interpolated file record when I use this
command :
> plot_data_plane out.grb aod.ps 'name="AOD"; level="A0";' and
generated
> the
> attached figure.
>
> Attached are original and interpolated grib1 files.
>
> Thanks.
> Partha
>
>
------------------------------------------------
Subject: Multiple records of AOD in a grib file
From: Partha Bhattacharjee - NOAA Affiliate
Time: Wed Jan 17 12:05:49 2018
John,
That is strange. When I did wgrib to out.grib I am getting the
following :
wgrib out.grb
1:0:d=15062500:AOD:kpds5=129:kpds6=200:kpds7=0:TR=10:P1=0:P2=0:TimeU=1:atmos
col:anl:NAve=0
2:129684:d=15062500:var133:kpds5=133:kpds6=200:kpds7=0:TR=10:P1=0:P2=0:TimeU=1:atmos
col:anl:NAve=0
3:259368:d=15062500:var134:kpds5=134:kpds6=200:kpds7=0:TR=10:P1=0:P2=0:TimeU=1:atmos
col:anl:NAve=0
4:389052:d=15062500:var135:kpds5=135:kpds6=200:kpds7=0:TR=10:P1=0:P2=0:TimeU=1:atmos
col:anl:NAve=0
5:518736:d=15062500:CWAT:kpds5=76:kpds6=200:kpds7=0:TR=10:P1=0:P2=0:TimeU=1:atmos
col:anl:NAve=0
6:648420:d=15062500:var137:kpds5=137:kpds6=200:kpds7=0:TR=10:P1=0:P2=0:TimeU=1:atmos
col:anl:NAve=0
and I am running in wcoss and wgrib points to
which wgrib
alias wgrib='/nwprod/util/exec/wgrib'
/nwprod/util/exec/wgrib
So, I used "AOD" in the first record to plot like :
plot_data_plane out.grb aod.ps 'name="AOD"; level="L0";' -v 6
DEBUG 1: Opening data file: out.grb
DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() -> created new
Met2dDataFile object of type "FileType_Gb1".
DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo object
of
type "FileType_Gb1".
DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found range match for
VarInfo "AOD/L0" in GRIB record 1 of GRIB file "out.grb".
DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 1 GRIB records
matching VarInfo "AOD/L0" in GRIB file "out.grb".
DEBUG 1: Creating postscript file: aod.ps
DEBUG 4: draw_map() -> Plotting map data file:
/awc/noscrub/Julie.Prestopnik/met/6.0/share/met/map/country_data
DEBUG 4: draw_map() -> Plotting map data file:
/awc/noscrub/Julie.Prestopnik/met/6.0/share/met/map/usa_state_data
DEBUG 4: draw_map() -> Plotting map data file:
/awc/noscrub/Julie.Prestopnik/met/6.0/share/met/map/major_lakes_data
I do not think I have $MET_GRIB_TABLES environment variable set.
Partha
On Wed, Jan 17, 2018 at 1:33 PM, John Halley Gotway via RT <
met_help at ucar.edu> wrote:
> Partha,
>
> I'm confused, when I run wgrib on the data you sent, I don't see any
> mention of AOD. But I do get a warning message about an undefined
> parameter table:
>
> [johnhg at number5]% wgrib out.grb
>
> Undefined parameter table (center 7-0 table 141), using NCEP-opn
> 1:0:d=15062500:MSLMA:kpds5=129:kpds6=200:kpds7=0:TR=10:
> P1=0:P2=0:TimeU=1:atmos
> col:anl:NAve=0
> 2:129684:d=15062500:KX:kpds5=133:kpds6=200:kpds7=0:TR=10:
> P1=0:P2=0:TimeU=1:atmos
> col:anl:NAve=0
> 3:259368:d=15062500:SX:kpds5=134:kpds6=200:kpds7=0:TR=10:
> P1=0:P2=0:TimeU=1:atmos
> col:anl:NAve=0
> 4:389052:d=15062500:MCONV:kpds5=135:kpds6=200:kpds7=0:
> TR=10:P1=0:P2=0:TimeU=1:atmos
> col:anl:NAve=0
> 5:518736:d=15062500:CWAT:kpds5=76:kpds6=200:kpds7=0:TR=
> 10:P1=0:P2=0:TimeU=1:atmos
> col:anl:NAve=0
> 6:648420:d=15062500:TSLSA:kpds5=137:kpds6=200:kpds7=0:
> TR=10:P1=0:P2=0:TimeU=1:atmos
> col:anl:NAve=0
>
> All of these 6 records appear to have different GRIB codes. Based
on the
> question you asked, I was expecting to find multiple records with
the
> *same* GRIB code but were different in some other way (i.e.
different level
> type or different statistical processing type).
>
> You say these are all AOD, but the plot_data_plane command you sent
does
> *not* work for me:
>
> plot_data_plane out.grb aod.ps 'name="AOD"; level="A0";'
> ERROR : plot_data_plane -> trouble getting field "name="AOD";
level="A0";"
> from file "out.grb"
>
> Please check to see if you have the $MET_GRIB_TABLES environment
variable
> set. If so, please send me the table(s) to which it points.
>
> Thanks,
> John
>
>
> On Wed, Jan 17, 2018 at 10:35 AM, Partha Bhattacharjee - NOAA
Affiliate via
> RT <met_help at ucar.edu> wrote:
>
> >
> > Wed Jan 17 10:35:05 2018: Request 83618 was acted upon.
> > Transaction: Ticket created by partha.bhattacharjee at noaa.gov
> > Queue: met_help
> > Subject: Multiple records of AOD in a grib file
> > Owner: Nobody
> > Requestors: partha.bhattacharjee at noaa.gov
> > Status: new
> > Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=83618 >
> >
> >
> > Hi
> >
> > I have an input GRIB1 file with multiple records of AOD,before I
used
> > copygb to interpolate that into (-90:90;-180:180} grid to use for
> > Grid-Stat. I am wondering how to read multiple AOD records in that
file,
> is
> > there any Grib table specific number or template I can use ?
> >
> > MET able to read the interpolated file record when I use this
command :
> > plot_data_plane out.grb aod.ps 'name="AOD"; level="A0";' and
generated
> > the
> > attached figure.
> >
> > Attached are original and interpolated grib1 files.
> >
> > Thanks.
> > Partha
> >
> >
>
>
--
Partha S. Bhattacharjee, Ph.D.
I.M. Systems Group
Environmental Modeling Centre
NCWCP
5830 University Research Court.
College Park, MD 20740
------------------------------------------------
Subject: Multiple records of AOD in a grib file
From: John Halley Gotway
Time: Wed Jan 17 14:26:40 2018
Partha,
Ah, OK. So your file contains 6 GRIB records, each of which has a
different GRIB code (129, 133, 134, 135, 76, and 137). Those are the
kpds5
values listed in the output of wgrib. Running wgrib -V reveals that 5
of
the 6 records are encoded as:
center 7 subcenter 0 process 82 Table 141 scan: WE:SN winds(N/S)
So this is NCEP GRIB1 table number 141, and that can be found here:
met-6.1/share/met/table_files/grib1_ncep_141_7.txt
Here's a selection of lines from that file:
129 141 7 -1 "AOD" "Aerosol Optical Depth" ""
133 141 7 -1 "var133" "undefined" ""
134 141 7 -1 "var134" "undefined" ""
135 141 7 -1 "var135" "undefined" ""
137 141 7 -1 "var137" "undefined" ""
And this info is taken from this GRIB table:
http://www.nco.ncep.noaa.gov/pmb/docs/on388/table2.html#TABLE141
So you could just refer to this data using those "var133" (and so on)
string:
plot_data_plane out.grb plot_rec1.ps 'name="var133"; level="L0";'
Alternatively, you could specifying the "GRIB1_code" value instead of
the
name:
plot_data_plane out.grb plot_rec1.ps 'GRIB1_code=133; level="L0";'
Both approaches would result in the same plot. Presumably, since you
have
data for codes 133, 134, 135, and 137, perhaps these codes are no
longer
"reserved for future use", as indicated on the NCO website.
Hope that helps clarify.
Thanks,
John
On Wed, Jan 17, 2018 at 12:05 PM, Partha Bhattacharjee - NOAA
Affiliate via
RT <met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=83618 >
>
> John,
>
> That is strange. When I did wgrib to out.grib I am getting the
following :
>
> wgrib out.grb
> 1:0:d=15062500:AOD:kpds5=129:kpds6=200:kpds7=0:TR=10:P1=0:
> P2=0:TimeU=1:atmos
> col:anl:NAve=0
> 2:129684:d=15062500:var133:kpds5=133:kpds6=200:kpds7=0:
> TR=10:P1=0:P2=0:TimeU=1:atmos
> col:anl:NAve=0
> 3:259368:d=15062500:var134:kpds5=134:kpds6=200:kpds7=0:
> TR=10:P1=0:P2=0:TimeU=1:atmos
> col:anl:NAve=0
> 4:389052:d=15062500:var135:kpds5=135:kpds6=200:kpds7=0:
> TR=10:P1=0:P2=0:TimeU=1:atmos
> col:anl:NAve=0
> 5:518736:d=15062500:CWAT:kpds5=76:kpds6=200:kpds7=0:TR=
> 10:P1=0:P2=0:TimeU=1:atmos
> col:anl:NAve=0
> 6:648420:d=15062500:var137:kpds5=137:kpds6=200:kpds7=0:
> TR=10:P1=0:P2=0:TimeU=1:atmos
> col:anl:NAve=0
>
> and I am running in wcoss and wgrib points to
> which wgrib
> alias wgrib='/nwprod/util/exec/wgrib'
> /nwprod/util/exec/wgrib
>
> So, I used "AOD" in the first record to plot like :
>
> plot_data_plane out.grb aod.ps 'name="AOD"; level="L0";' -v 6
> DEBUG 1: Opening data file: out.grb
> DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() -> created new
> Met2dDataFile object of type "FileType_Gb1".
> DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo
object of
> type "FileType_Gb1".
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found range match
for
> VarInfo "AOD/L0" in GRIB record 1 of GRIB file "out.grb".
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 1 GRIB
records
> matching VarInfo "AOD/L0" in GRIB file "out.grb".
> DEBUG 1: Creating postscript file: aod.ps
> DEBUG 4: draw_map() -> Plotting map data file:
> /awc/noscrub/Julie.Prestopnik/met/6.0/share/met/map/country_data
> DEBUG 4: draw_map() -> Plotting map data file:
> /awc/noscrub/Julie.Prestopnik/met/6.0/share/met/map/usa_state_data
> DEBUG 4: draw_map() -> Plotting map data file:
> /awc/noscrub/Julie.Prestopnik/met/6.0/share/met/map/major_lakes_data
>
> I do not think I have $MET_GRIB_TABLES environment variable set.
>
> Partha
>
> On Wed, Jan 17, 2018 at 1:33 PM, John Halley Gotway via RT <
> met_help at ucar.edu> wrote:
>
> > Partha,
> >
> > I'm confused, when I run wgrib on the data you sent, I don't see
any
> > mention of AOD. But I do get a warning message about an undefined
> > parameter table:
> >
> > [johnhg at number5]% wgrib out.grb
> >
> > Undefined parameter table (center 7-0 table 141), using NCEP-opn
> > 1:0:d=15062500:MSLMA:kpds5=129:kpds6=200:kpds7=0:TR=10:
> > P1=0:P2=0:TimeU=1:atmos
> > col:anl:NAve=0
> > 2:129684:d=15062500:KX:kpds5=133:kpds6=200:kpds7=0:TR=10:
> > P1=0:P2=0:TimeU=1:atmos
> > col:anl:NAve=0
> > 3:259368:d=15062500:SX:kpds5=134:kpds6=200:kpds7=0:TR=10:
> > P1=0:P2=0:TimeU=1:atmos
> > col:anl:NAve=0
> > 4:389052:d=15062500:MCONV:kpds5=135:kpds6=200:kpds7=0:
> > TR=10:P1=0:P2=0:TimeU=1:atmos
> > col:anl:NAve=0
> > 5:518736:d=15062500:CWAT:kpds5=76:kpds6=200:kpds7=0:TR=
> > 10:P1=0:P2=0:TimeU=1:atmos
> > col:anl:NAve=0
> > 6:648420:d=15062500:TSLSA:kpds5=137:kpds6=200:kpds7=0:
> > TR=10:P1=0:P2=0:TimeU=1:atmos
> > col:anl:NAve=0
> >
> > All of these 6 records appear to have different GRIB codes. Based
on the
> > question you asked, I was expecting to find multiple records with
the
> > *same* GRIB code but were different in some other way (i.e.
different
> level
> > type or different statistical processing type).
> >
> > You say these are all AOD, but the plot_data_plane command you
sent does
> > *not* work for me:
> >
> > plot_data_plane out.grb aod.ps 'name="AOD"; level="A0";'
> > ERROR : plot_data_plane -> trouble getting field "name="AOD";
> level="A0";"
> > from file "out.grb"
> >
> > Please check to see if you have the $MET_GRIB_TABLES environment
variable
> > set. If so, please send me the table(s) to which it points.
> >
> > Thanks,
> > John
> >
> >
> > On Wed, Jan 17, 2018 at 10:35 AM, Partha Bhattacharjee - NOAA
Affiliate
> via
> > RT <met_help at ucar.edu> wrote:
> >
> > >
> > > Wed Jan 17 10:35:05 2018: Request 83618 was acted upon.
> > > Transaction: Ticket created by partha.bhattacharjee at noaa.gov
> > > Queue: met_help
> > > Subject: Multiple records of AOD in a grib file
> > > Owner: Nobody
> > > Requestors: partha.bhattacharjee at noaa.gov
> > > Status: new
> > > Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=83618
> >
> > >
> > >
> > > Hi
> > >
> > > I have an input GRIB1 file with multiple records of AOD,before I
used
> > > copygb to interpolate that into (-90:90;-180:180} grid to use
for
> > > Grid-Stat. I am wondering how to read multiple AOD records in
that
> file,
> > is
> > > there any Grib table specific number or template I can use ?
> > >
> > > MET able to read the interpolated file record when I use this
command :
> > > plot_data_plane out.grb aod.ps 'name="AOD"; level="A0";' and
> generated
> > > the
> > > attached figure.
> > >
> > > Attached are original and interpolated grib1 files.
> > >
> > > Thanks.
> > > Partha
> > >
> > >
> >
> >
>
>
> --
> Partha S. Bhattacharjee, Ph.D.
> I.M. Systems Group
> Environmental Modeling Centre
> NCWCP
> 5830 University Research Court.
> College Park, MD 20740
>
>
------------------------------------------------
Subject: Multiple records of AOD in a grib file
From: Partha Bhattacharjee - NOAA Affiliate
Time: Wed Jan 17 14:34:40 2018
Thanks John. Related to my earlier email, I am getting error message
when
trying to use MET Series Stat. I am trying to compare regridded
forecast
file and observational file but getting the following error :
GSL_RNG_TYPE=mt19937
GSL_RNG_SEED=18446744073587810138
terminate called after throwing an instance of 'std::logic_error'
what(): basic_string::_S_construct NULL not valid
I have attached the output.txt along with this email. My files are in
WCOSS
at : /global/noscrub/Partha.Bhattacharjee/MET_Test. Original Grib
files are
in the Fcst directory : "AOD550.Grb.t00z.f00" and interpolated ones
are
ones with "grb" at the end. I have used copygb to interpolate these
input
files which are at 384 X 190 to 360 X 180 to make them same as Obs
files
are. But not sure if the error is due to not matching grids or from
something else.
Thanks.
Partha
On Wed, Jan 17, 2018 at 4:26 PM, John Halley Gotway via RT <
met_help at ucar.edu> wrote:
> Partha,
>
> Ah, OK. So your file contains 6 GRIB records, each of which has a
> different GRIB code (129, 133, 134, 135, 76, and 137). Those are
the kpds5
> values listed in the output of wgrib. Running wgrib -V reveals that
5 of
> the 6 records are encoded as:
> center 7 subcenter 0 process 82 Table 141 scan: WE:SN winds(N/S)
>
> So this is NCEP GRIB1 table number 141, and that can be found here:
> met-6.1/share/met/table_files/grib1_ncep_141_7.txt
>
> Here's a selection of lines from that file:
>
> 129 141 7 -1 "AOD" "Aerosol Optical Depth" ""
> 133 141 7 -1 "var133" "undefined" ""
> 134 141 7 -1 "var134" "undefined" ""
> 135 141 7 -1 "var135" "undefined" ""
> 137 141 7 -1 "var137" "undefined" ""
>
> And this info is taken from this GRIB table:
> http://www.nco.ncep.noaa.gov/pmb/docs/on388/table2.html#TABLE141
>
> So you could just refer to this data using those "var133" (and so
on)
> string:
>
> plot_data_plane out.grb plot_rec1.ps 'name="var133";
level="L0";'
>
> Alternatively, you could specifying the "GRIB1_code" value instead
of the
> name:
>
> plot_data_plane out.grb plot_rec1.ps 'GRIB1_code=133;
level="L0";'
>
> Both approaches would result in the same plot. Presumably, since
you have
> data for codes 133, 134, 135, and 137, perhaps these codes are no
longer
> "reserved for future use", as indicated on the NCO website.
>
> Hope that helps clarify.
>
> Thanks,
> John
>
> On Wed, Jan 17, 2018 at 12:05 PM, Partha Bhattacharjee - NOAA
Affiliate via
> RT <met_help at ucar.edu> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=83618 >
> >
> > John,
> >
> > That is strange. When I did wgrib to out.grib I am getting the
following
> :
> >
> > wgrib out.grb
> > 1:0:d=15062500:AOD:kpds5=129:kpds6=200:kpds7=0:TR=10:P1=0:
> > P2=0:TimeU=1:atmos
> > col:anl:NAve=0
> > 2:129684:d=15062500:var133:kpds5=133:kpds6=200:kpds7=0:
> > TR=10:P1=0:P2=0:TimeU=1:atmos
> > col:anl:NAve=0
> > 3:259368:d=15062500:var134:kpds5=134:kpds6=200:kpds7=0:
> > TR=10:P1=0:P2=0:TimeU=1:atmos
> > col:anl:NAve=0
> > 4:389052:d=15062500:var135:kpds5=135:kpds6=200:kpds7=0:
> > TR=10:P1=0:P2=0:TimeU=1:atmos
> > col:anl:NAve=0
> > 5:518736:d=15062500:CWAT:kpds5=76:kpds6=200:kpds7=0:TR=
> > 10:P1=0:P2=0:TimeU=1:atmos
> > col:anl:NAve=0
> > 6:648420:d=15062500:var137:kpds5=137:kpds6=200:kpds7=0:
> > TR=10:P1=0:P2=0:TimeU=1:atmos
> > col:anl:NAve=0
> >
> > and I am running in wcoss and wgrib points to
> > which wgrib
> > alias wgrib='/nwprod/util/exec/wgrib'
> > /nwprod/util/exec/wgrib
> >
> > So, I used "AOD" in the first record to plot like :
> >
> > plot_data_plane out.grb aod.ps 'name="AOD"; level="L0";' -v 6
> > DEBUG 1: Opening data file: out.grb
> > DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() -> created
new
> > Met2dDataFile object of type "FileType_Gb1".
> > DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo
object of
> > type "FileType_Gb1".
> > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found range match
for
> > VarInfo "AOD/L0" in GRIB record 1 of GRIB file "out.grb".
> > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 1 GRIB
records
> > matching VarInfo "AOD/L0" in GRIB file "out.grb".
> > DEBUG 1: Creating postscript file: aod.ps
> > DEBUG 4: draw_map() -> Plotting map data file:
> > /awc/noscrub/Julie.Prestopnik/met/6.0/share/met/map/country_data
> > DEBUG 4: draw_map() -> Plotting map data file:
> > /awc/noscrub/Julie.Prestopnik/met/6.0/share/met/map/usa_state_data
> > DEBUG 4: draw_map() -> Plotting map data file:
> >
/awc/noscrub/Julie.Prestopnik/met/6.0/share/met/map/major_lakes_data
> >
> > I do not think I have $MET_GRIB_TABLES environment variable set.
> >
> > Partha
> >
> > On Wed, Jan 17, 2018 at 1:33 PM, John Halley Gotway via RT <
> > met_help at ucar.edu> wrote:
> >
> > > Partha,
> > >
> > > I'm confused, when I run wgrib on the data you sent, I don't see
any
> > > mention of AOD. But I do get a warning message about an
undefined
> > > parameter table:
> > >
> > > [johnhg at number5]% wgrib out.grb
> > >
> > > Undefined parameter table (center 7-0 table 141), using NCEP-opn
> > > 1:0:d=15062500:MSLMA:kpds5=129:kpds6=200:kpds7=0:TR=10:
> > > P1=0:P2=0:TimeU=1:atmos
> > > col:anl:NAve=0
> > > 2:129684:d=15062500:KX:kpds5=133:kpds6=200:kpds7=0:TR=10:
> > > P1=0:P2=0:TimeU=1:atmos
> > > col:anl:NAve=0
> > > 3:259368:d=15062500:SX:kpds5=134:kpds6=200:kpds7=0:TR=10:
> > > P1=0:P2=0:TimeU=1:atmos
> > > col:anl:NAve=0
> > > 4:389052:d=15062500:MCONV:kpds5=135:kpds6=200:kpds7=0:
> > > TR=10:P1=0:P2=0:TimeU=1:atmos
> > > col:anl:NAve=0
> > > 5:518736:d=15062500:CWAT:kpds5=76:kpds6=200:kpds7=0:TR=
> > > 10:P1=0:P2=0:TimeU=1:atmos
> > > col:anl:NAve=0
> > > 6:648420:d=15062500:TSLSA:kpds5=137:kpds6=200:kpds7=0:
> > > TR=10:P1=0:P2=0:TimeU=1:atmos
> > > col:anl:NAve=0
> > >
> > > All of these 6 records appear to have different GRIB codes.
Based on
> the
> > > question you asked, I was expecting to find multiple records
with the
> > > *same* GRIB code but were different in some other way (i.e.
different
> > level
> > > type or different statistical processing type).
> > >
> > > You say these are all AOD, but the plot_data_plane command you
sent
> does
> > > *not* work for me:
> > >
> > > plot_data_plane out.grb aod.ps 'name="AOD"; level="A0";'
> > > ERROR : plot_data_plane -> trouble getting field "name="AOD";
> > level="A0";"
> > > from file "out.grb"
> > >
> > > Please check to see if you have the $MET_GRIB_TABLES environment
> variable
> > > set. If so, please send me the table(s) to which it points.
> > >
> > > Thanks,
> > > John
> > >
> > >
> > > On Wed, Jan 17, 2018 at 10:35 AM, Partha Bhattacharjee - NOAA
Affiliate
> > via
> > > RT <met_help at ucar.edu> wrote:
> > >
> > > >
> > > > Wed Jan 17 10:35:05 2018: Request 83618 was acted upon.
> > > > Transaction: Ticket created by partha.bhattacharjee at noaa.gov
> > > > Queue: met_help
> > > > Subject: Multiple records of AOD in a grib file
> > > > Owner: Nobody
> > > > Requestors: partha.bhattacharjee at noaa.gov
> > > > Status: new
> > > > Ticket <URL: https://rt.rap.ucar.edu/rt/
> Ticket/Display.html?id=83618
> > >
> > > >
> > > >
> > > > Hi
> > > >
> > > > I have an input GRIB1 file with multiple records of AOD,before
I used
> > > > copygb to interpolate that into (-90:90;-180:180} grid to use
for
> > > > Grid-Stat. I am wondering how to read multiple AOD records in
that
> > file,
> > > is
> > > > there any Grib table specific number or template I can use ?
> > > >
> > > > MET able to read the interpolated file record when I use this
> command :
> > > > plot_data_plane out.grb aod.ps 'name="AOD"; level="A0";' and
> > generated
> > > > the
> > > > attached figure.
> > > >
> > > > Attached are original and interpolated grib1 files.
> > > >
> > > > Thanks.
> > > > Partha
> > > >
> > > >
> > >
> > >
> >
> >
> > --
> > Partha S. Bhattacharjee, Ph.D.
> > I.M. Systems Group
> > Environmental Modeling Centre
> > NCWCP
> > 5830 University Research Court.
> > College Park, MD 20740
> >
> >
>
>
--
Partha S. Bhattacharjee, Ph.D.
I.M. Systems Group
Environmental Modeling Centre
NCWCP
5830 University Research Court.
College Park, MD 20740
------------------------------------------------
Subject: Multiple records of AOD in a grib file
From: Partha Bhattacharjee - NOAA Affiliate
Time: Wed Jan 17 14:34:40 2018
DEBUG 1: Default Config File:
/usrx/local/dev/MET/6.0/share/met/config/SeriesAnalysisConfig_default
DEBUG 1: User Config File: SeriesAnlConfig_fh018Test
DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() -> created new
Met2dDataFile object of type "FileType_Gb1".
DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() -> created new
Met2dDataFile object of type "FileType_NcCF".
DEBUG 4: NcCfFile::open() -> parsing units for the time variable
"seconds since 2015-06-28 00:00:00 UTC"
DEBUG 4: parse_cf_time_string() -> parsed NetCDF CF convention time
unit string "seconds since 2015-06-28 00:00:00 UTC" as a reference
time of 20150628_000000 and 1 second(s) per time step.
DEBUG 4: NcCfFile::open() -> could not extract init time from the
"forecast_reference_time" variable.
DEBUG 4: NcCfFile::open() -> could not extract init time from file
name.
DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() -> created new
Met2dDataFile object of type "FileType_Gb1".
DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() -> created new
Met2dDataFile object of type "FileType_NcCF".
DEBUG 4: NcCfFile::open() -> parsing units for the time variable
"seconds since 2015-06-28 00:00:00 UTC"
DEBUG 4: parse_cf_time_string() -> parsed NetCDF CF convention time
unit string "seconds since 2015-06-28 00:00:00 UTC" as a reference
time of 20150628_000000 and 1 second(s) per time step.
DEBUG 4: NcCfFile::open() -> could not extract init time from the
"forecast_reference_time" variable.
DEBUG 4: NcCfFile::open() -> could not extract init time from file
name.
Filename ...
| Nelements = 3
| Nalloc = 50
| MaxLength = 69
| IgnoreCase = 0
| | Element # 0 =
"/usrx/local/dev/MET/6.0/share/met/config/ConfigConstants"
| | Element # 1 =
"/usrx/local/dev/MET/6.0/share/met/config/SeriesAnalysisConfig_default"
| | Element # 2 = "SeriesAnlConfig_fh018Test"
Debug = false
Config File Entries ...
| Nentries = 79
| Address = 0x8a9f88
| Parent = 0
| IsArray = false
| Dictionary Entry[0] ...
| | Name = exit_on_warning
| | Type = BooleanType
| | Boolean Value = false
| Dictionary Entry[1] ...
| | Name = nc_compression
| | Type = IntegerType
| | Integer Value = 0
| Dictionary Entry[2] ...
| | Name = output_precision
| | Type = IntegerType
| | Integer Value = 5
| Dictionary Entry[3] ...
| | Name = tmp_dir
| | Type = StringType
| | String Value = "/tmp"
| Dictionary Entry[4] ...
| | Name = GRIB1
| | Type = IntegerType
| | Integer Value = 1
| Dictionary Entry[5] ...
| | Name = GRIB2
| | Type = IntegerType
| | Integer Value = 2
| Dictionary Entry[6] ...
| | Name = NETCDF_MET
| | Type = IntegerType
| | Integer Value = 3
| Dictionary Entry[7] ...
| | Name = NETCDF_PINT
| | Type = IntegerType
| | Integer Value = 4
| Dictionary Entry[8] ...
| | Name = NETCDF_NCCF
| | Type = IntegerType
| | Integer Value = 5
| Dictionary Entry[9] ...
| | Name = NONE
| | Type = IntegerType
| | Integer Value = 1
| Dictionary Entry[10] ...
| | Name = UNION
| | Type = IntegerType
| | Integer Value = 2
| Dictionary Entry[11] ...
| | Name = INTERSECTION
| | Type = IntegerType
| | Integer Value = 3
| Dictionary Entry[12] ...
| | Name = SYMDIFF
| | Type = IntegerType
| | Integer Value = 4
| Dictionary Entry[13] ...
| | Name = NEAREST
| | Type = IntegerType
| | Integer Value = 9
| Dictionary Entry[14] ...
| | Name = MEDIAN
| | Type = IntegerType
| | Integer Value = 3
| Dictionary Entry[15] ...
| | Name = UW_MEAN
| | Type = IntegerType
| | Integer Value = 4
| Dictionary Entry[16] ...
| | Name = DW_MEAN
| | Type = IntegerType
| | Integer Value = 5
| Dictionary Entry[17] ...
| | Name = LS_FIT
| | Type = IntegerType
| | Integer Value = 6
| Dictionary Entry[18] ...
| | Name = BILIN
| | Type = IntegerType
| | Integer Value = 7
| Dictionary Entry[19] ...
| | Name = NBRHD
| | Type = IntegerType
| | Integer Value = 8
| Dictionary Entry[20] ...
| | Name = BUDGET
| | Type = IntegerType
| | Integer Value = 10
| Dictionary Entry[21] ...
| | Name = FORCE
| | Type = IntegerType
| | Integer Value = 11
| Dictionary Entry[22] ...
| | Name = BEST
| | Type = IntegerType
| | Integer Value = 12
| Dictionary Entry[23] ...
| | Name = UPPER_LEFT
| | Type = IntegerType
| | Integer Value = 13
| Dictionary Entry[24] ...
| | Name = UPPER_RIGHT
| | Type = IntegerType
| | Integer Value = 14
| Dictionary Entry[25] ...
| | Name = LOWER_RIGHT
| | Type = IntegerType
| | Integer Value = 15
| Dictionary Entry[26] ...
| | Name = LOWER_LEFT
| | Type = IntegerType
| | Integer Value = 16
| Dictionary Entry[27] ...
| | Name = PERC
| | Type = IntegerType
| | Integer Value = 17
| Dictionary Entry[28] ...
| | Name = MIN
| | Type = IntegerType
| | Integer Value = 18
| Dictionary Entry[29] ...
| | Name = MAX
| | Type = IntegerType
| | Integer Value = 19
| Dictionary Entry[30] ...
| | Name = BOTH
| | Type = IntegerType
| | Integer Value = 2
| Dictionary Entry[31] ...
| | Name = FCST
| | Type = IntegerType
| | Integer Value = 3
| Dictionary Entry[32] ...
| | Name = OBS
| | Type = IntegerType
| | Integer Value = 4
| Dictionary Entry[33] ...
| | Name = SQUARE
| | Type = IntegerType
| | Integer Value = 1
| Dictionary Entry[34] ...
| | Name = CIRCLE
| | Type = IntegerType
| | Integer Value = 2
| Dictionary Entry[35] ...
| | Name = ADECK
| | Type = IntegerType
| | Integer Value = 3
| Dictionary Entry[36] ...
| | Name = BDECK
| | Type = IntegerType
| | Integer Value = 4
| Dictionary Entry[37] ...
| | Name = FILL
| | Type = IntegerType
| | Integer Value = 2
| Dictionary Entry[38] ...
| | Name = REPLACE
| | Type = IntegerType
| | Integer Value = 3
| Dictionary Entry[39] ...
| | Name = STAT
| | Type = IntegerType
| | Integer Value = 3
| Dictionary Entry[40] ...
| | Name = PCTILE
| | Type = IntegerType
| | Integer Value = 1
| Dictionary Entry[41] ...
| | Name = BCA
| | Type = IntegerType
| | Integer Value = 2
| Dictionary Entry[42] ...
| | Name = COS_LAT
| | Type = IntegerType
| | Integer Value = 2
| Dictionary Entry[43] ...
| | Name = AREA
| | Type = IntegerType
| | Integer Value = 3
| Dictionary Entry[44] ...
| | Name = UNIQUE
| | Type = IntegerType
| | Integer Value = 2
| Dictionary Entry[45] ...
| | Name = AUTO
| | Type = IntegerType
| | Integer Value = 2
| Dictionary Entry[46] ...
| | Name = TILE
| | Type = IntegerType
| | Integer Value = 3
| Dictionary Entry[47] ...
| | Name = PAD
| | Type = IntegerType
| | Integer Value = 4
| Dictionary Entry[48] ...
| | Name = HAAR
| | Type = IntegerType
| | Integer Value = 2
| Dictionary Entry[49] ...
| | Name = HAAR_CNTR
| | Type = IntegerType
| | Integer Value = 3
| Dictionary Entry[50] ...
| | Name = DAUB
| | Type = IntegerType
| | Integer Value = 4
| Dictionary Entry[51] ...
| | Name = DAUB_CNTR
| | Type = IntegerType
| | Integer Value = 5
| Dictionary Entry[52] ...
| | Name = BSPLINE
| | Type = IntegerType
| | Integer Value = 6
| Dictionary Entry[53] ...
| | Name = BSPLINE_CNTR
| | Type = IntegerType
| | Integer Value = 7
| Dictionary Entry[54] ...
| | Name = THRESH
| | Type = IntegerType
| | Integer Value = 3
| Dictionary Entry[55] ...
| | Name = ENGINE
| | Type = IntegerType
| | Integer Value = 4
| Dictionary Entry[56] ...
| | Name = MERGE_BOTH
| | Type = IntegerType
| | Integer Value = 2
| Dictionary Entry[57] ...
| | Name = MERGE_FCST
| | Type = IntegerType
| | Integer Value = 3
| Dictionary Entry[58] ...
| | Name = NO_MERGE
| | Type = IntegerType
| | Integer Value = 4
| Dictionary Entry[59] ...
| | Name = TOP
| | Type = BooleanType
| | Boolean Value = true
| Dictionary Entry[60] ...
| | Name = BOTTOM
| | Type = BooleanType
| | Boolean Value = false
| Dictionary Entry[61] ...
| | Name = model
| | Type = StringType
| | String Value = "WRF"
| Dictionary Entry[62] ...
| | Name = desc
| | Type = StringType
| | String Value = "NA"
| Dictionary Entry[63] ...
| | Name = obtype
| | Type = StringType
| | String Value = "ANALYS"
| Dictionary Entry[64] ...
| | Name = regrid
| | Type = DictionaryType
| | Dict/Array Value =
| | | Nentries = 4
| | | Address = 0x2a34bc0
| | | Parent = 0x8a9f88
| | | IsArray = false
| | | Dictionary Entry[0] ...
| | | | Name = to_grid
| | | | Type = IntegerType
| | | | Integer Value = 1
| | | Dictionary Entry[1] ...
| | | | Name = method
| | | | Type = IntegerType
| | | | Integer Value = 9
| | | Dictionary Entry[2] ...
| | | | Name = width
| | | | Type = IntegerType
| | | | Integer Value = 1
| | | Dictionary Entry[3] ...
| | | | Name = vld_thresh
| | | | Type = FloatType
| | | | Float Value = 0.5
| Dictionary Entry[65] ...
| | Name = cat_thresh
| | Type = ArrayType
| | Dict/Array Value =
| | | Nentries = 1
| | | Address = 0x2d6d130
| | | Parent = 0x8a9f88
| | | IsArray = true
| | | Dictionary Entry[0] ...
| | | | Name = (nul)
| | | | Type = ThresholdType
| | | | Thresh Value =
| | | | | type = 0 (NA)
| | | | | value = -9999
| Dictionary Entry[66] ...
| | Name = cnt_thresh
| | Type = ArrayType
| | Dict/Array Value =
| | | Nentries = 1
| | | Address = 0x2771e30
| | | Parent = 0x8a9f88
| | | IsArray = true
| | | Dictionary Entry[0] ...
| | | | Name = (nul)
| | | | Type = ThresholdType
| | | | Thresh Value =
| | | | | type = 0 (NA)
| | | | | value = -9999
| Dictionary Entry[67] ...
| | Name = cnt_logic
| | Type = IntegerType
| | Integer Value = 3
| Dictionary Entry[68] ...
| | Name = fcst
| | Type = DictionaryType
| | Dict/Array Value =
| | | Nentries = 1
| | | Address = 0x2dfde70
| | | Parent = 0x8a9f88
| | | IsArray = false
| | | Dictionary Entry[0] ...
| | | | Name = field
| | | | Type = ArrayType
| | | | Dict/Array Value =
| | | | | Nentries = 1
| | | | | Address = 0x27726d0
| | | | | Parent = 0x2dfde70
| | | | | IsArray = true
| | | | | Dictionary Entry[0] ...
| | | | | | Name = (nul)
| | | | | | Type = DictionaryType
| | | | | | Dict/Array Value =
| | | | | | | Nentries = 2
| | | | | | | Address = 0x2947d30
| | | | | | | Parent = 0x27726d0
| | | | | | | IsArray = false
| | | | | | | Dictionary Entry[0] ...
| | | | | | | | Name = name
| | | | | | | | Type = StringType
| | | | | | | | String Value = "AOD"
| | | | | | | Dictionary Entry[1] ...
| | | | | | | | Name = level
| | | | | | | | Type = ArrayType
| | | | | | | | Dict/Array Value =
| | | | | | | | | Nentries = 12
| | | | | | | | | Address = 0x28fbb70
| | | | | | | | | Parent = 0x2947d30
| | | | | | | | | IsArray = true
| | | | | | | | | Dictionary Entry[0] ...
| | | | | | | | | | Name = (nul)
| | | | | | | | | | Type = StringType
| | | | | | | | | | String Value = "L0"
| | | | | | | | | Dictionary Entry[1] ...
| | | | | | | | | | Name = (nul)
| | | | | | | | | | Type = StringType
| | | | | | | | | | String Value = "L0"
| | | | | | | | | Dictionary Entry[2] ...
| | | | | | | | | | Name = (nul)
| | | | | | | | | | Type = StringType
| | | | | | | | | | String Value = "L0"
| | | | | | | | | Dictionary Entry[3] ...
| | | | | | | | | | Name = (nul)
| | | | | | | | | | Type = StringType
| | | | | | | | | | String Value = "L0"
| | | | | | | | | Dictionary Entry[4] ...
| | | | | | | | | | Name = (nul)
| | | | | | | | | | Type = StringType
| | | | | | | | | | String Value = "L0"
| | | | | | | | | Dictionary Entry[5] ...
| | | | | | | | | | Name = (nul)
| | | | | | | | | | Type = StringType
| | | | | | | | | | String Value = "L0"
| | | | | | | | | Dictionary Entry[6] ...
| | | | | | | | | | Name = (nul)
| | | | | | | | | | Type = StringType
| | | | | | | | | | String Value = "L0"
| | | | | | | | | Dictionary Entry[7] ...
| | | | | | | | | | Name = (nul)
| | | | | | | | | | Type = StringType
| | | | | | | | | | String Value = "L0"
| | | | | | | | | Dictionary Entry[8] ...
| | | | | | | | | | Name = (nul)
| | | | | | | | | | Type = StringType
| | | | | | | | | | String Value = "L0"
| | | | | | | | | Dictionary Entry[9] ...
| | | | | | | | | | Name = (nul)
| | | | | | | | | | Type = StringType
| | | | | | | | | | String Value = "L0"
| | | | | | | | | Dictionary Entry[10] ...
| | | | | | | | | | Name = (nul)
| | | | | | | | | | Type = StringType
| | | | | | | | | | String Value = "L0"
| | | | | | | | | Dictionary Entry[11] ...
| | | | | | | | | | Name = (nul)
| | | | | | | | | | Type = StringType
| | | | | | | | | | String Value = "L0"
| Dictionary Entry[69] ...
| | Name = obs
| | Type = DictionaryType
| | Dict/Array Value =
| | | Nentries = 1
| | | Address = 0x2790ff0
| | | Parent = 0x8a9f88
| | | IsArray = false
| | | Dictionary Entry[0] ...
| | | | Name = field
| | | | Type = ArrayType
| | | | Dict/Array Value =
| | | | | Nentries = 1
| | | | | Address = 0x2b50f30
| | | | | Parent = 0x2790ff0
| | | | | IsArray = true
| | | | | Dictionary Entry[0] ...
| | | | | | Name = (nul)
| | | | | | Type = DictionaryType
| | | | | | Dict/Array Value =
| | | | | | | Nentries = 2
| | | | | | | Address = 0x2dfd1c0
| | | | | | | Parent = 0x2b50f30
| | | | | | | IsArray = false
| | | | | | | Dictionary Entry[0] ...
| | | | | | | | Name = name
| | | | | | | | Type = StringType
| | | | | | | | String Value =
"total_aod_mean"
| | | | | | | Dictionary Entry[1] ...
| | | | | | | | Name = level
| | | | | | | | Type = ArrayType
| | | | | | | | Dict/Array Value =
| | | | | | | | | Nentries = 12
| | | | | | | | | Address = 0x2795780
| | | | | | | | | Parent = 0x2dfd1c0
| | | | | | | | | IsArray = true
| | | | | | | | | Dictionary Entry[0] ...
| | | | | | | | | | Name = (nul)
| | | | | | | | | | Type = StringType
| | | | | | | | | | String Value =
"(0,*,*)"
| | | | | | | | | Dictionary Entry[1] ...
| | | | | | | | | | Name = (nul)
| | | | | | | | | | Type = StringType
| | | | | | | | | | String Value =
"(1,*,*)"
| | | | | | | | | Dictionary Entry[2] ...
| | | | | | | | | | Name = (nul)
| | | | | | | | | | Type = StringType
| | | | | | | | | | String Value =
"(2,*,*)"
| | | | | | | | | Dictionary Entry[3] ...
| | | | | | | | | | Name = (nul)
| | | | | | | | | | Type = StringType
| | | | | | | | | | String Value =
"(3,*,*)"
| | | | | | | | | Dictionary Entry[4] ...
| | | | | | | | | | Name = (nul)
| | | | | | | | | | Type = StringType
| | | | | | | | | | String Value =
"(0,*,*)"
| | | | | | | | | Dictionary Entry[5] ...
| | | | | | | | | | Name = (nul)
| | | | | | | | | | Type = StringType
| | | | | | | | | | String Value =
"(1,*,*)"
| | | | | | | | | Dictionary Entry[6] ...
| | | | | | | | | | Name = (nul)
| | | | | | | | | | Type = StringType
| | | | | | | | | | String Value =
"(2,*,*)"
| | | | | | | | | Dictionary Entry[7] ...
| | | | | | | | | | Name = (nul)
| | | | | | | | | | Type = StringType
| | | | | | | | | | String Value =
"(3,*,*)"
| | | | | | | | | Dictionary Entry[8] ...
| | | | | | | | | | Name = (nul)
| | | | | | | | | | Type = StringType
| | | | | | | | | | String Value =
"(0,*,*)"
| | | | | | | | | Dictionary Entry[9] ...
| | | | | | | | | | Name = (nul)
| | | | | | | | | | Type = StringType
| | | | | | | | | | String Value =
"(1,*,*)"
| | | | | | | | | Dictionary Entry[10] ...
| | | | | | | | | | Name = (nul)
| | | | | | | | | | Type = StringType
| | | | | | | | | | String Value =
"(2,*,*)"
| | | | | | | | | Dictionary Entry[11] ...
| | | | | | | | | | Name = (nul)
| | | | | | | | | | Type = StringType
| | | | | | | | | | String Value =
"(3,*,*)"
| Dictionary Entry[70] ...
| | Name = climo_mean
| | Type = DictionaryType
| | Dict/Array Value =
| | | Nentries = 6
| | | Address = 0x2ca3970
| | | Parent = 0x8a9f88
| | | IsArray = false
| | | Dictionary Entry[0] ...
| | | | Name = file_name
| | | | Type = ArrayType
| | | | Dict/Array Value =
| | | | | Nentries = 0
| | | | | Address = 0x2990d80
| | | | | Parent = 0x2ca3970
| | | | | IsArray = true
| | | Dictionary Entry[1] ...
| | | | Name = field
| | | | Type = ArrayType
| | | | Dict/Array Value =
| | | | | Nentries = 0
| | | | | Address = 0x291fdd0
| | | | | Parent = 0x2ca3970
| | | | | IsArray = true
| | | Dictionary Entry[2] ...
| | | | Name = regrid
| | | | Type = DictionaryType
| | | | Dict/Array Value =
| | | | | Nentries = 3
| | | | | Address = 0x2866260
| | | | | Parent = 0x2ca3970
| | | | | IsArray = false
| | | | | Dictionary Entry[0] ...
| | | | | | Name = method
| | | | | | Type = IntegerType
| | | | | | Integer Value = 9
| | | | | Dictionary Entry[1] ...
| | | | | | Name = width
| | | | | | Type = IntegerType
| | | | | | Integer Value = 1
| | | | | Dictionary Entry[2] ...
| | | | | | Name = vld_thresh
| | | | | | Type = FloatType
| | | | | | Float Value = 0.5
| | | Dictionary Entry[3] ...
| | | | Name = time_interp_method
| | | | Type = IntegerType
| | | | Integer Value = 5
| | | Dictionary Entry[4] ...
| | | | Name = match_day
| | | | Type = BooleanType
| | | | Boolean Value = false
| | | Dictionary Entry[5] ...
| | | | Name = time_step
| | | | Type = IntegerType
| | | | Integer Value = 21600
| Dictionary Entry[71] ...
| | Name = ci_alpha
| | Type = ArrayType
| | Dict/Array Value =
| | | Nentries = 1
| | | Address = 0x2de5720
| | | Parent = 0x8a9f88
| | | IsArray = true
| | | Dictionary Entry[0] ...
| | | | Name = (nul)
| | | | Type = FloatType
| | | | Float Value = 0.05
| Dictionary Entry[72] ...
| | Name = boot
| | Type = DictionaryType
| | Dict/Array Value =
| | | Nentries = 5
| | | Address = 0x288c970
| | | Parent = 0x8a9f88
| | | IsArray = false
| | | Dictionary Entry[0] ...
| | | | Name = interval
| | | | Type = IntegerType
| | | | Integer Value = 1
| | | Dictionary Entry[1] ...
| | | | Name = rep_prop
| | | | Type = FloatType
| | | | Float Value = 1
| | | Dictionary Entry[2] ...
| | | | Name = n_rep
| | | | Type = IntegerType
| | | | Integer Value = 0
| | | Dictionary Entry[3] ...
| | | | Name = rng
| | | | Type = StringType
| | | | String Value = "mt19937"
| | | Dictionary Entry[4] ...
| | | | Name = seed
| | | | Type = StringType
| | | | String Value = (nul)
| Dictionary Entry[73] ...
| | Name = mask
| | Type = DictionaryType
| | Dict/Array Value =
| | | Nentries = 2
| | | Address = 0x2b10c60
| | | Parent = 0x8a9f88
| | | IsArray = false
| | | Dictionary Entry[0] ...
| | | | Name = grid
| | | | Type = StringType
| | | | String Value = (nul)
| | | Dictionary Entry[1] ...
| | | | Name = poly
| | | | Type = StringType
| | | | String Value = (nul)
| Dictionary Entry[74] ...
| | Name = block_size
| | Type = IntegerType
| | Integer Value = 100000
| Dictionary Entry[75] ...
| | Name = vld_thresh
| | Type = FloatType
| | Float Value = 1
| Dictionary Entry[76] ...
| | Name = output_stats
| | Type = DictionaryType
| | Dict/Array Value =
| | | Nentries = 12
| | | Address = 0x28b85c0
| | | Parent = 0x8a9f88
| | | IsArray = false
| | | Dictionary Entry[0] ...
| | | | Name = fho
| | | | Type = ArrayType
| | | | Dict/Array Value =
| | | | | Nentries = 3
| | | | | Address = 0x2f64430
| | | | | Parent = 0x28b85c0
| | | | | IsArray = true
| | | | | Dictionary Entry[0] ...
| | | | | | Name = (nul)
| | | | | | Type = StringType
| | | | | | String Value = "F_RATE"
| | | | | Dictionary Entry[1] ...
| | | | | | Name = (nul)
| | | | | | Type = StringType
| | | | | | String Value = "O_RATE"
| | | | | Dictionary Entry[2] ...
| | | | | | Name = (nul)
| | | | | | Type = StringType
| | | | | | String Value = "TOTAL"
| | | Dictionary Entry[1] ...
| | | | Name = ctc
| | | | Type = ArrayType
| | | | Dict/Array Value =
| | | | | Nentries = 0
| | | | | Address = 0x2e045d0
| | | | | Parent = 0x28b85c0
| | | | | IsArray = true
| | | Dictionary Entry[2] ...
| | | | Name = cts
| | | | Type = ArrayType
| | | | Dict/Array Value =
| | | | | Nentries = 3
| | | | | Address = 0x2e04ad0
| | | | | Parent = 0x28b85c0
| | | | | IsArray = true
| | | | | Dictionary Entry[0] ...
| | | | | | Name = (nul)
| | | | | | Type = StringType
| | | | | | String Value = "CSI"
| | | | | Dictionary Entry[1] ...
| | | | | | Name = (nul)
| | | | | | Type = StringType
| | | | | | String Value = "GSS"
| | | | | Dictionary Entry[2] ...
| | | | | | Name = (nul)
| | | | | | Type = StringType
| | | | | | String Value = "FAR"
| | | Dictionary Entry[3] ...
| | | | Name = mctc
| | | | Type = ArrayType
| | | | Dict/Array Value =
| | | | | Nentries = 0
| | | | | Address = 0x2e04fd0
| | | | | Parent = 0x28b85c0
| | | | | IsArray = true
| | | Dictionary Entry[4] ...
| | | | Name = mcts
| | | | Type = ArrayType
| | | | Dict/Array Value =
| | | | | Nentries = 0
| | | | | Address = 0x2e054d0
| | | | | Parent = 0x28b85c0
| | | | | IsArray = true
| | | Dictionary Entry[5] ...
| | | | Name = cnt
| | | | Type = ArrayType
| | | | Dict/Array Value =
| | | | | Nentries = 18
| | | | | Address = 0x2e059d0
| | | | | Parent = 0x28b85c0
| | | | | IsArray = true
| | | | | Dictionary Entry[0] ...
| | | | | | Name = (nul)
| | | | | | Type = StringType
| | | | | | String Value = "TOTAL"
| | | | | Dictionary Entry[1] ...
| | | | | | Name = (nul)
| | | | | | Type = StringType
| | | | | | String Value = "FBAR"
| | | | | Dictionary Entry[2] ...
| | | | | | Name = (nul)
| | | | | | Type = StringType
| | | | | | String Value = "OBAR"
| | | | | Dictionary Entry[3] ...
| | | | | | Name = (nul)
| | | | | | Type = StringType
| | | | | | String Value = "PR_CORR"
| | | | | Dictionary Entry[4] ...
| | | | | | Name = (nul)
| | | | | | Type = StringType
| | | | | | String Value = "ME"
| | | | | Dictionary Entry[5] ...
| | | | | | Name = (nul)
| | | | | | Type = StringType
| | | | | | String Value = "MAE"
| | | | | Dictionary Entry[6] ...
| | | | | | Name = (nul)
| | | | | | Type = StringType
| | | | | | String Value = "MSE"
| | | | | Dictionary Entry[7] ...
| | | | | | Name = (nul)
| | | | | | Type = StringType
| | | | | | String Value = "RMSE"
| | | | | Dictionary Entry[8] ...
| | | | | | Name = (nul)
| | | | | | Type = StringType
| | | | | | String Value = "BCMSE"
| | | | | Dictionary Entry[9] ...
| | | | | | Name = (nul)
| | | | | | Type = StringType
| | | | | | String Value = "MAD"
| | | | | Dictionary Entry[10] ...
| | | | | | Name = (nul)
| | | | | | Type = StringType
| | | | | | String Value = "E10"
| | | | | Dictionary Entry[11] ...
| | | | | | Name = (nul)
| | | | | | Type = StringType
| | | | | | String Value = "E25"
| | | | | Dictionary Entry[12] ...
| | | | | | Name = (nul)
| | | | | | Type = StringType
| | | | | | String Value = "E50"
| | | | | Dictionary Entry[13] ...
| | | | | | Name = (nul)
| | | | | | Type = StringType
| | | | | | String Value = "E75"
| | | | | Dictionary Entry[14] ...
| | | | | | Name = (nul)
| | | | | | Type = StringType
| | | | | | String Value = "E90"
| | | | | Dictionary Entry[15] ...
| | | | | | Name = (nul)
| | | | | | Type = StringType
| | | | | | String Value = "KT_CORR"
| | | | | Dictionary Entry[16] ...
| | | | | | Name = (nul)
| | | | | | Type = StringType
| | | | | | String Value = "ESTDEV"
| | | | | Dictionary Entry[17] ...
| | | | | | Name = (nul)
| | | | | | Type = StringType
| | | | | | String Value = "SP_CORR"
| | | Dictionary Entry[6] ...
| | | | Name = sl1l2
| | | | Type = ArrayType
| | | | Dict/Array Value =
| | | | | Nentries = 0
| | | | | Address = 0x2e05ed0
| | | | | Parent = 0x28b85c0
| | | | | IsArray = true
| | | Dictionary Entry[7] ...
| | | | Name = sal1l2
| | | | Type = ArrayType
| | | | Dict/Array Value =
| | | | | Nentries = 0
| | | | | Address = 0x2e063d0
| | | | | Parent = 0x28b85c0
| | | | | IsArray = true
| | | Dictionary Entry[8] ...
| | | | Name = pct
| | | | Type = ArrayType
| | | | Dict/Array Value =
| | | | | Nentries = 0
| | | | | Address = 0x2e068d0
| | | | | Parent = 0x28b85c0
| | | | | IsArray = true
| | | Dictionary Entry[9] ...
| | | | Name = pstd
| | | | Type = ArrayType
| | | | Dict/Array Value =
| | | | | Nentries = 0
| | | | | Address = 0x2e06dd0
| | | | | Parent = 0x28b85c0
| | | | | IsArray = true
| | | Dictionary Entry[10] ...
| | | | Name = pjc
| | | | Type = ArrayType
| | | | Dict/Array Value =
| | | | | Nentries = 0
| | | | | Address = 0x2f585b0
| | | | | Parent = 0x28b85c0
| | | | | IsArray = true
| | | Dictionary Entry[11] ...
| | | | Name = prc
| | | | Type = ArrayType
| | | | Dict/Array Value =
| | | | | Nentries = 0
| | | | | Address = 0x2f58ab0
| | | | | Parent = 0x28b85c0
| | | | | IsArray = true
| Dictionary Entry[77] ...
| | Name = rank_corr_flag
| | Type = BooleanType
| | Boolean Value = true
| Dictionary Entry[78] ...
| | Name = version
| | Type = StringType
| | String Value = "V6.0"
DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo object
of type "FileType_Gb1".
DEBUG 5: Parsed forecast field number 1:
VarInfoGrib::dump():
PTV = -9999
Code = -9999
LvlType = -9999
PCode = -9999
Center = -9999
Subcenter = -9999
DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo object
of type "FileType_Gb1".
DEBUG 5: Parsed forecast field number 2:
VarInfoGrib::dump():
PTV = -9999
Code = -9999
LvlType = -9999
PCode = -9999
Center = -9999
Subcenter = -9999
DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo object
of type "FileType_Gb1".
DEBUG 5: Parsed forecast field number 3:
VarInfoGrib::dump():
PTV = -9999
Code = -9999
LvlType = -9999
PCode = -9999
Center = -9999
Subcenter = -9999
DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo object
of type "FileType_Gb1".
DEBUG 5: Parsed forecast field number 4:
VarInfoGrib::dump():
PTV = -9999
Code = -9999
LvlType = -9999
PCode = -9999
Center = -9999
Subcenter = -9999
DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo object
of type "FileType_Gb1".
DEBUG 5: Parsed forecast field number 5:
VarInfoGrib::dump():
PTV = -9999
Code = -9999
LvlType = -9999
PCode = -9999
Center = -9999
Subcenter = -9999
DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo object
of type "FileType_Gb1".
DEBUG 5: Parsed forecast field number 6:
VarInfoGrib::dump():
PTV = -9999
Code = -9999
LvlType = -9999
PCode = -9999
Center = -9999
Subcenter = -9999
DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo object
of type "FileType_Gb1".
DEBUG 5: Parsed forecast field number 7:
VarInfoGrib::dump():
PTV = -9999
Code = -9999
LvlType = -9999
PCode = -9999
Center = -9999
Subcenter = -9999
DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo object
of type "FileType_Gb1".
DEBUG 5: Parsed forecast field number 8:
VarInfoGrib::dump():
PTV = -9999
Code = -9999
LvlType = -9999
PCode = -9999
Center = -9999
Subcenter = -9999
DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo object
of type "FileType_Gb1".
DEBUG 5: Parsed forecast field number 9:
VarInfoGrib::dump():
PTV = -9999
Code = -9999
LvlType = -9999
PCode = -9999
Center = -9999
Subcenter = -9999
DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo object
of type "FileType_Gb1".
DEBUG 5: Parsed forecast field number 10:
VarInfoGrib::dump():
PTV = -9999
Code = -9999
LvlType = -9999
PCode = -9999
Center = -9999
Subcenter = -9999
DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo object
of type "FileType_Gb1".
DEBUG 5: Parsed forecast field number 11:
VarInfoGrib::dump():
PTV = -9999
Code = -9999
LvlType = -9999
PCode = -9999
Center = -9999
Subcenter = -9999
DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo object
of type "FileType_Gb1".
DEBUG 5: Parsed forecast field number 12:
VarInfoGrib::dump():
PTV = -9999
Code = -9999
LvlType = -9999
PCode = -9999
Center = -9999
Subcenter = -9999
DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo object
of type "FileType_NcCF".
DEBUG 5: Parsed observation field number 1:
VarInfoNcCF::dump():
Dimension:
Nelements = 3
Nalloc = 10
AllocInc = 10
Element # 0 = 0
Element # 1 = -12345
Element # 2 = -12345
DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo object
of type "FileType_NcCF".
DEBUG 5: Parsed observation field number 2:
VarInfoNcCF::dump():
Dimension:
Nelements = 3
Nalloc = 10
AllocInc = 10
Element # 0 = 1
Element # 1 = -12345
Element # 2 = -12345
DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo object
of type "FileType_NcCF".
DEBUG 5: Parsed observation field number 3:
VarInfoNcCF::dump():
Dimension:
Nelements = 3
Nalloc = 10
AllocInc = 10
Element # 0 = 2
Element # 1 = -12345
Element # 2 = -12345
DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo object
of type "FileType_NcCF".
DEBUG 5: Parsed observation field number 4:
VarInfoNcCF::dump():
Dimension:
Nelements = 3
Nalloc = 10
AllocInc = 10
Element # 0 = 3
Element # 1 = -12345
Element # 2 = -12345
DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo object
of type "FileType_NcCF".
DEBUG 5: Parsed observation field number 5:
VarInfoNcCF::dump():
Dimension:
Nelements = 3
Nalloc = 10
AllocInc = 10
Element # 0 = 0
Element # 1 = -12345
Element # 2 = -12345
DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo object
of type "FileType_NcCF".
DEBUG 5: Parsed observation field number 6:
VarInfoNcCF::dump():
Dimension:
Nelements = 3
Nalloc = 10
AllocInc = 10
Element # 0 = 1
Element # 1 = -12345
Element # 2 = -12345
DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo object
of type "FileType_NcCF".
DEBUG 5: Parsed observation field number 7:
VarInfoNcCF::dump():
Dimension:
Nelements = 3
Nalloc = 10
AllocInc = 10
Element # 0 = 2
Element # 1 = -12345
Element # 2 = -12345
DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo object
of type "FileType_NcCF".
DEBUG 5: Parsed observation field number 8:
VarInfoNcCF::dump():
Dimension:
Nelements = 3
Nalloc = 10
AllocInc = 10
Element # 0 = 3
Element # 1 = -12345
Element # 2 = -12345
DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo object
of type "FileType_NcCF".
DEBUG 5: Parsed observation field number 9:
VarInfoNcCF::dump():
Dimension:
Nelements = 3
Nalloc = 10
AllocInc = 10
Element # 0 = 0
Element # 1 = -12345
Element # 2 = -12345
DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo object
of type "FileType_NcCF".
DEBUG 5: Parsed observation field number 10:
VarInfoNcCF::dump():
Dimension:
Nelements = 3
Nalloc = 10
AllocInc = 10
Element # 0 = 1
Element # 1 = -12345
Element # 2 = -12345
DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo object
of type "FileType_NcCF".
DEBUG 5: Parsed observation field number 11:
VarInfoNcCF::dump():
Dimension:
Nelements = 3
Nalloc = 10
AllocInc = 10
Element # 0 = 2
Element # 1 = -12345
Element # 2 = -12345
DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo object
of type "FileType_NcCF".
DEBUG 5: Parsed observation field number 12:
VarInfoNcCF::dump():
Dimension:
Nelements = 3
Nalloc = 10
AllocInc = 10
Element # 0 = 3
Element # 1 = -12345
Element # 2 = -12345
DEBUG 5: Parsed forecast categorical thresholds: NA
DEBUG 5: Parsed observed categorical thresholds: NA
DEBUG 5: Parsed forecast continuous thresholds: NA
DEBUG 5: Parsed observed continuous thresholds: NA
DEBUG 5: Parsed continuous threshold logic: INTERSECTION
DEBUG 3: Use the matching forecast and observation grids.
DEBUG 3: Grid Definition: Projection: Lat/Lon Nx: 360 Ny: 180 lat_ll:
-89.500 lon_ll: 179.500 delta_lat: 1.000 delta_lon: 1.000
DEBUG 2: Processing masking regions.
DEBUG 1: Length of configuration "fcst.field" = 12
DEBUG 1: Length of configuration "obs.field" = 12
DEBUG 1: Length of forecast file list = 12
DEBUG 1: Length of observation file list = 12
DEBUG 1: Series defined by the "fcst.field" configuration entry of
length 12.
DEBUG 2: Computing statistics using a block size of 100000, requiring
1 pass(es) through the 360 x 180 grid.
DEBUG 2: Processing data pass number 1 of 1 for grid points 1 to
64800.
DEBUG 2: Processing series entry 1 of 12: AOD/L0 versus
total_aod_mean(0,*,*)
DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() -> created new
Met2dDataFile object of type "FileType_Gb1".
DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found range match for
VarInfo "AOD/L0" in GRIB record 1 of GRIB file
"/global/noscrub/Partha.Bhattacharjee/MET_Test/Fcst/20150628/AOD550.t00z.f00.grb".
DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 1 GRIB records
matching VarInfo "AOD/L0" in GRIB file
"/global/noscrub/Partha.Bhattacharjee/MET_Test/Fcst/20150628/AOD550.t00z.f00.grb".
DEBUG 2: Found data for AOD/L0 in file:
/global/noscrub/Partha.Bhattacharjee/MET_Test/Fcst/20150628/AOD550.t00z.f00.grb
DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() -> created new
Met2dDataFile object of type "FileType_NcCF".
DEBUG 4: NcCfFile::open() -> parsing units for the time variable
"seconds since 2015-06-28 00:00:00 UTC"
DEBUG 4: parse_cf_time_string() -> parsed NetCDF CF convention time
unit string "seconds since 2015-06-28 00:00:00 UTC" as a reference
time of 20150628_000000 and 1 second(s) per time step.
DEBUG 4: NcCfFile::open() -> could not extract init time from the
"forecast_reference_time" variable.
DEBUG 4: NcCfFile::open() -> could not extract init time from file
name.
DEBUG 4: NcCfFile::getData() -> setting the unset init time to the
valid time of 20150628_000000.
DEBUG 4:
DEBUG 4: Data plane information:
DEBUG 4: plane min: 0.0089742
DEBUG 4: plane max: 6.77107
DEBUG 4: valid time: 20150628_000000
DEBUG 4: lead time: 000000
DEBUG 4: init time: 20150628_000000
DEBUG 4: accum time: 000000
DEBUG 2: Found data for total_aod_mean(0,*,*) in file:
/global/noscrub/Partha.Bhattacharjee/MET_Test/Obs/icap_2015062800_MME_modeaod550.nc
DEBUG 3: Found 0 climatology mean field(s) for forecast AOD/L0.
------------------------------------------------
Subject: Multiple records of AOD in a grib file
From: John Halley Gotway
Time: Wed Jan 17 15:43:28 2018
Partha,
Unfortunately, I don't have access to WCOSS, but my colleague, Julie
Prestopnik, does. I'll ask her to log on and try to replicate the
behavior
you describe. And then we can debug it.
Thanks,
John
On Wed, Jan 17, 2018 at 2:34 PM, Partha Bhattacharjee - NOAA Affiliate
via
RT <met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=83618 >
>
> Thanks John. Related to my earlier email, I am getting error message
when
> trying to use MET Series Stat. I am trying to compare regridded
forecast
> file and observational file but getting the following error :
>
> GSL_RNG_TYPE=mt19937
> GSL_RNG_SEED=18446744073587810138
> terminate called after throwing an instance of 'std::logic_error'
> what(): basic_string::_S_construct NULL not valid
>
> I have attached the output.txt along with this email. My files are
in WCOSS
> at : /global/noscrub/Partha.Bhattacharjee/MET_Test. Original Grib
files
> are
> in the Fcst directory : "AOD550.Grb.t00z.f00" and interpolated ones
are
> ones with "grb" at the end. I have used copygb to interpolate these
input
> files which are at 384 X 190 to 360 X 180 to make them same as Obs
files
> are. But not sure if the error is due to not matching grids or from
> something else.
>
> Thanks.
> Partha
>
> On Wed, Jan 17, 2018 at 4:26 PM, John Halley Gotway via RT <
> met_help at ucar.edu> wrote:
>
> > Partha,
> >
> > Ah, OK. So your file contains 6 GRIB records, each of which has a
> > different GRIB code (129, 133, 134, 135, 76, and 137). Those are
the
> kpds5
> > values listed in the output of wgrib. Running wgrib -V reveals
that 5 of
> > the 6 records are encoded as:
> > center 7 subcenter 0 process 82 Table 141 scan: WE:SN
winds(N/S)
> >
> > So this is NCEP GRIB1 table number 141, and that can be found
here:
> > met-6.1/share/met/table_files/grib1_ncep_141_7.txt
> >
> > Here's a selection of lines from that file:
> >
> > 129 141 7 -1 "AOD" "Aerosol Optical Depth" ""
> > 133 141 7 -1 "var133" "undefined" ""
> > 134 141 7 -1 "var134" "undefined" ""
> > 135 141 7 -1 "var135" "undefined" ""
> > 137 141 7 -1 "var137" "undefined" ""
> >
> > And this info is taken from this GRIB table:
> >
http://www.nco.ncep.noaa.gov/pmb/docs/on388/table2.html#TABLE141
> >
> > So you could just refer to this data using those "var133" (and so
on)
> > string:
> >
> > plot_data_plane out.grb plot_rec1.ps 'name="var133";
level="L0";'
> >
> > Alternatively, you could specifying the "GRIB1_code" value instead
of the
> > name:
> >
> > plot_data_plane out.grb plot_rec1.ps 'GRIB1_code=133;
level="L0";'
> >
> > Both approaches would result in the same plot. Presumably, since
you
> have
> > data for codes 133, 134, 135, and 137, perhaps these codes are no
longer
> > "reserved for future use", as indicated on the NCO website.
> >
> > Hope that helps clarify.
> >
> > Thanks,
> > John
> >
> > On Wed, Jan 17, 2018 at 12:05 PM, Partha Bhattacharjee - NOAA
Affiliate
> via
> > RT <met_help at ucar.edu> wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=83618 >
> > >
> > > John,
> > >
> > > That is strange. When I did wgrib to out.grib I am getting the
> following
> > :
> > >
> > > wgrib out.grb
> > > 1:0:d=15062500:AOD:kpds5=129:kpds6=200:kpds7=0:TR=10:P1=0:
> > > P2=0:TimeU=1:atmos
> > > col:anl:NAve=0
> > > 2:129684:d=15062500:var133:kpds5=133:kpds6=200:kpds7=0:
> > > TR=10:P1=0:P2=0:TimeU=1:atmos
> > > col:anl:NAve=0
> > > 3:259368:d=15062500:var134:kpds5=134:kpds6=200:kpds7=0:
> > > TR=10:P1=0:P2=0:TimeU=1:atmos
> > > col:anl:NAve=0
> > > 4:389052:d=15062500:var135:kpds5=135:kpds6=200:kpds7=0:
> > > TR=10:P1=0:P2=0:TimeU=1:atmos
> > > col:anl:NAve=0
> > > 5:518736:d=15062500:CWAT:kpds5=76:kpds6=200:kpds7=0:TR=
> > > 10:P1=0:P2=0:TimeU=1:atmos
> > > col:anl:NAve=0
> > > 6:648420:d=15062500:var137:kpds5=137:kpds6=200:kpds7=0:
> > > TR=10:P1=0:P2=0:TimeU=1:atmos
> > > col:anl:NAve=0
> > >
> > > and I am running in wcoss and wgrib points to
> > > which wgrib
> > > alias wgrib='/nwprod/util/exec/wgrib'
> > > /nwprod/util/exec/wgrib
> > >
> > > So, I used "AOD" in the first record to plot like :
> > >
> > > plot_data_plane out.grb aod.ps 'name="AOD"; level="L0";' -v 6
> > > DEBUG 1: Opening data file: out.grb
> > > DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() -> created
new
> > > Met2dDataFile object of type "FileType_Gb1".
> > > DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo
object
> of
> > > type "FileType_Gb1".
> > > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found range
match for
> > > VarInfo "AOD/L0" in GRIB record 1 of GRIB file "out.grb".
> > > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 1 GRIB
records
> > > matching VarInfo "AOD/L0" in GRIB file "out.grb".
> > > DEBUG 1: Creating postscript file: aod.ps
> > > DEBUG 4: draw_map() -> Plotting map data file:
> > > /awc/noscrub/Julie.Prestopnik/met/6.0/share/met/map/country_data
> > > DEBUG 4: draw_map() -> Plotting map data file:
> > >
/awc/noscrub/Julie.Prestopnik/met/6.0/share/met/map/usa_state_data
> > > DEBUG 4: draw_map() -> Plotting map data file:
> > >
/awc/noscrub/Julie.Prestopnik/met/6.0/share/met/map/major_lakes_data
> > >
> > > I do not think I have $MET_GRIB_TABLES environment variable set.
> > >
> > > Partha
> > >
> > > On Wed, Jan 17, 2018 at 1:33 PM, John Halley Gotway via RT <
> > > met_help at ucar.edu> wrote:
> > >
> > > > Partha,
> > > >
> > > > I'm confused, when I run wgrib on the data you sent, I don't
see any
> > > > mention of AOD. But I do get a warning message about an
undefined
> > > > parameter table:
> > > >
> > > > [johnhg at number5]% wgrib out.grb
> > > >
> > > > Undefined parameter table (center 7-0 table 141), using NCEP-
opn
> > > > 1:0:d=15062500:MSLMA:kpds5=129:kpds6=200:kpds7=0:TR=10:
> > > > P1=0:P2=0:TimeU=1:atmos
> > > > col:anl:NAve=0
> > > > 2:129684:d=15062500:KX:kpds5=133:kpds6=200:kpds7=0:TR=10:
> > > > P1=0:P2=0:TimeU=1:atmos
> > > > col:anl:NAve=0
> > > > 3:259368:d=15062500:SX:kpds5=134:kpds6=200:kpds7=0:TR=10:
> > > > P1=0:P2=0:TimeU=1:atmos
> > > > col:anl:NAve=0
> > > > 4:389052:d=15062500:MCONV:kpds5=135:kpds6=200:kpds7=0:
> > > > TR=10:P1=0:P2=0:TimeU=1:atmos
> > > > col:anl:NAve=0
> > > > 5:518736:d=15062500:CWAT:kpds5=76:kpds6=200:kpds7=0:TR=
> > > > 10:P1=0:P2=0:TimeU=1:atmos
> > > > col:anl:NAve=0
> > > > 6:648420:d=15062500:TSLSA:kpds5=137:kpds6=200:kpds7=0:
> > > > TR=10:P1=0:P2=0:TimeU=1:atmos
> > > > col:anl:NAve=0
> > > >
> > > > All of these 6 records appear to have different GRIB codes.
Based on
> > the
> > > > question you asked, I was expecting to find multiple records
with the
> > > > *same* GRIB code but were different in some other way (i.e.
different
> > > level
> > > > type or different statistical processing type).
> > > >
> > > > You say these are all AOD, but the plot_data_plane command you
sent
> > does
> > > > *not* work for me:
> > > >
> > > > plot_data_plane out.grb aod.ps 'name="AOD"; level="A0";'
> > > > ERROR : plot_data_plane -> trouble getting field "name="AOD";
> > > level="A0";"
> > > > from file "out.grb"
> > > >
> > > > Please check to see if you have the $MET_GRIB_TABLES
environment
> > variable
> > > > set. If so, please send me the table(s) to which it points.
> > > >
> > > > Thanks,
> > > > John
> > > >
> > > >
> > > > On Wed, Jan 17, 2018 at 10:35 AM, Partha Bhattacharjee - NOAA
> Affiliate
> > > via
> > > > RT <met_help at ucar.edu> wrote:
> > > >
> > > > >
> > > > > Wed Jan 17 10:35:05 2018: Request 83618 was acted upon.
> > > > > Transaction: Ticket created by partha.bhattacharjee at noaa.gov
> > > > > Queue: met_help
> > > > > Subject: Multiple records of AOD in a grib file
> > > > > Owner: Nobody
> > > > > Requestors: partha.bhattacharjee at noaa.gov
> > > > > Status: new
> > > > > Ticket <URL: https://rt.rap.ucar.edu/rt/
> > Ticket/Display.html?id=83618
> > > >
> > > > >
> > > > >
> > > > > Hi
> > > > >
> > > > > I have an input GRIB1 file with multiple records of
AOD,before I
> used
> > > > > copygb to interpolate that into (-90:90;-180:180} grid to
use for
> > > > > Grid-Stat. I am wondering how to read multiple AOD records
in that
> > > file,
> > > > is
> > > > > there any Grib table specific number or template I can use ?
> > > > >
> > > > > MET able to read the interpolated file record when I use
this
> > command :
> > > > > plot_data_plane out.grb aod.ps 'name="AOD"; level="A0";'
and
> > > generated
> > > > > the
> > > > > attached figure.
> > > > >
> > > > > Attached are original and interpolated grib1 files.
> > > > >
> > > > > Thanks.
> > > > > Partha
> > > > >
> > > > >
> > > >
> > > >
> > >
> > >
> > > --
> > > Partha S. Bhattacharjee, Ph.D.
> > > I.M. Systems Group
> > > Environmental Modeling Centre
> > > NCWCP
> > > 5830 University Research Court.
> > > College Park, MD 20740
> > >
> > >
> >
> >
>
>
> --
> Partha S. Bhattacharjee, Ph.D.
> I.M. Systems Group
> Environmental Modeling Centre
> NCWCP
> 5830 University Research Court.
> College Park, MD 20740
>
>
------------------------------------------------
Subject: Multiple records of AOD in a grib file
From: Partha Bhattacharjee - NOAA Affiliate
Time: Wed Jan 17 17:07:16 2018
Thanks John. I realized WCOSS will be out next 2 days, so I copied the
files over Theia
at /scratch4/NCEPDEV/global/noscrub/Partha.Bhattacharjee/MET_Test.
Partha
On Wed, Jan 17, 2018 at 5:43 PM, John Halley Gotway via RT <
met_help at ucar.edu> wrote:
> Partha,
>
> Unfortunately, I don't have access to WCOSS, but my colleague, Julie
> Prestopnik, does. I'll ask her to log on and try to replicate the
behavior
> you describe. And then we can debug it.
>
> Thanks,
> John
>
> On Wed, Jan 17, 2018 at 2:34 PM, Partha Bhattacharjee - NOAA
Affiliate via
> RT <met_help at ucar.edu> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=83618 >
> >
> > Thanks John. Related to my earlier email, I am getting error
message when
> > trying to use MET Series Stat. I am trying to compare regridded
forecast
> > file and observational file but getting the following error :
> >
> > GSL_RNG_TYPE=mt19937
> > GSL_RNG_SEED=18446744073587810138
> > terminate called after throwing an instance of 'std::logic_error'
> > what(): basic_string::_S_construct NULL not valid
> >
> > I have attached the output.txt along with this email. My files are
in
> WCOSS
> > at : /global/noscrub/Partha.Bhattacharjee/MET_Test. Original Grib
files
> > are
> > in the Fcst directory : "AOD550.Grb.t00z.f00" and interpolated
ones are
> > ones with "grb" at the end. I have used copygb to interpolate
these input
> > files which are at 384 X 190 to 360 X 180 to make them same as Obs
files
> > are. But not sure if the error is due to not matching grids or
from
> > something else.
> >
> > Thanks.
> > Partha
> >
> > On Wed, Jan 17, 2018 at 4:26 PM, John Halley Gotway via RT <
> > met_help at ucar.edu> wrote:
> >
> > > Partha,
> > >
> > > Ah, OK. So your file contains 6 GRIB records, each of which has
a
> > > different GRIB code (129, 133, 134, 135, 76, and 137). Those
are the
> > kpds5
> > > values listed in the output of wgrib. Running wgrib -V reveals
that 5
> of
> > > the 6 records are encoded as:
> > > center 7 subcenter 0 process 82 Table 141 scan: WE:SN
winds(N/S)
> > >
> > > So this is NCEP GRIB1 table number 141, and that can be found
here:
> > > met-6.1/share/met/table_files/grib1_ncep_141_7.txt
> > >
> > > Here's a selection of lines from that file:
> > >
> > > 129 141 7 -1 "AOD" "Aerosol Optical Depth" ""
> > > 133 141 7 -1 "var133" "undefined" ""
> > > 134 141 7 -1 "var134" "undefined" ""
> > > 135 141 7 -1 "var135" "undefined" ""
> > > 137 141 7 -1 "var137" "undefined" ""
> > >
> > > And this info is taken from this GRIB table:
> > >
http://www.nco.ncep.noaa.gov/pmb/docs/on388/table2.html#TABLE141
> > >
> > > So you could just refer to this data using those "var133" (and
so on)
> > > string:
> > >
> > > plot_data_plane out.grb plot_rec1.ps 'name="var133";
level="L0";'
> > >
> > > Alternatively, you could specifying the "GRIB1_code" value
instead of
> the
> > > name:
> > >
> > > plot_data_plane out.grb plot_rec1.ps 'GRIB1_code=133;
level="L0";'
> > >
> > > Both approaches would result in the same plot. Presumably,
since you
> > have
> > > data for codes 133, 134, 135, and 137, perhaps these codes are
no
> longer
> > > "reserved for future use", as indicated on the NCO website.
> > >
> > > Hope that helps clarify.
> > >
> > > Thanks,
> > > John
> > >
> > > On Wed, Jan 17, 2018 at 12:05 PM, Partha Bhattacharjee - NOAA
Affiliate
> > via
> > > RT <met_help at ucar.edu> wrote:
> > >
> > > >
> > > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=83618
>
> > > >
> > > > John,
> > > >
> > > > That is strange. When I did wgrib to out.grib I am getting the
> > following
> > > :
> > > >
> > > > wgrib out.grb
> > > > 1:0:d=15062500:AOD:kpds5=129:kpds6=200:kpds7=0:TR=10:P1=0:
> > > > P2=0:TimeU=1:atmos
> > > > col:anl:NAve=0
> > > > 2:129684:d=15062500:var133:kpds5=133:kpds6=200:kpds7=0:
> > > > TR=10:P1=0:P2=0:TimeU=1:atmos
> > > > col:anl:NAve=0
> > > > 3:259368:d=15062500:var134:kpds5=134:kpds6=200:kpds7=0:
> > > > TR=10:P1=0:P2=0:TimeU=1:atmos
> > > > col:anl:NAve=0
> > > > 4:389052:d=15062500:var135:kpds5=135:kpds6=200:kpds7=0:
> > > > TR=10:P1=0:P2=0:TimeU=1:atmos
> > > > col:anl:NAve=0
> > > > 5:518736:d=15062500:CWAT:kpds5=76:kpds6=200:kpds7=0:TR=
> > > > 10:P1=0:P2=0:TimeU=1:atmos
> > > > col:anl:NAve=0
> > > > 6:648420:d=15062500:var137:kpds5=137:kpds6=200:kpds7=0:
> > > > TR=10:P1=0:P2=0:TimeU=1:atmos
> > > > col:anl:NAve=0
> > > >
> > > > and I am running in wcoss and wgrib points to
> > > > which wgrib
> > > > alias wgrib='/nwprod/util/exec/wgrib'
> > > > /nwprod/util/exec/wgrib
> > > >
> > > > So, I used "AOD" in the first record to plot like :
> > > >
> > > > plot_data_plane out.grb aod.ps 'name="AOD"; level="L0";' -v 6
> > > > DEBUG 1: Opening data file: out.grb
> > > > DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() ->
created new
> > > > Met2dDataFile object of type "FileType_Gb1".
> > > > DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo
object
> > of
> > > > type "FileType_Gb1".
> > > > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found range
match
> for
> > > > VarInfo "AOD/L0" in GRIB record 1 of GRIB file "out.grb".
> > > > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 1 GRIB
> records
> > > > matching VarInfo "AOD/L0" in GRIB file "out.grb".
> > > > DEBUG 1: Creating postscript file: aod.ps
> > > > DEBUG 4: draw_map() -> Plotting map data file:
> > > >
/awc/noscrub/Julie.Prestopnik/met/6.0/share/met/map/country_data
> > > > DEBUG 4: draw_map() -> Plotting map data file:
> > > >
/awc/noscrub/Julie.Prestopnik/met/6.0/share/met/map/usa_state_data
> > > > DEBUG 4: draw_map() -> Plotting map data file:
> > > >
/awc/noscrub/Julie.Prestopnik/met/6.0/share/met/map/major_lakes_data
> > > >
> > > > I do not think I have $MET_GRIB_TABLES environment variable
set.
> > > >
> > > > Partha
> > > >
> > > > On Wed, Jan 17, 2018 at 1:33 PM, John Halley Gotway via RT <
> > > > met_help at ucar.edu> wrote:
> > > >
> > > > > Partha,
> > > > >
> > > > > I'm confused, when I run wgrib on the data you sent, I don't
see
> any
> > > > > mention of AOD. But I do get a warning message about an
undefined
> > > > > parameter table:
> > > > >
> > > > > [johnhg at number5]% wgrib out.grb
> > > > >
> > > > > Undefined parameter table (center 7-0 table 141), using
NCEP-opn
> > > > > 1:0:d=15062500:MSLMA:kpds5=129:kpds6=200:kpds7=0:TR=10:
> > > > > P1=0:P2=0:TimeU=1:atmos
> > > > > col:anl:NAve=0
> > > > > 2:129684:d=15062500:KX:kpds5=133:kpds6=200:kpds7=0:TR=10:
> > > > > P1=0:P2=0:TimeU=1:atmos
> > > > > col:anl:NAve=0
> > > > > 3:259368:d=15062500:SX:kpds5=134:kpds6=200:kpds7=0:TR=10:
> > > > > P1=0:P2=0:TimeU=1:atmos
> > > > > col:anl:NAve=0
> > > > > 4:389052:d=15062500:MCONV:kpds5=135:kpds6=200:kpds7=0:
> > > > > TR=10:P1=0:P2=0:TimeU=1:atmos
> > > > > col:anl:NAve=0
> > > > > 5:518736:d=15062500:CWAT:kpds5=76:kpds6=200:kpds7=0:TR=
> > > > > 10:P1=0:P2=0:TimeU=1:atmos
> > > > > col:anl:NAve=0
> > > > > 6:648420:d=15062500:TSLSA:kpds5=137:kpds6=200:kpds7=0:
> > > > > TR=10:P1=0:P2=0:TimeU=1:atmos
> > > > > col:anl:NAve=0
> > > > >
> > > > > All of these 6 records appear to have different GRIB codes.
Based
> on
> > > the
> > > > > question you asked, I was expecting to find multiple records
with
> the
> > > > > *same* GRIB code but were different in some other way (i.e.
> different
> > > > level
> > > > > type or different statistical processing type).
> > > > >
> > > > > You say these are all AOD, but the plot_data_plane command
you sent
> > > does
> > > > > *not* work for me:
> > > > >
> > > > > plot_data_plane out.grb aod.ps 'name="AOD"; level="A0";'
> > > > > ERROR : plot_data_plane -> trouble getting field
"name="AOD";
> > > > level="A0";"
> > > > > from file "out.grb"
> > > > >
> > > > > Please check to see if you have the $MET_GRIB_TABLES
environment
> > > variable
> > > > > set. If so, please send me the table(s) to which it points.
> > > > >
> > > > > Thanks,
> > > > > John
> > > > >
> > > > >
> > > > > On Wed, Jan 17, 2018 at 10:35 AM, Partha Bhattacharjee -
NOAA
> > Affiliate
> > > > via
> > > > > RT <met_help at ucar.edu> wrote:
> > > > >
> > > > > >
> > > > > > Wed Jan 17 10:35:05 2018: Request 83618 was acted upon.
> > > > > > Transaction: Ticket created by
partha.bhattacharjee at noaa.gov
> > > > > > Queue: met_help
> > > > > > Subject: Multiple records of AOD in a grib file
> > > > > > Owner: Nobody
> > > > > > Requestors: partha.bhattacharjee at noaa.gov
> > > > > > Status: new
> > > > > > Ticket <URL: https://rt.rap.ucar.edu/rt/
> > > Ticket/Display.html?id=83618
> > > > >
> > > > > >
> > > > > >
> > > > > > Hi
> > > > > >
> > > > > > I have an input GRIB1 file with multiple records of
AOD,before I
> > used
> > > > > > copygb to interpolate that into (-90:90;-180:180} grid to
use for
> > > > > > Grid-Stat. I am wondering how to read multiple AOD records
in
> that
> > > > file,
> > > > > is
> > > > > > there any Grib table specific number or template I can use
?
> > > > > >
> > > > > > MET able to read the interpolated file record when I use
this
> > > command :
> > > > > > plot_data_plane out.grb aod.ps 'name="AOD"; level="A0";'
and
> > > > generated
> > > > > > the
> > > > > > attached figure.
> > > > > >
> > > > > > Attached are original and interpolated grib1 files.
> > > > > >
> > > > > > Thanks.
> > > > > > Partha
> > > > > >
> > > > > >
> > > > >
> > > > >
> > > >
> > > >
> > > > --
> > > > Partha S. Bhattacharjee, Ph.D.
> > > > I.M. Systems Group
> > > > Environmental Modeling Centre
> > > > NCWCP
> > > > 5830 University Research Court.
> > > > College Park, MD 20740
> > > >
> > > >
> > >
> > >
> >
> >
> > --
> > Partha S. Bhattacharjee, Ph.D.
> > I.M. Systems Group
> > Environmental Modeling Centre
> > NCWCP
> > 5830 University Research Court.
> > College Park, MD 20740
> >
> >
>
>
--
Partha S. Bhattacharjee, Ph.D.
I.M. Systems Group
Environmental Modeling Centre
NCWCP
5830 University Research Court.
College Park, MD 20740
------------------------------------------------
Subject: Multiple records of AOD in a grib file
From: Julie Prestopnik
Time: Thu Jan 18 14:35:04 2018
Hi Partha. Thank you for remembering the WCOSS outage and for coyping
the
files over to theia! We were able to reproduce the problem you are
experiencing.
This error is occurring because of the empty string for the units for
AOD
in the table file grib1_ncep_141_7.txt. Thank you for bringing this
to our
attention. We have added a ticket to resolve this issue. In the
meantime,
you can work around this issue by copying over the file
grib1_ncep_141_7.txt from your met-6.0/share/met/table_files/
directory
into your current directory. Once you do that, please modify the
line:
129 141 7 -1 "AOD" "Aerosol Optical Depth" ""
changing the empty string at the end to be "NA" or whatever you'd
like.
For example,
129 141 7 -1 "AOD" "Aerosol Optical Depth" "NA"
Then, please set the environment variable MET_GRIB_TABLES to point to
the
location of this new table. For example:
export MET_GRIB_TABLES= /scratch4/NCEPDEV/global/
noscrub/Partha.Bhattacharjee/MET_Test/grib1_ncep_141_7.txt
And, then please try running your script again. We were able to get a
successful run once we made these changes. Please let us know if you
encounter any problems or have any further questions.
Thanks,
Julie
On Wed, Jan 17, 2018 at 5:07 PM, Partha Bhattacharjee - NOAA Affiliate
via
RT <met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=83618 >
>
> Thanks John. I realized WCOSS will be out next 2 days, so I copied
the
> files over Theia
> at /scratch4/NCEPDEV/global/noscrub/Partha.Bhattacharjee/MET_Test.
>
> Partha
>
> On Wed, Jan 17, 2018 at 5:43 PM, John Halley Gotway via RT <
> met_help at ucar.edu> wrote:
>
> > Partha,
> >
> > Unfortunately, I don't have access to WCOSS, but my colleague,
Julie
> > Prestopnik, does. I'll ask her to log on and try to replicate the
> behavior
> > you describe. And then we can debug it.
> >
> > Thanks,
> > John
> >
> > On Wed, Jan 17, 2018 at 2:34 PM, Partha Bhattacharjee - NOAA
Affiliate
> via
> > RT <met_help at ucar.edu> wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=83618 >
> > >
> > > Thanks John. Related to my earlier email, I am getting error
message
> when
> > > trying to use MET Series Stat. I am trying to compare regridded
> forecast
> > > file and observational file but getting the following error :
> > >
> > > GSL_RNG_TYPE=mt19937
> > > GSL_RNG_SEED=18446744073587810138
> > > terminate called after throwing an instance of
'std::logic_error'
> > > what(): basic_string::_S_construct NULL not valid
> > >
> > > I have attached the output.txt along with this email. My files
are in
> > WCOSS
> > > at : /global/noscrub/Partha.Bhattacharjee/MET_Test. Original
Grib
> files
> > > are
> > > in the Fcst directory : "AOD550.Grb.t00z.f00" and interpolated
ones are
> > > ones with "grb" at the end. I have used copygb to interpolate
these
> input
> > > files which are at 384 X 190 to 360 X 180 to make them same as
Obs
> files
> > > are. But not sure if the error is due to not matching grids or
from
> > > something else.
> > >
> > > Thanks.
> > > Partha
> > >
> > > On Wed, Jan 17, 2018 at 4:26 PM, John Halley Gotway via RT <
> > > met_help at ucar.edu> wrote:
> > >
> > > > Partha,
> > > >
> > > > Ah, OK. So your file contains 6 GRIB records, each of which
has a
> > > > different GRIB code (129, 133, 134, 135, 76, and 137). Those
are the
> > > kpds5
> > > > values listed in the output of wgrib. Running wgrib -V
reveals that
> 5
> > of
> > > > the 6 records are encoded as:
> > > > center 7 subcenter 0 process 82 Table 141 scan: WE:SN
winds(N/S)
> > > >
> > > > So this is NCEP GRIB1 table number 141, and that can be found
here:
> > > > met-6.1/share/met/table_files/grib1_ncep_141_7.txt
> > > >
> > > > Here's a selection of lines from that file:
> > > >
> > > > 129 141 7 -1 "AOD" "Aerosol Optical Depth" ""
> > > > 133 141 7 -1 "var133" "undefined" ""
> > > > 134 141 7 -1 "var134" "undefined" ""
> > > > 135 141 7 -1 "var135" "undefined" ""
> > > > 137 141 7 -1 "var137" "undefined" ""
> > > >
> > > > And this info is taken from this GRIB table:
> > > >
http://www.nco.ncep.noaa.gov/pmb/docs/on388/table2.html#TABLE141
> > > >
> > > > So you could just refer to this data using those "var133" (and
so on)
> > > > string:
> > > >
> > > > plot_data_plane out.grb plot_rec1.ps 'name="var133";
> level="L0";'
> > > >
> > > > Alternatively, you could specifying the "GRIB1_code" value
instead of
> > the
> > > > name:
> > > >
> > > > plot_data_plane out.grb plot_rec1.ps 'GRIB1_code=133;
> level="L0";'
> > > >
> > > > Both approaches would result in the same plot. Presumably,
since you
> > > have
> > > > data for codes 133, 134, 135, and 137, perhaps these codes are
no
> > longer
> > > > "reserved for future use", as indicated on the NCO website.
> > > >
> > > > Hope that helps clarify.
> > > >
> > > > Thanks,
> > > > John
> > > >
> > > > On Wed, Jan 17, 2018 at 12:05 PM, Partha Bhattacharjee - NOAA
> Affiliate
> > > via
> > > > RT <met_help at ucar.edu> wrote:
> > > >
> > > > >
> > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=83618 >
> > > > >
> > > > > John,
> > > > >
> > > > > That is strange. When I did wgrib to out.grib I am getting
the
> > > following
> > > > :
> > > > >
> > > > > wgrib out.grb
> > > > > 1:0:d=15062500:AOD:kpds5=129:kpds6=200:kpds7=0:TR=10:P1=0:
> > > > > P2=0:TimeU=1:atmos
> > > > > col:anl:NAve=0
> > > > > 2:129684:d=15062500:var133:kpds5=133:kpds6=200:kpds7=0:
> > > > > TR=10:P1=0:P2=0:TimeU=1:atmos
> > > > > col:anl:NAve=0
> > > > > 3:259368:d=15062500:var134:kpds5=134:kpds6=200:kpds7=0:
> > > > > TR=10:P1=0:P2=0:TimeU=1:atmos
> > > > > col:anl:NAve=0
> > > > > 4:389052:d=15062500:var135:kpds5=135:kpds6=200:kpds7=0:
> > > > > TR=10:P1=0:P2=0:TimeU=1:atmos
> > > > > col:anl:NAve=0
> > > > > 5:518736:d=15062500:CWAT:kpds5=76:kpds6=200:kpds7=0:TR=
> > > > > 10:P1=0:P2=0:TimeU=1:atmos
> > > > > col:anl:NAve=0
> > > > > 6:648420:d=15062500:var137:kpds5=137:kpds6=200:kpds7=0:
> > > > > TR=10:P1=0:P2=0:TimeU=1:atmos
> > > > > col:anl:NAve=0
> > > > >
> > > > > and I am running in wcoss and wgrib points to
> > > > > which wgrib
> > > > > alias wgrib='/nwprod/util/exec/wgrib'
> > > > > /nwprod/util/exec/wgrib
> > > > >
> > > > > So, I used "AOD" in the first record to plot like :
> > > > >
> > > > > plot_data_plane out.grb aod.ps 'name="AOD"; level="L0";' -v
6
> > > > > DEBUG 1: Opening data file: out.grb
> > > > > DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() ->
created
> new
> > > > > Met2dDataFile object of type "FileType_Gb1".
> > > > > DEBUG 4: VarInfoFactory::new_var_info() -> created new
VarInfo
> object
> > > of
> > > > > type "FileType_Gb1".
> > > > > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found range
match
> > for
> > > > > VarInfo "AOD/L0" in GRIB record 1 of GRIB file "out.grb".
> > > > > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 1
GRIB
> > records
> > > > > matching VarInfo "AOD/L0" in GRIB file "out.grb".
> > > > > DEBUG 1: Creating postscript file: aod.ps
> > > > > DEBUG 4: draw_map() -> Plotting map data file:
> > > > >
/awc/noscrub/Julie.Prestopnik/met/6.0/share/met/map/country_data
> > > > > DEBUG 4: draw_map() -> Plotting map data file:
> > > > >
/awc/noscrub/Julie.Prestopnik/met/6.0/share/met/map/usa_state_data
> > > > > DEBUG 4: draw_map() -> Plotting map data file:
> > > > > /awc/noscrub/Julie.Prestopnik/met/6.0/share/met/map/major_
> lakes_data
> > > > >
> > > > > I do not think I have $MET_GRIB_TABLES environment variable
set.
> > > > >
> > > > > Partha
> > > > >
> > > > > On Wed, Jan 17, 2018 at 1:33 PM, John Halley Gotway via RT <
> > > > > met_help at ucar.edu> wrote:
> > > > >
> > > > > > Partha,
> > > > > >
> > > > > > I'm confused, when I run wgrib on the data you sent, I
don't see
> > any
> > > > > > mention of AOD. But I do get a warning message about an
> undefined
> > > > > > parameter table:
> > > > > >
> > > > > > [johnhg at number5]% wgrib out.grb
> > > > > >
> > > > > > Undefined parameter table (center 7-0 table 141), using
NCEP-opn
> > > > > > 1:0:d=15062500:MSLMA:kpds5=129:kpds6=200:kpds7=0:TR=10:
> > > > > > P1=0:P2=0:TimeU=1:atmos
> > > > > > col:anl:NAve=0
> > > > > > 2:129684:d=15062500:KX:kpds5=133:kpds6=200:kpds7=0:TR=10:
> > > > > > P1=0:P2=0:TimeU=1:atmos
> > > > > > col:anl:NAve=0
> > > > > > 3:259368:d=15062500:SX:kpds5=134:kpds6=200:kpds7=0:TR=10:
> > > > > > P1=0:P2=0:TimeU=1:atmos
> > > > > > col:anl:NAve=0
> > > > > > 4:389052:d=15062500:MCONV:kpds5=135:kpds6=200:kpds7=0:
> > > > > > TR=10:P1=0:P2=0:TimeU=1:atmos
> > > > > > col:anl:NAve=0
> > > > > > 5:518736:d=15062500:CWAT:kpds5=76:kpds6=200:kpds7=0:TR=
> > > > > > 10:P1=0:P2=0:TimeU=1:atmos
> > > > > > col:anl:NAve=0
> > > > > > 6:648420:d=15062500:TSLSA:kpds5=137:kpds6=200:kpds7=0:
> > > > > > TR=10:P1=0:P2=0:TimeU=1:atmos
> > > > > > col:anl:NAve=0
> > > > > >
> > > > > > All of these 6 records appear to have different GRIB
codes.
> Based
> > on
> > > > the
> > > > > > question you asked, I was expecting to find multiple
records with
> > the
> > > > > > *same* GRIB code but were different in some other way
(i.e.
> > different
> > > > > level
> > > > > > type or different statistical processing type).
> > > > > >
> > > > > > You say these are all AOD, but the plot_data_plane command
you
> sent
> > > > does
> > > > > > *not* work for me:
> > > > > >
> > > > > > plot_data_plane out.grb aod.ps 'name="AOD"; level="A0";'
> > > > > > ERROR : plot_data_plane -> trouble getting field
"name="AOD";
> > > > > level="A0";"
> > > > > > from file "out.grb"
> > > > > >
> > > > > > Please check to see if you have the $MET_GRIB_TABLES
environment
> > > > variable
> > > > > > set. If so, please send me the table(s) to which it
points.
> > > > > >
> > > > > > Thanks,
> > > > > > John
> > > > > >
> > > > > >
> > > > > > On Wed, Jan 17, 2018 at 10:35 AM, Partha Bhattacharjee -
NOAA
> > > Affiliate
> > > > > via
> > > > > > RT <met_help at ucar.edu> wrote:
> > > > > >
> > > > > > >
> > > > > > > Wed Jan 17 10:35:05 2018: Request 83618 was acted upon.
> > > > > > > Transaction: Ticket created by
partha.bhattacharjee at noaa.gov
> > > > > > > Queue: met_help
> > > > > > > Subject: Multiple records of AOD in a grib file
> > > > > > > Owner: Nobody
> > > > > > > Requestors: partha.bhattacharjee at noaa.gov
> > > > > > > Status: new
> > > > > > > Ticket <URL: https://rt.rap.ucar.edu/rt/
> > > > Ticket/Display.html?id=83618
> > > > > >
> > > > > > >
> > > > > > >
> > > > > > > Hi
> > > > > > >
> > > > > > > I have an input GRIB1 file with multiple records of
AOD,before
> I
> > > used
> > > > > > > copygb to interpolate that into (-90:90;-180:180} grid
to use
> for
> > > > > > > Grid-Stat. I am wondering how to read multiple AOD
records in
> > that
> > > > > file,
> > > > > > is
> > > > > > > there any Grib table specific number or template I can
use ?
> > > > > > >
> > > > > > > MET able to read the interpolated file record when I use
this
> > > > command :
> > > > > > > plot_data_plane out.grb aod.ps 'name="AOD";
level="A0";' and
> > > > > generated
> > > > > > > the
> > > > > > > attached figure.
> > > > > > >
> > > > > > > Attached are original and interpolated grib1 files.
> > > > > > >
> > > > > > > Thanks.
> > > > > > > Partha
> > > > > > >
> > > > > > >
> > > > > >
> > > > > >
> > > > >
> > > > >
> > > > > --
> > > > > Partha S. Bhattacharjee, Ph.D.
> > > > > I.M. Systems Group
> > > > > Environmental Modeling Centre
> > > > > NCWCP
> > > > > 5830 University Research Court.
> > > > > College Park, MD 20740
> > > > >
> > > > >
> > > >
> > > >
> > >
> > >
> > > --
> > > Partha S. Bhattacharjee, Ph.D.
> > > I.M. Systems Group
> > > Environmental Modeling Centre
> > > NCWCP
> > > 5830 University Research Court.
> > > College Park, MD 20740
> > >
> > >
> >
> >
>
>
> --
> Partha S. Bhattacharjee, Ph.D.
> I.M. Systems Group
> Environmental Modeling Centre
> NCWCP
> 5830 University Research Court.
> College Park, MD 20740
>
>
------------------------------------------------
Subject: Multiple records of AOD in a grib file
From: Partha Bhattacharjee - NOAA Affiliate
Time: Thu Jan 18 21:40:51 2018
Thanks Julie. I am able to modify accordingly and run Grid-Stat
without
any problem. Do I need to apply same change in WCOSS as well (when it
comes
back) ? And is this problem only with version 6 ?
Partha
On Thu, Jan 18, 2018 at 4:35 PM, Julie Prestopnik via RT
<met_help at ucar.edu>
wrote:
> Hi Partha. Thank you for remembering the WCOSS outage and for
coyping the
> files over to theia! We were able to reproduce the problem you are
> experiencing.
>
> This error is occurring because of the empty string for the units
for AOD
> in the table file grib1_ncep_141_7.txt. Thank you for bringing this
to our
> attention. We have added a ticket to resolve this issue. In the
meantime,
> you can work around this issue by copying over the file
> grib1_ncep_141_7.txt from your met-6.0/share/met/table_files/
directory
> into your current directory. Once you do that, please modify the
line:
>
> 129 141 7 -1 "AOD" "Aerosol Optical Depth" ""
>
> changing the empty string at the end to be "NA" or whatever you'd
like.
> For example,
>
> 129 141 7 -1 "AOD" "Aerosol Optical Depth" "NA"
>
> Then, please set the environment variable MET_GRIB_TABLES to point
to the
> location of this new table. For example:
>
> export MET_GRIB_TABLES= /scratch4/NCEPDEV/global/
> noscrub/Partha.Bhattacharjee/MET_Test/grib1_ncep_141_7.txt
>
> And, then please try running your script again. We were able to get
a
> successful run once we made these changes. Please let us know if
you
> encounter any problems or have any further questions.
>
> Thanks,
> Julie
>
> On Wed, Jan 17, 2018 at 5:07 PM, Partha Bhattacharjee - NOAA
Affiliate via
> RT <met_help at ucar.edu> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=83618 >
> >
> > Thanks John. I realized WCOSS will be out next 2 days, so I copied
the
> > files over Theia
> > at /scratch4/NCEPDEV/global/noscrub/Partha.Bhattacharjee/MET_Test.
> >
> > Partha
> >
> > On Wed, Jan 17, 2018 at 5:43 PM, John Halley Gotway via RT <
> > met_help at ucar.edu> wrote:
> >
> > > Partha,
> > >
> > > Unfortunately, I don't have access to WCOSS, but my colleague,
Julie
> > > Prestopnik, does. I'll ask her to log on and try to replicate
the
> > behavior
> > > you describe. And then we can debug it.
> > >
> > > Thanks,
> > > John
> > >
> > > On Wed, Jan 17, 2018 at 2:34 PM, Partha Bhattacharjee - NOAA
Affiliate
> > via
> > > RT <met_help at ucar.edu> wrote:
> > >
> > > >
> > > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=83618
>
> > > >
> > > > Thanks John. Related to my earlier email, I am getting error
message
> > when
> > > > trying to use MET Series Stat. I am trying to compare
regridded
> > forecast
> > > > file and observational file but getting the following error :
> > > >
> > > > GSL_RNG_TYPE=mt19937
> > > > GSL_RNG_SEED=18446744073587810138
> > > > terminate called after throwing an instance of
'std::logic_error'
> > > > what(): basic_string::_S_construct NULL not valid
> > > >
> > > > I have attached the output.txt along with this email. My files
are in
> > > WCOSS
> > > > at : /global/noscrub/Partha.Bhattacharjee/MET_Test. Original
Grib
> > files
> > > > are
> > > > in the Fcst directory : "AOD550.Grb.t00z.f00" and interpolated
ones
> are
> > > > ones with "grb" at the end. I have used copygb to interpolate
these
> > input
> > > > files which are at 384 X 190 to 360 X 180 to make them same as
Obs
> > files
> > > > are. But not sure if the error is due to not matching grids or
from
> > > > something else.
> > > >
> > > > Thanks.
> > > > Partha
> > > >
> > > > On Wed, Jan 17, 2018 at 4:26 PM, John Halley Gotway via RT <
> > > > met_help at ucar.edu> wrote:
> > > >
> > > > > Partha,
> > > > >
> > > > > Ah, OK. So your file contains 6 GRIB records, each of which
has a
> > > > > different GRIB code (129, 133, 134, 135, 76, and 137).
Those are
> the
> > > > kpds5
> > > > > values listed in the output of wgrib. Running wgrib -V
reveals
> that
> > 5
> > > of
> > > > > the 6 records are encoded as:
> > > > > center 7 subcenter 0 process 82 Table 141 scan: WE:SN
winds(N/S)
> > > > >
> > > > > So this is NCEP GRIB1 table number 141, and that can be
found here:
> > > > > met-6.1/share/met/table_files/grib1_ncep_141_7.txt
> > > > >
> > > > > Here's a selection of lines from that file:
> > > > >
> > > > > 129 141 7 -1 "AOD" "Aerosol Optical Depth" ""
> > > > > 133 141 7 -1 "var133" "undefined" ""
> > > > > 134 141 7 -1 "var134" "undefined" ""
> > > > > 135 141 7 -1 "var135" "undefined" ""
> > > > > 137 141 7 -1 "var137" "undefined" ""
> > > > >
> > > > > And this info is taken from this GRIB table:
> > > > > http://www.nco.ncep.noaa.gov/pmb/docs/on388/table2.html#
> TABLE141
> > > > >
> > > > > So you could just refer to this data using those "var133"
(and so
> on)
> > > > > string:
> > > > >
> > > > > plot_data_plane out.grb plot_rec1.ps 'name="var133";
> > level="L0";'
> > > > >
> > > > > Alternatively, you could specifying the "GRIB1_code" value
instead
> of
> > > the
> > > > > name:
> > > > >
> > > > > plot_data_plane out.grb plot_rec1.ps 'GRIB1_code=133;
> > level="L0";'
> > > > >
> > > > > Both approaches would result in the same plot. Presumably,
since
> you
> > > > have
> > > > > data for codes 133, 134, 135, and 137, perhaps these codes
are no
> > > longer
> > > > > "reserved for future use", as indicated on the NCO website.
> > > > >
> > > > > Hope that helps clarify.
> > > > >
> > > > > Thanks,
> > > > > John
> > > > >
> > > > > On Wed, Jan 17, 2018 at 12:05 PM, Partha Bhattacharjee -
NOAA
> > Affiliate
> > > > via
> > > > > RT <met_help at ucar.edu> wrote:
> > > > >
> > > > > >
> > > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=83618 >
> > > > > >
> > > > > > John,
> > > > > >
> > > > > > That is strange. When I did wgrib to out.grib I am getting
the
> > > > following
> > > > > :
> > > > > >
> > > > > > wgrib out.grb
> > > > > > 1:0:d=15062500:AOD:kpds5=129:kpds6=200:kpds7=0:TR=10:P1=0:
> > > > > > P2=0:TimeU=1:atmos
> > > > > > col:anl:NAve=0
> > > > > > 2:129684:d=15062500:var133:kpds5=133:kpds6=200:kpds7=0:
> > > > > > TR=10:P1=0:P2=0:TimeU=1:atmos
> > > > > > col:anl:NAve=0
> > > > > > 3:259368:d=15062500:var134:kpds5=134:kpds6=200:kpds7=0:
> > > > > > TR=10:P1=0:P2=0:TimeU=1:atmos
> > > > > > col:anl:NAve=0
> > > > > > 4:389052:d=15062500:var135:kpds5=135:kpds6=200:kpds7=0:
> > > > > > TR=10:P1=0:P2=0:TimeU=1:atmos
> > > > > > col:anl:NAve=0
> > > > > > 5:518736:d=15062500:CWAT:kpds5=76:kpds6=200:kpds7=0:TR=
> > > > > > 10:P1=0:P2=0:TimeU=1:atmos
> > > > > > col:anl:NAve=0
> > > > > > 6:648420:d=15062500:var137:kpds5=137:kpds6=200:kpds7=0:
> > > > > > TR=10:P1=0:P2=0:TimeU=1:atmos
> > > > > > col:anl:NAve=0
> > > > > >
> > > > > > and I am running in wcoss and wgrib points to
> > > > > > which wgrib
> > > > > > alias wgrib='/nwprod/util/exec/wgrib'
> > > > > > /nwprod/util/exec/wgrib
> > > > > >
> > > > > > So, I used "AOD" in the first record to plot like :
> > > > > >
> > > > > > plot_data_plane out.grb aod.ps 'name="AOD"; level="L0";'
-v 6
> > > > > > DEBUG 1: Opening data file: out.grb
> > > > > > DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() ->
created
> > new
> > > > > > Met2dDataFile object of type "FileType_Gb1".
> > > > > > DEBUG 4: VarInfoFactory::new_var_info() -> created new
VarInfo
> > object
> > > > of
> > > > > > type "FileType_Gb1".
> > > > > > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found
range
> match
> > > for
> > > > > > VarInfo "AOD/L0" in GRIB record 1 of GRIB file "out.grb".
> > > > > > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 1
GRIB
> > > records
> > > > > > matching VarInfo "AOD/L0" in GRIB file "out.grb".
> > > > > > DEBUG 1: Creating postscript file: aod.ps
> > > > > > DEBUG 4: draw_map() -> Plotting map data file:
> > > > > >
/awc/noscrub/Julie.Prestopnik/met/6.0/share/met/map/country_data
> > > > > > DEBUG 4: draw_map() -> Plotting map data file:
> > > > > > /awc/noscrub/Julie.Prestopnik/met/6.0/share/met/map/usa_
> state_data
> > > > > > DEBUG 4: draw_map() -> Plotting map data file:
> > > > > > /awc/noscrub/Julie.Prestopnik/met/6.0/share/met/map/major_
> > lakes_data
> > > > > >
> > > > > > I do not think I have $MET_GRIB_TABLES environment
variable set.
> > > > > >
> > > > > > Partha
> > > > > >
> > > > > > On Wed, Jan 17, 2018 at 1:33 PM, John Halley Gotway via RT
<
> > > > > > met_help at ucar.edu> wrote:
> > > > > >
> > > > > > > Partha,
> > > > > > >
> > > > > > > I'm confused, when I run wgrib on the data you sent, I
don't
> see
> > > any
> > > > > > > mention of AOD. But I do get a warning message about an
> > undefined
> > > > > > > parameter table:
> > > > > > >
> > > > > > > [johnhg at number5]% wgrib out.grb
> > > > > > >
> > > > > > > Undefined parameter table (center 7-0 table 141), using
> NCEP-opn
> > > > > > > 1:0:d=15062500:MSLMA:kpds5=129:kpds6=200:kpds7=0:TR=10:
> > > > > > > P1=0:P2=0:TimeU=1:atmos
> > > > > > > col:anl:NAve=0
> > > > > > >
2:129684:d=15062500:KX:kpds5=133:kpds6=200:kpds7=0:TR=10:
> > > > > > > P1=0:P2=0:TimeU=1:atmos
> > > > > > > col:anl:NAve=0
> > > > > > >
3:259368:d=15062500:SX:kpds5=134:kpds6=200:kpds7=0:TR=10:
> > > > > > > P1=0:P2=0:TimeU=1:atmos
> > > > > > > col:anl:NAve=0
> > > > > > > 4:389052:d=15062500:MCONV:kpds5=135:kpds6=200:kpds7=0:
> > > > > > > TR=10:P1=0:P2=0:TimeU=1:atmos
> > > > > > > col:anl:NAve=0
> > > > > > > 5:518736:d=15062500:CWAT:kpds5=76:kpds6=200:kpds7=0:TR=
> > > > > > > 10:P1=0:P2=0:TimeU=1:atmos
> > > > > > > col:anl:NAve=0
> > > > > > > 6:648420:d=15062500:TSLSA:kpds5=137:kpds6=200:kpds7=0:
> > > > > > > TR=10:P1=0:P2=0:TimeU=1:atmos
> > > > > > > col:anl:NAve=0
> > > > > > >
> > > > > > > All of these 6 records appear to have different GRIB
codes.
> > Based
> > > on
> > > > > the
> > > > > > > question you asked, I was expecting to find multiple
records
> with
> > > the
> > > > > > > *same* GRIB code but were different in some other way
(i.e.
> > > different
> > > > > > level
> > > > > > > type or different statistical processing type).
> > > > > > >
> > > > > > > You say these are all AOD, but the plot_data_plane
command you
> > sent
> > > > > does
> > > > > > > *not* work for me:
> > > > > > >
> > > > > > > plot_data_plane out.grb aod.ps 'name="AOD";
level="A0";'
> > > > > > > ERROR : plot_data_plane -> trouble getting field
"name="AOD";
> > > > > > level="A0";"
> > > > > > > from file "out.grb"
> > > > > > >
> > > > > > > Please check to see if you have the $MET_GRIB_TABLES
> environment
> > > > > variable
> > > > > > > set. If so, please send me the table(s) to which it
points.
> > > > > > >
> > > > > > > Thanks,
> > > > > > > John
> > > > > > >
> > > > > > >
> > > > > > > On Wed, Jan 17, 2018 at 10:35 AM, Partha Bhattacharjee -
NOAA
> > > > Affiliate
> > > > > > via
> > > > > > > RT <met_help at ucar.edu> wrote:
> > > > > > >
> > > > > > > >
> > > > > > > > Wed Jan 17 10:35:05 2018: Request 83618 was acted
upon.
> > > > > > > > Transaction: Ticket created by
partha.bhattacharjee at noaa.gov
> > > > > > > > Queue: met_help
> > > > > > > > Subject: Multiple records of AOD in a grib file
> > > > > > > > Owner: Nobody
> > > > > > > > Requestors: partha.bhattacharjee at noaa.gov
> > > > > > > > Status: new
> > > > > > > > Ticket <URL: https://rt.rap.ucar.edu/rt/
> > > > > Ticket/Display.html?id=83618
> > > > > > >
> > > > > > > >
> > > > > > > >
> > > > > > > > Hi
> > > > > > > >
> > > > > > > > I have an input GRIB1 file with multiple records of
> AOD,before
> > I
> > > > used
> > > > > > > > copygb to interpolate that into (-90:90;-180:180} grid
to use
> > for
> > > > > > > > Grid-Stat. I am wondering how to read multiple AOD
records in
> > > that
> > > > > > file,
> > > > > > > is
> > > > > > > > there any Grib table specific number or template I can
use ?
> > > > > > > >
> > > > > > > > MET able to read the interpolated file record when I
use this
> > > > > command :
> > > > > > > > plot_data_plane out.grb aod.ps 'name="AOD";
level="A0";'
> and
> > > > > > generated
> > > > > > > > the
> > > > > > > > attached figure.
> > > > > > > >
> > > > > > > > Attached are original and interpolated grib1 files.
> > > > > > > >
> > > > > > > > Thanks.
> > > > > > > > Partha
> > > > > > > >
> > > > > > > >
> > > > > > >
> > > > > > >
> > > > > >
> > > > > >
> > > > > > --
> > > > > > Partha S. Bhattacharjee, Ph.D.
> > > > > > I.M. Systems Group
> > > > > > Environmental Modeling Centre
> > > > > > NCWCP
> > > > > > 5830 University Research Court.
> > > > > > College Park, MD 20740
> > > > > >
> > > > > >
> > > > >
> > > > >
> > > >
> > > >
> > > > --
> > > > Partha S. Bhattacharjee, Ph.D.
> > > > I.M. Systems Group
> > > > Environmental Modeling Centre
> > > > NCWCP
> > > > 5830 University Research Court.
> > > > College Park, MD 20740
> > > >
> > > >
> > >
> > >
> >
> >
> > --
> > Partha S. Bhattacharjee, Ph.D.
> > I.M. Systems Group
> > Environmental Modeling Centre
> > NCWCP
> > 5830 University Research Court.
> > College Park, MD 20740
> >
> >
>
>
--
Partha S. Bhattacharjee, Ph.D.
I.M. Systems Group
Environmental Modeling Centre
NCWCP
5830 University Research Court.
College Park, MD 20740
------------------------------------------------
Subject: Multiple records of AOD in a grib file
From: Julie Prestopnik
Time: Fri Jan 19 08:34:22 2018
Hi Partha.
Yes, this same change will apply on WCOSS. This problem exists in 6.1
as
well. We will put in a fix in a patched version.
Julie
On Thu, Jan 18, 2018 at 9:40 PM, Partha Bhattacharjee - NOAA Affiliate
via
RT <met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=83618 >
>
> Thanks Julie. I am able to modify accordingly and run Grid-Stat
without
> any problem. Do I need to apply same change in WCOSS as well (when
it comes
> back) ? And is this problem only with version 6 ?
>
> Partha
>
> On Thu, Jan 18, 2018 at 4:35 PM, Julie Prestopnik via RT <
> met_help at ucar.edu>
> wrote:
>
> > Hi Partha. Thank you for remembering the WCOSS outage and for
coyping the
> > files over to theia! We were able to reproduce the problem you
are
> > experiencing.
> >
> > This error is occurring because of the empty string for the units
for AOD
> > in the table file grib1_ncep_141_7.txt. Thank you for bringing
this to
> our
> > attention. We have added a ticket to resolve this issue. In the
> meantime,
> > you can work around this issue by copying over the file
> > grib1_ncep_141_7.txt from your met-6.0/share/met/table_files/
directory
> > into your current directory. Once you do that, please modify the
line:
> >
> > 129 141 7 -1 "AOD" "Aerosol Optical Depth" ""
> >
> > changing the empty string at the end to be "NA" or whatever you'd
like.
> > For example,
> >
> > 129 141 7 -1 "AOD" "Aerosol Optical Depth" "NA"
> >
> > Then, please set the environment variable MET_GRIB_TABLES to point
to the
> > location of this new table. For example:
> >
> > export MET_GRIB_TABLES= /scratch4/NCEPDEV/global/
> > noscrub/Partha.Bhattacharjee/MET_Test/grib1_ncep_141_7.txt
> >
> > And, then please try running your script again. We were able to
get a
> > successful run once we made these changes. Please let us know if
you
> > encounter any problems or have any further questions.
> >
> > Thanks,
> > Julie
> >
> > On Wed, Jan 17, 2018 at 5:07 PM, Partha Bhattacharjee - NOAA
Affiliate
> via
> > RT <met_help at ucar.edu> wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=83618 >
> > >
> > > Thanks John. I realized WCOSS will be out next 2 days, so I
copied the
> > > files over Theia
> > > at
/scratch4/NCEPDEV/global/noscrub/Partha.Bhattacharjee/MET_Test.
> > >
> > > Partha
> > >
> > > On Wed, Jan 17, 2018 at 5:43 PM, John Halley Gotway via RT <
> > > met_help at ucar.edu> wrote:
> > >
> > > > Partha,
> > > >
> > > > Unfortunately, I don't have access to WCOSS, but my colleague,
Julie
> > > > Prestopnik, does. I'll ask her to log on and try to replicate
the
> > > behavior
> > > > you describe. And then we can debug it.
> > > >
> > > > Thanks,
> > > > John
> > > >
> > > > On Wed, Jan 17, 2018 at 2:34 PM, Partha Bhattacharjee - NOAA
> Affiliate
> > > via
> > > > RT <met_help at ucar.edu> wrote:
> > > >
> > > > >
> > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=83618 >
> > > > >
> > > > > Thanks John. Related to my earlier email, I am getting error
> message
> > > when
> > > > > trying to use MET Series Stat. I am trying to compare
regridded
> > > forecast
> > > > > file and observational file but getting the following error
:
> > > > >
> > > > > GSL_RNG_TYPE=mt19937
> > > > > GSL_RNG_SEED=18446744073587810138
> > > > > terminate called after throwing an instance of
'std::logic_error'
> > > > > what(): basic_string::_S_construct NULL not valid
> > > > >
> > > > > I have attached the output.txt along with this email. My
files are
> in
> > > > WCOSS
> > > > > at : /global/noscrub/Partha.Bhattacharjee/MET_Test. Original
Grib
> > > files
> > > > > are
> > > > > in the Fcst directory : "AOD550.Grb.t00z.f00" and
interpolated ones
> > are
> > > > > ones with "grb" at the end. I have used copygb to
interpolate these
> > > input
> > > > > files which are at 384 X 190 to 360 X 180 to make them same
as Obs
> > > files
> > > > > are. But not sure if the error is due to not matching grids
or from
> > > > > something else.
> > > > >
> > > > > Thanks.
> > > > > Partha
> > > > >
> > > > > On Wed, Jan 17, 2018 at 4:26 PM, John Halley Gotway via RT <
> > > > > met_help at ucar.edu> wrote:
> > > > >
> > > > > > Partha,
> > > > > >
> > > > > > Ah, OK. So your file contains 6 GRIB records, each of
which has
> a
> > > > > > different GRIB code (129, 133, 134, 135, 76, and 137).
Those are
> > the
> > > > > kpds5
> > > > > > values listed in the output of wgrib. Running wgrib -V
reveals
> > that
> > > 5
> > > > of
> > > > > > the 6 records are encoded as:
> > > > > > center 7 subcenter 0 process 82 Table 141 scan: WE:SN
> winds(N/S)
> > > > > >
> > > > > > So this is NCEP GRIB1 table number 141, and that can be
found
> here:
> > > > > > met-6.1/share/met/table_files/grib1_ncep_141_7.txt
> > > > > >
> > > > > > Here's a selection of lines from that file:
> > > > > >
> > > > > > 129 141 7 -1 "AOD" "Aerosol Optical Depth" ""
> > > > > > 133 141 7 -1 "var133" "undefined" ""
> > > > > > 134 141 7 -1 "var134" "undefined" ""
> > > > > > 135 141 7 -1 "var135" "undefined" ""
> > > > > > 137 141 7 -1 "var137" "undefined" ""
> > > > > >
> > > > > > And this info is taken from this GRIB table:
> > > > > >
http://www.nco.ncep.noaa.gov/pmb/docs/on388/table2.html#
> > TABLE141
> > > > > >
> > > > > > So you could just refer to this data using those "var133"
(and so
> > on)
> > > > > > string:
> > > > > >
> > > > > > plot_data_plane out.grb plot_rec1.ps 'name="var133";
> > > level="L0";'
> > > > > >
> > > > > > Alternatively, you could specifying the "GRIB1_code" value
> instead
> > of
> > > > the
> > > > > > name:
> > > > > >
> > > > > > plot_data_plane out.grb plot_rec1.ps 'GRIB1_code=133;
> > > level="L0";'
> > > > > >
> > > > > > Both approaches would result in the same plot.
Presumably, since
> > you
> > > > > have
> > > > > > data for codes 133, 134, 135, and 137, perhaps these codes
are no
> > > > longer
> > > > > > "reserved for future use", as indicated on the NCO
website.
> > > > > >
> > > > > > Hope that helps clarify.
> > > > > >
> > > > > > Thanks,
> > > > > > John
> > > > > >
> > > > > > On Wed, Jan 17, 2018 at 12:05 PM, Partha Bhattacharjee -
NOAA
> > > Affiliate
> > > > > via
> > > > > > RT <met_help at ucar.edu> wrote:
> > > > > >
> > > > > > >
> > > > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=83618
> >
> > > > > > >
> > > > > > > John,
> > > > > > >
> > > > > > > That is strange. When I did wgrib to out.grib I am
getting the
> > > > > following
> > > > > > :
> > > > > > >
> > > > > > > wgrib out.grb
> > > > > > >
1:0:d=15062500:AOD:kpds5=129:kpds6=200:kpds7=0:TR=10:P1=0:
> > > > > > > P2=0:TimeU=1:atmos
> > > > > > > col:anl:NAve=0
> > > > > > > 2:129684:d=15062500:var133:kpds5=133:kpds6=200:kpds7=0:
> > > > > > > TR=10:P1=0:P2=0:TimeU=1:atmos
> > > > > > > col:anl:NAve=0
> > > > > > > 3:259368:d=15062500:var134:kpds5=134:kpds6=200:kpds7=0:
> > > > > > > TR=10:P1=0:P2=0:TimeU=1:atmos
> > > > > > > col:anl:NAve=0
> > > > > > > 4:389052:d=15062500:var135:kpds5=135:kpds6=200:kpds7=0:
> > > > > > > TR=10:P1=0:P2=0:TimeU=1:atmos
> > > > > > > col:anl:NAve=0
> > > > > > > 5:518736:d=15062500:CWAT:kpds5=76:kpds6=200:kpds7=0:TR=
> > > > > > > 10:P1=0:P2=0:TimeU=1:atmos
> > > > > > > col:anl:NAve=0
> > > > > > > 6:648420:d=15062500:var137:kpds5=137:kpds6=200:kpds7=0:
> > > > > > > TR=10:P1=0:P2=0:TimeU=1:atmos
> > > > > > > col:anl:NAve=0
> > > > > > >
> > > > > > > and I am running in wcoss and wgrib points to
> > > > > > > which wgrib
> > > > > > > alias wgrib='/nwprod/util/exec/wgrib'
> > > > > > > /nwprod/util/exec/wgrib
> > > > > > >
> > > > > > > So, I used "AOD" in the first record to plot like :
> > > > > > >
> > > > > > > plot_data_plane out.grb aod.ps 'name="AOD";
level="L0";' -v 6
> > > > > > > DEBUG 1: Opening data file: out.grb
> > > > > > > DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() ->
> created
> > > new
> > > > > > > Met2dDataFile object of type "FileType_Gb1".
> > > > > > > DEBUG 4: VarInfoFactory::new_var_info() -> created new
VarInfo
> > > object
> > > > > of
> > > > > > > type "FileType_Gb1".
> > > > > > > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found
range
> > match
> > > > for
> > > > > > > VarInfo "AOD/L0" in GRIB record 1 of GRIB file
"out.grb".
> > > > > > > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 1
GRIB
> > > > records
> > > > > > > matching VarInfo "AOD/L0" in GRIB file "out.grb".
> > > > > > > DEBUG 1: Creating postscript file: aod.ps
> > > > > > > DEBUG 4: draw_map() -> Plotting map data file:
> > > > > > >
/awc/noscrub/Julie.Prestopnik/met/6.0/share/met/map/country_
> data
> > > > > > > DEBUG 4: draw_map() -> Plotting map data file:
> > > > > > > /awc/noscrub/Julie.Prestopnik/met/6.0/share/met/map/usa_
> > state_data
> > > > > > > DEBUG 4: draw_map() -> Plotting map data file:
> > > > > > >
/awc/noscrub/Julie.Prestopnik/met/6.0/share/met/map/major_
> > > lakes_data
> > > > > > >
> > > > > > > I do not think I have $MET_GRIB_TABLES environment
variable
> set.
> > > > > > >
> > > > > > > Partha
> > > > > > >
> > > > > > > On Wed, Jan 17, 2018 at 1:33 PM, John Halley Gotway via
RT <
> > > > > > > met_help at ucar.edu> wrote:
> > > > > > >
> > > > > > > > Partha,
> > > > > > > >
> > > > > > > > I'm confused, when I run wgrib on the data you sent, I
don't
> > see
> > > > any
> > > > > > > > mention of AOD. But I do get a warning message about
an
> > > undefined
> > > > > > > > parameter table:
> > > > > > > >
> > > > > > > > [johnhg at number5]% wgrib out.grb
> > > > > > > >
> > > > > > > > Undefined parameter table (center 7-0 table 141),
using
> > NCEP-opn
> > > > > > > >
1:0:d=15062500:MSLMA:kpds5=129:kpds6=200:kpds7=0:TR=10:
> > > > > > > > P1=0:P2=0:TimeU=1:atmos
> > > > > > > > col:anl:NAve=0
> > > > > > > >
2:129684:d=15062500:KX:kpds5=133:kpds6=200:kpds7=0:TR=10:
> > > > > > > > P1=0:P2=0:TimeU=1:atmos
> > > > > > > > col:anl:NAve=0
> > > > > > > >
3:259368:d=15062500:SX:kpds5=134:kpds6=200:kpds7=0:TR=10:
> > > > > > > > P1=0:P2=0:TimeU=1:atmos
> > > > > > > > col:anl:NAve=0
> > > > > > > > 4:389052:d=15062500:MCONV:kpds5=135:kpds6=200:kpds7=0:
> > > > > > > > TR=10:P1=0:P2=0:TimeU=1:atmos
> > > > > > > > col:anl:NAve=0
> > > > > > > >
5:518736:d=15062500:CWAT:kpds5=76:kpds6=200:kpds7=0:TR=
> > > > > > > > 10:P1=0:P2=0:TimeU=1:atmos
> > > > > > > > col:anl:NAve=0
> > > > > > > > 6:648420:d=15062500:TSLSA:kpds5=137:kpds6=200:kpds7=0:
> > > > > > > > TR=10:P1=0:P2=0:TimeU=1:atmos
> > > > > > > > col:anl:NAve=0
> > > > > > > >
> > > > > > > > All of these 6 records appear to have different GRIB
codes.
> > > Based
> > > > on
> > > > > > the
> > > > > > > > question you asked, I was expecting to find multiple
records
> > with
> > > > the
> > > > > > > > *same* GRIB code but were different in some other way
(i.e.
> > > > different
> > > > > > > level
> > > > > > > > type or different statistical processing type).
> > > > > > > >
> > > > > > > > You say these are all AOD, but the plot_data_plane
command
> you
> > > sent
> > > > > > does
> > > > > > > > *not* work for me:
> > > > > > > >
> > > > > > > > plot_data_plane out.grb aod.ps 'name="AOD";
level="A0";'
> > > > > > > > ERROR : plot_data_plane -> trouble getting field
> "name="AOD";
> > > > > > > level="A0";"
> > > > > > > > from file "out.grb"
> > > > > > > >
> > > > > > > > Please check to see if you have the $MET_GRIB_TABLES
> > environment
> > > > > > variable
> > > > > > > > set. If so, please send me the table(s) to which it
points.
> > > > > > > >
> > > > > > > > Thanks,
> > > > > > > > John
> > > > > > > >
> > > > > > > >
> > > > > > > > On Wed, Jan 17, 2018 at 10:35 AM, Partha Bhattacharjee
- NOAA
> > > > > Affiliate
> > > > > > > via
> > > > > > > > RT <met_help at ucar.edu> wrote:
> > > > > > > >
> > > > > > > > >
> > > > > > > > > Wed Jan 17 10:35:05 2018: Request 83618 was acted
upon.
> > > > > > > > > Transaction: Ticket created by
> partha.bhattacharjee at noaa.gov
> > > > > > > > > Queue: met_help
> > > > > > > > > Subject: Multiple records of AOD in a grib file
> > > > > > > > > Owner: Nobody
> > > > > > > > > Requestors: partha.bhattacharjee at noaa.gov
> > > > > > > > > Status: new
> > > > > > > > > Ticket <URL: https://rt.rap.ucar.edu/rt/
> > > > > > Ticket/Display.html?id=83618
> > > > > > > >
> > > > > > > > >
> > > > > > > > >
> > > > > > > > > Hi
> > > > > > > > >
> > > > > > > > > I have an input GRIB1 file with multiple records of
> > AOD,before
> > > I
> > > > > used
> > > > > > > > > copygb to interpolate that into (-90:90;-180:180}
grid to
> use
> > > for
> > > > > > > > > Grid-Stat. I am wondering how to read multiple AOD
records
> in
> > > > that
> > > > > > > file,
> > > > > > > > is
> > > > > > > > > there any Grib table specific number or template I
can use
> ?
> > > > > > > > >
> > > > > > > > > MET able to read the interpolated file record when I
use
> this
> > > > > > command :
> > > > > > > > > plot_data_plane out.grb aod.ps 'name="AOD";
level="A0";'
> > and
> > > > > > > generated
> > > > > > > > > the
> > > > > > > > > attached figure.
> > > > > > > > >
> > > > > > > > > Attached are original and interpolated grib1 files.
> > > > > > > > >
> > > > > > > > > Thanks.
> > > > > > > > > Partha
> > > > > > > > >
> > > > > > > > >
> > > > > > > >
> > > > > > > >
> > > > > > >
> > > > > > >
> > > > > > > --
> > > > > > > Partha S. Bhattacharjee, Ph.D.
> > > > > > > I.M. Systems Group
> > > > > > > Environmental Modeling Centre
> > > > > > > NCWCP
> > > > > > > 5830 University Research Court.
> > > > > > > College Park, MD 20740
> > > > > > >
> > > > > > >
> > > > > >
> > > > > >
> > > > >
> > > > >
> > > > > --
> > > > > Partha S. Bhattacharjee, Ph.D.
> > > > > I.M. Systems Group
> > > > > Environmental Modeling Centre
> > > > > NCWCP
> > > > > 5830 University Research Court.
> > > > > College Park, MD 20740
> > > > >
> > > > >
> > > >
> > > >
> > >
> > >
> > > --
> > > Partha S. Bhattacharjee, Ph.D.
> > > I.M. Systems Group
> > > Environmental Modeling Centre
> > > NCWCP
> > > 5830 University Research Court.
> > > College Park, MD 20740
> > >
> > >
> >
> >
>
>
> --
> Partha S. Bhattacharjee, Ph.D.
> I.M. Systems Group
> Environmental Modeling Centre
> NCWCP
> 5830 University Research Court.
> College Park, MD 20740
>
>
------------------------------------------------
Subject: Multiple records of AOD in a grib file
From: Partha Bhattacharjee - NOAA Affiliate
Time: Mon Jan 22 07:09:34 2018
Thanks Julie. I did it successfully in WCOSS as well. You may close
this
ticket.
Partha
On Fri, Jan 19, 2018 at 10:34 AM, Julie Prestopnik via RT
<met_help at ucar.edu
> wrote:
> Hi Partha.
>
> Yes, this same change will apply on WCOSS. This problem exists in
6.1 as
> well. We will put in a fix in a patched version.
>
> Julie
>
> On Thu, Jan 18, 2018 at 9:40 PM, Partha Bhattacharjee - NOAA
Affiliate via
> RT <met_help at ucar.edu> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=83618 >
> >
> > Thanks Julie. I am able to modify accordingly and run Grid-Stat
without
> > any problem. Do I need to apply same change in WCOSS as well (when
it
> comes
> > back) ? And is this problem only with version 6 ?
> >
> > Partha
> >
> > On Thu, Jan 18, 2018 at 4:35 PM, Julie Prestopnik via RT <
> > met_help at ucar.edu>
> > wrote:
> >
> > > Hi Partha. Thank you for remembering the WCOSS outage and for
coyping
> the
> > > files over to theia! We were able to reproduce the problem you
are
> > > experiencing.
> > >
> > > This error is occurring because of the empty string for the
units for
> AOD
> > > in the table file grib1_ncep_141_7.txt. Thank you for bringing
this to
> > our
> > > attention. We have added a ticket to resolve this issue. In
the
> > meantime,
> > > you can work around this issue by copying over the file
> > > grib1_ncep_141_7.txt from your met-6.0/share/met/table_files/
directory
> > > into your current directory. Once you do that, please modify
the line:
> > >
> > > 129 141 7 -1 "AOD" "Aerosol Optical Depth" ""
> > >
> > > changing the empty string at the end to be "NA" or whatever
you'd like.
> > > For example,
> > >
> > > 129 141 7 -1 "AOD" "Aerosol Optical Depth" "NA"
> > >
> > > Then, please set the environment variable MET_GRIB_TABLES to
point to
> the
> > > location of this new table. For example:
> > >
> > > export MET_GRIB_TABLES= /scratch4/NCEPDEV/global/
> > > noscrub/Partha.Bhattacharjee/MET_Test/grib1_ncep_141_7.txt
> > >
> > > And, then please try running your script again. We were able to
get a
> > > successful run once we made these changes. Please let us know
if you
> > > encounter any problems or have any further questions.
> > >
> > > Thanks,
> > > Julie
> > >
> > > On Wed, Jan 17, 2018 at 5:07 PM, Partha Bhattacharjee - NOAA
Affiliate
> > via
> > > RT <met_help at ucar.edu> wrote:
> > >
> > > >
> > > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=83618
>
> > > >
> > > > Thanks John. I realized WCOSS will be out next 2 days, so I
copied
> the
> > > > files over Theia
> > > > at
/scratch4/NCEPDEV/global/noscrub/Partha.Bhattacharjee/MET_Test.
> > > >
> > > > Partha
> > > >
> > > > On Wed, Jan 17, 2018 at 5:43 PM, John Halley Gotway via RT <
> > > > met_help at ucar.edu> wrote:
> > > >
> > > > > Partha,
> > > > >
> > > > > Unfortunately, I don't have access to WCOSS, but my
colleague,
> Julie
> > > > > Prestopnik, does. I'll ask her to log on and try to
replicate the
> > > > behavior
> > > > > you describe. And then we can debug it.
> > > > >
> > > > > Thanks,
> > > > > John
> > > > >
> > > > > On Wed, Jan 17, 2018 at 2:34 PM, Partha Bhattacharjee - NOAA
> > Affiliate
> > > > via
> > > > > RT <met_help at ucar.edu> wrote:
> > > > >
> > > > > >
> > > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=83618 >
> > > > > >
> > > > > > Thanks John. Related to my earlier email, I am getting
error
> > message
> > > > when
> > > > > > trying to use MET Series Stat. I am trying to compare
regridded
> > > > forecast
> > > > > > file and observational file but getting the following
error :
> > > > > >
> > > > > > GSL_RNG_TYPE=mt19937
> > > > > > GSL_RNG_SEED=18446744073587810138
> > > > > > terminate called after throwing an instance of
'std::logic_error'
> > > > > > what(): basic_string::_S_construct NULL not valid
> > > > > >
> > > > > > I have attached the output.txt along with this email. My
files
> are
> > in
> > > > > WCOSS
> > > > > > at : /global/noscrub/Partha.Bhattacharjee/MET_Test.
Original
> Grib
> > > > files
> > > > > > are
> > > > > > in the Fcst directory : "AOD550.Grb.t00z.f00" and
interpolated
> ones
> > > are
> > > > > > ones with "grb" at the end. I have used copygb to
interpolate
> these
> > > > input
> > > > > > files which are at 384 X 190 to 360 X 180 to make them
same as
> Obs
> > > > files
> > > > > > are. But not sure if the error is due to not matching
grids or
> from
> > > > > > something else.
> > > > > >
> > > > > > Thanks.
> > > > > > Partha
> > > > > >
> > > > > > On Wed, Jan 17, 2018 at 4:26 PM, John Halley Gotway via RT
<
> > > > > > met_help at ucar.edu> wrote:
> > > > > >
> > > > > > > Partha,
> > > > > > >
> > > > > > > Ah, OK. So your file contains 6 GRIB records, each of
which
> has
> > a
> > > > > > > different GRIB code (129, 133, 134, 135, 76, and 137).
Those
> are
> > > the
> > > > > > kpds5
> > > > > > > values listed in the output of wgrib. Running wgrib -V
reveals
> > > that
> > > > 5
> > > > > of
> > > > > > > the 6 records are encoded as:
> > > > > > > center 7 subcenter 0 process 82 Table 141 scan: WE:SN
> > winds(N/S)
> > > > > > >
> > > > > > > So this is NCEP GRIB1 table number 141, and that can be
found
> > here:
> > > > > > > met-6.1/share/met/table_files/grib1_ncep_141_7.txt
> > > > > > >
> > > > > > > Here's a selection of lines from that file:
> > > > > > >
> > > > > > > 129 141 7 -1 "AOD" "Aerosol Optical Depth" ""
> > > > > > > 133 141 7 -1 "var133" "undefined" ""
> > > > > > > 134 141 7 -1 "var134" "undefined" ""
> > > > > > > 135 141 7 -1 "var135" "undefined" ""
> > > > > > > 137 141 7 -1 "var137" "undefined" ""
> > > > > > >
> > > > > > > And this info is taken from this GRIB table:
> > > > > > >
http://www.nco.ncep.noaa.gov/pmb/docs/on388/table2.html#
> > > TABLE141
> > > > > > >
> > > > > > > So you could just refer to this data using those
"var133" (and
> so
> > > on)
> > > > > > > string:
> > > > > > >
> > > > > > > plot_data_plane out.grb plot_rec1.ps 'name="var133";
> > > > level="L0";'
> > > > > > >
> > > > > > > Alternatively, you could specifying the "GRIB1_code"
value
> > instead
> > > of
> > > > > the
> > > > > > > name:
> > > > > > >
> > > > > > > plot_data_plane out.grb plot_rec1.ps
'GRIB1_code=133;
> > > > level="L0";'
> > > > > > >
> > > > > > > Both approaches would result in the same plot.
Presumably,
> since
> > > you
> > > > > > have
> > > > > > > data for codes 133, 134, 135, and 137, perhaps these
codes are
> no
> > > > > longer
> > > > > > > "reserved for future use", as indicated on the NCO
website.
> > > > > > >
> > > > > > > Hope that helps clarify.
> > > > > > >
> > > > > > > Thanks,
> > > > > > > John
> > > > > > >
> > > > > > > On Wed, Jan 17, 2018 at 12:05 PM, Partha Bhattacharjee -
NOAA
> > > > Affiliate
> > > > > > via
> > > > > > > RT <met_help at ucar.edu> wrote:
> > > > > > >
> > > > > > > >
> > > > > > > > <URL: https://rt.rap.ucar.edu/rt/
> Ticket/Display.html?id=83618
> > >
> > > > > > > >
> > > > > > > > John,
> > > > > > > >
> > > > > > > > That is strange. When I did wgrib to out.grib I am
getting
> the
> > > > > > following
> > > > > > > :
> > > > > > > >
> > > > > > > > wgrib out.grb
> > > > > > > >
1:0:d=15062500:AOD:kpds5=129:kpds6=200:kpds7=0:TR=10:P1=0:
> > > > > > > > P2=0:TimeU=1:atmos
> > > > > > > > col:anl:NAve=0
> > > > > > > >
2:129684:d=15062500:var133:kpds5=133:kpds6=200:kpds7=0:
> > > > > > > > TR=10:P1=0:P2=0:TimeU=1:atmos
> > > > > > > > col:anl:NAve=0
> > > > > > > >
3:259368:d=15062500:var134:kpds5=134:kpds6=200:kpds7=0:
> > > > > > > > TR=10:P1=0:P2=0:TimeU=1:atmos
> > > > > > > > col:anl:NAve=0
> > > > > > > >
4:389052:d=15062500:var135:kpds5=135:kpds6=200:kpds7=0:
> > > > > > > > TR=10:P1=0:P2=0:TimeU=1:atmos
> > > > > > > > col:anl:NAve=0
> > > > > > > >
5:518736:d=15062500:CWAT:kpds5=76:kpds6=200:kpds7=0:TR=
> > > > > > > > 10:P1=0:P2=0:TimeU=1:atmos
> > > > > > > > col:anl:NAve=0
> > > > > > > >
6:648420:d=15062500:var137:kpds5=137:kpds6=200:kpds7=0:
> > > > > > > > TR=10:P1=0:P2=0:TimeU=1:atmos
> > > > > > > > col:anl:NAve=0
> > > > > > > >
> > > > > > > > and I am running in wcoss and wgrib points to
> > > > > > > > which wgrib
> > > > > > > > alias wgrib='/nwprod/util/exec/wgrib'
> > > > > > > > /nwprod/util/exec/wgrib
> > > > > > > >
> > > > > > > > So, I used "AOD" in the first record to plot like :
> > > > > > > >
> > > > > > > > plot_data_plane out.grb aod.ps 'name="AOD";
level="L0";'
> -v 6
> > > > > > > > DEBUG 1: Opening data file: out.grb
> > > > > > > > DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file()
->
> > created
> > > > new
> > > > > > > > Met2dDataFile object of type "FileType_Gb1".
> > > > > > > > DEBUG 4: VarInfoFactory::new_var_info() -> created new
> VarInfo
> > > > object
> > > > > > of
> > > > > > > > type "FileType_Gb1".
> > > > > > > > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found
range
> > > match
> > > > > for
> > > > > > > > VarInfo "AOD/L0" in GRIB record 1 of GRIB file
"out.grb".
> > > > > > > > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found
1
> GRIB
> > > > > records
> > > > > > > > matching VarInfo "AOD/L0" in GRIB file "out.grb".
> > > > > > > > DEBUG 1: Creating postscript file: aod.ps
> > > > > > > > DEBUG 4: draw_map() -> Plotting map data file:
> > > > > > > >
/awc/noscrub/Julie.Prestopnik/met/6.0/share/met/map/country_
> > data
> > > > > > > > DEBUG 4: draw_map() -> Plotting map data file:
> > > > > > > >
/awc/noscrub/Julie.Prestopnik/met/6.0/share/met/map/usa_
> > > state_data
> > > > > > > > DEBUG 4: draw_map() -> Plotting map data file:
> > > > > > > >
/awc/noscrub/Julie.Prestopnik/met/6.0/share/met/map/major_
> > > > lakes_data
> > > > > > > >
> > > > > > > > I do not think I have $MET_GRIB_TABLES environment
variable
> > set.
> > > > > > > >
> > > > > > > > Partha
> > > > > > > >
> > > > > > > > On Wed, Jan 17, 2018 at 1:33 PM, John Halley Gotway
via RT <
> > > > > > > > met_help at ucar.edu> wrote:
> > > > > > > >
> > > > > > > > > Partha,
> > > > > > > > >
> > > > > > > > > I'm confused, when I run wgrib on the data you sent,
I
> don't
> > > see
> > > > > any
> > > > > > > > > mention of AOD. But I do get a warning message
about an
> > > > undefined
> > > > > > > > > parameter table:
> > > > > > > > >
> > > > > > > > > [johnhg at number5]% wgrib out.grb
> > > > > > > > >
> > > > > > > > > Undefined parameter table (center 7-0 table 141),
using
> > > NCEP-opn
> > > > > > > > >
1:0:d=15062500:MSLMA:kpds5=129:kpds6=200:kpds7=0:TR=10:
> > > > > > > > > P1=0:P2=0:TimeU=1:atmos
> > > > > > > > > col:anl:NAve=0
> > > > > > > > >
2:129684:d=15062500:KX:kpds5=133:kpds6=200:kpds7=0:TR=10:
> > > > > > > > > P1=0:P2=0:TimeU=1:atmos
> > > > > > > > > col:anl:NAve=0
> > > > > > > > >
3:259368:d=15062500:SX:kpds5=134:kpds6=200:kpds7=0:TR=10:
> > > > > > > > > P1=0:P2=0:TimeU=1:atmos
> > > > > > > > > col:anl:NAve=0
> > > > > > > > >
4:389052:d=15062500:MCONV:kpds5=135:kpds6=200:kpds7=0:
> > > > > > > > > TR=10:P1=0:P2=0:TimeU=1:atmos
> > > > > > > > > col:anl:NAve=0
> > > > > > > > >
5:518736:d=15062500:CWAT:kpds5=76:kpds6=200:kpds7=0:TR=
> > > > > > > > > 10:P1=0:P2=0:TimeU=1:atmos
> > > > > > > > > col:anl:NAve=0
> > > > > > > > >
6:648420:d=15062500:TSLSA:kpds5=137:kpds6=200:kpds7=0:
> > > > > > > > > TR=10:P1=0:P2=0:TimeU=1:atmos
> > > > > > > > > col:anl:NAve=0
> > > > > > > > >
> > > > > > > > > All of these 6 records appear to have different GRIB
codes.
> > > > Based
> > > > > on
> > > > > > > the
> > > > > > > > > question you asked, I was expecting to find multiple
> records
> > > with
> > > > > the
> > > > > > > > > *same* GRIB code but were different in some other
way (i.e.
> > > > > different
> > > > > > > > level
> > > > > > > > > type or different statistical processing type).
> > > > > > > > >
> > > > > > > > > You say these are all AOD, but the plot_data_plane
command
> > you
> > > > sent
> > > > > > > does
> > > > > > > > > *not* work for me:
> > > > > > > > >
> > > > > > > > > plot_data_plane out.grb aod.ps 'name="AOD";
level="A0";'
> > > > > > > > > ERROR : plot_data_plane -> trouble getting field
> > "name="AOD";
> > > > > > > > level="A0";"
> > > > > > > > > from file "out.grb"
> > > > > > > > >
> > > > > > > > > Please check to see if you have the $MET_GRIB_TABLES
> > > environment
> > > > > > > variable
> > > > > > > > > set. If so, please send me the table(s) to which it
> points.
> > > > > > > > >
> > > > > > > > > Thanks,
> > > > > > > > > John
> > > > > > > > >
> > > > > > > > >
> > > > > > > > > On Wed, Jan 17, 2018 at 10:35 AM, Partha
Bhattacharjee -
> NOAA
> > > > > > Affiliate
> > > > > > > > via
> > > > > > > > > RT <met_help at ucar.edu> wrote:
> > > > > > > > >
> > > > > > > > > >
> > > > > > > > > > Wed Jan 17 10:35:05 2018: Request 83618 was acted
upon.
> > > > > > > > > > Transaction: Ticket created by
> > partha.bhattacharjee at noaa.gov
> > > > > > > > > > Queue: met_help
> > > > > > > > > > Subject: Multiple records of AOD in a grib
file
> > > > > > > > > > Owner: Nobody
> > > > > > > > > > Requestors: partha.bhattacharjee at noaa.gov
> > > > > > > > > > Status: new
> > > > > > > > > > Ticket <URL: https://rt.rap.ucar.edu/rt/
> > > > > > > Ticket/Display.html?id=83618
> > > > > > > > >
> > > > > > > > > >
> > > > > > > > > >
> > > > > > > > > > Hi
> > > > > > > > > >
> > > > > > > > > > I have an input GRIB1 file with multiple records
of
> > > AOD,before
> > > > I
> > > > > > used
> > > > > > > > > > copygb to interpolate that into (-90:90;-180:180}
grid to
> > use
> > > > for
> > > > > > > > > > Grid-Stat. I am wondering how to read multiple AOD
> records
> > in
> > > > > that
> > > > > > > > file,
> > > > > > > > > is
> > > > > > > > > > there any Grib table specific number or template I
can
> use
> > ?
> > > > > > > > > >
> > > > > > > > > > MET able to read the interpolated file record when
I use
> > this
> > > > > > > command :
> > > > > > > > > > plot_data_plane out.grb aod.ps 'name="AOD";
> level="A0";'
> > > and
> > > > > > > > generated
> > > > > > > > > > the
> > > > > > > > > > attached figure.
> > > > > > > > > >
> > > > > > > > > > Attached are original and interpolated grib1
files.
> > > > > > > > > >
> > > > > > > > > > Thanks.
> > > > > > > > > > Partha
> > > > > > > > > >
> > > > > > > > > >
> > > > > > > > >
> > > > > > > > >
> > > > > > > >
> > > > > > > >
> > > > > > > > --
> > > > > > > > Partha S. Bhattacharjee, Ph.D.
> > > > > > > > I.M. Systems Group
> > > > > > > > Environmental Modeling Centre
> > > > > > > > NCWCP
> > > > > > > > 5830 University Research Court.
> > > > > > > > College Park, MD 20740
> > > > > > > >
> > > > > > > >
> > > > > > >
> > > > > > >
> > > > > >
> > > > > >
> > > > > > --
> > > > > > Partha S. Bhattacharjee, Ph.D.
> > > > > > I.M. Systems Group
> > > > > > Environmental Modeling Centre
> > > > > > NCWCP
> > > > > > 5830 University Research Court.
> > > > > > College Park, MD 20740
> > > > > >
> > > > > >
> > > > >
> > > > >
> > > >
> > > >
> > > > --
> > > > Partha S. Bhattacharjee, Ph.D.
> > > > I.M. Systems Group
> > > > Environmental Modeling Centre
> > > > NCWCP
> > > > 5830 University Research Court.
> > > > College Park, MD 20740
> > > >
> > > >
> > >
> > >
> >
> >
> > --
> > Partha S. Bhattacharjee, Ph.D.
> > I.M. Systems Group
> > Environmental Modeling Centre
> > NCWCP
> > 5830 University Research Court.
> > College Park, MD 20740
> >
> >
>
>
--
Partha S. Bhattacharjee, Ph.D.
I.M. Systems Group
Environmental Modeling Centre
NCWCP
5830 University Research Court.
College Park, MD 20740
------------------------------------------------
Subject: Multiple records of AOD in a grib file
From: Julie Prestopnik
Time: Mon Jan 22 10:43:37 2018
Great, Partha! Thank you for letting me know. I will close the
ticket.
Julie
On Mon, Jan 22, 2018 at 7:09 AM, Partha Bhattacharjee - NOAA Affiliate
via
RT <met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=83618 >
>
> Thanks Julie. I did it successfully in WCOSS as well. You may close
this
> ticket.
>
> Partha
>
> On Fri, Jan 19, 2018 at 10:34 AM, Julie Prestopnik via RT <
> met_help at ucar.edu
> > wrote:
>
> > Hi Partha.
> >
> > Yes, this same change will apply on WCOSS. This problem exists in
6.1 as
> > well. We will put in a fix in a patched version.
> >
> > Julie
> >
> > On Thu, Jan 18, 2018 at 9:40 PM, Partha Bhattacharjee - NOAA
Affiliate
> via
> > RT <met_help at ucar.edu> wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=83618 >
> > >
> > > Thanks Julie. I am able to modify accordingly and run Grid-Stat
> without
> > > any problem. Do I need to apply same change in WCOSS as well
(when it
> > comes
> > > back) ? And is this problem only with version 6 ?
> > >
> > > Partha
> > >
> > > On Thu, Jan 18, 2018 at 4:35 PM, Julie Prestopnik via RT <
> > > met_help at ucar.edu>
> > > wrote:
> > >
> > > > Hi Partha. Thank you for remembering the WCOSS outage and for
coyping
> > the
> > > > files over to theia! We were able to reproduce the problem
you are
> > > > experiencing.
> > > >
> > > > This error is occurring because of the empty string for the
units for
> > AOD
> > > > in the table file grib1_ncep_141_7.txt. Thank you for
bringing this
> to
> > > our
> > > > attention. We have added a ticket to resolve this issue. In
the
> > > meantime,
> > > > you can work around this issue by copying over the file
> > > > grib1_ncep_141_7.txt from your met-6.0/share/met/table_files/
> directory
> > > > into your current directory. Once you do that, please modify
the
> line:
> > > >
> > > > 129 141 7 -1 "AOD" "Aerosol Optical Depth" ""
> > > >
> > > > changing the empty string at the end to be "NA" or whatever
you'd
> like.
> > > > For example,
> > > >
> > > > 129 141 7 -1 "AOD" "Aerosol Optical Depth" "NA"
> > > >
> > > > Then, please set the environment variable MET_GRIB_TABLES to
point to
> > the
> > > > location of this new table. For example:
> > > >
> > > > export MET_GRIB_TABLES= /scratch4/NCEPDEV/global/
> > > > noscrub/Partha.Bhattacharjee/MET_Test/grib1_ncep_141_7.txt
> > > >
> > > > And, then please try running your script again. We were able
to get
> a
> > > > successful run once we made these changes. Please let us know
if you
> > > > encounter any problems or have any further questions.
> > > >
> > > > Thanks,
> > > > Julie
> > > >
> > > > On Wed, Jan 17, 2018 at 5:07 PM, Partha Bhattacharjee - NOAA
> Affiliate
> > > via
> > > > RT <met_help at ucar.edu> wrote:
> > > >
> > > > >
> > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=83618 >
> > > > >
> > > > > Thanks John. I realized WCOSS will be out next 2 days, so I
copied
> > the
> > > > > files over Theia
> > > > > at
/scratch4/NCEPDEV/global/noscrub/Partha.Bhattacharjee/MET_Test.
> > > > >
> > > > > Partha
> > > > >
> > > > > On Wed, Jan 17, 2018 at 5:43 PM, John Halley Gotway via RT <
> > > > > met_help at ucar.edu> wrote:
> > > > >
> > > > > > Partha,
> > > > > >
> > > > > > Unfortunately, I don't have access to WCOSS, but my
colleague,
> > Julie
> > > > > > Prestopnik, does. I'll ask her to log on and try to
replicate
> the
> > > > > behavior
> > > > > > you describe. And then we can debug it.
> > > > > >
> > > > > > Thanks,
> > > > > > John
> > > > > >
> > > > > > On Wed, Jan 17, 2018 at 2:34 PM, Partha Bhattacharjee -
NOAA
> > > Affiliate
> > > > > via
> > > > > > RT <met_help at ucar.edu> wrote:
> > > > > >
> > > > > > >
> > > > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=83618
> >
> > > > > > >
> > > > > > > Thanks John. Related to my earlier email, I am getting
error
> > > message
> > > > > when
> > > > > > > trying to use MET Series Stat. I am trying to compare
regridded
> > > > > forecast
> > > > > > > file and observational file but getting the following
error :
> > > > > > >
> > > > > > > GSL_RNG_TYPE=mt19937
> > > > > > > GSL_RNG_SEED=18446744073587810138
> > > > > > > terminate called after throwing an instance of
> 'std::logic_error'
> > > > > > > what(): basic_string::_S_construct NULL not valid
> > > > > > >
> > > > > > > I have attached the output.txt along with this email. My
files
> > are
> > > in
> > > > > > WCOSS
> > > > > > > at : /global/noscrub/Partha.Bhattacharjee/MET_Test.
Original
> > Grib
> > > > > files
> > > > > > > are
> > > > > > > in the Fcst directory : "AOD550.Grb.t00z.f00" and
interpolated
> > ones
> > > > are
> > > > > > > ones with "grb" at the end. I have used copygb to
interpolate
> > these
> > > > > input
> > > > > > > files which are at 384 X 190 to 360 X 180 to make them
same as
> > Obs
> > > > > files
> > > > > > > are. But not sure if the error is due to not matching
grids or
> > from
> > > > > > > something else.
> > > > > > >
> > > > > > > Thanks.
> > > > > > > Partha
> > > > > > >
> > > > > > > On Wed, Jan 17, 2018 at 4:26 PM, John Halley Gotway via
RT <
> > > > > > > met_help at ucar.edu> wrote:
> > > > > > >
> > > > > > > > Partha,
> > > > > > > >
> > > > > > > > Ah, OK. So your file contains 6 GRIB records, each of
which
> > has
> > > a
> > > > > > > > different GRIB code (129, 133, 134, 135, 76, and 137).
Those
> > are
> > > > the
> > > > > > > kpds5
> > > > > > > > values listed in the output of wgrib. Running wgrib
-V
> reveals
> > > > that
> > > > > 5
> > > > > > of
> > > > > > > > the 6 records are encoded as:
> > > > > > > > center 7 subcenter 0 process 82 Table 141 scan:
WE:SN
> > > winds(N/S)
> > > > > > > >
> > > > > > > > So this is NCEP GRIB1 table number 141, and that can
be found
> > > here:
> > > > > > > > met-6.1/share/met/table_files/grib1_ncep_141_7.txt
> > > > > > > >
> > > > > > > > Here's a selection of lines from that file:
> > > > > > > >
> > > > > > > > 129 141 7 -1 "AOD" "Aerosol Optical Depth" ""
> > > > > > > > 133 141 7 -1 "var133" "undefined" ""
> > > > > > > > 134 141 7 -1 "var134" "undefined" ""
> > > > > > > > 135 141 7 -1 "var135" "undefined" ""
> > > > > > > > 137 141 7 -1 "var137" "undefined" ""
> > > > > > > >
> > > > > > > > And this info is taken from this GRIB table:
> > > > > > > >
http://www.nco.ncep.noaa.gov/pmb/docs/on388/table2.html#
> > > > TABLE141
> > > > > > > >
> > > > > > > > So you could just refer to this data using those
"var133"
> (and
> > so
> > > > on)
> > > > > > > > string:
> > > > > > > >
> > > > > > > > plot_data_plane out.grb plot_rec1.ps
'name="var133";
> > > > > level="L0";'
> > > > > > > >
> > > > > > > > Alternatively, you could specifying the "GRIB1_code"
value
> > > instead
> > > > of
> > > > > > the
> > > > > > > > name:
> > > > > > > >
> > > > > > > > plot_data_plane out.grb plot_rec1.ps
'GRIB1_code=133;
> > > > > level="L0";'
> > > > > > > >
> > > > > > > > Both approaches would result in the same plot.
Presumably,
> > since
> > > > you
> > > > > > > have
> > > > > > > > data for codes 133, 134, 135, and 137, perhaps these
codes
> are
> > no
> > > > > > longer
> > > > > > > > "reserved for future use", as indicated on the NCO
website.
> > > > > > > >
> > > > > > > > Hope that helps clarify.
> > > > > > > >
> > > > > > > > Thanks,
> > > > > > > > John
> > > > > > > >
> > > > > > > > On Wed, Jan 17, 2018 at 12:05 PM, Partha Bhattacharjee
- NOAA
> > > > > Affiliate
> > > > > > > via
> > > > > > > > RT <met_help at ucar.edu> wrote:
> > > > > > > >
> > > > > > > > >
> > > > > > > > > <URL: https://rt.rap.ucar.edu/rt/
> > Ticket/Display.html?id=83618
> > > >
> > > > > > > > >
> > > > > > > > > John,
> > > > > > > > >
> > > > > > > > > That is strange. When I did wgrib to out.grib I am
getting
> > the
> > > > > > > following
> > > > > > > > :
> > > > > > > > >
> > > > > > > > > wgrib out.grb
> > > > > > > > >
1:0:d=15062500:AOD:kpds5=129:kpds6=200:kpds7=0:TR=10:P1=0:
> > > > > > > > > P2=0:TimeU=1:atmos
> > > > > > > > > col:anl:NAve=0
> > > > > > > > >
2:129684:d=15062500:var133:kpds5=133:kpds6=200:kpds7=0:
> > > > > > > > > TR=10:P1=0:P2=0:TimeU=1:atmos
> > > > > > > > > col:anl:NAve=0
> > > > > > > > >
3:259368:d=15062500:var134:kpds5=134:kpds6=200:kpds7=0:
> > > > > > > > > TR=10:P1=0:P2=0:TimeU=1:atmos
> > > > > > > > > col:anl:NAve=0
> > > > > > > > >
4:389052:d=15062500:var135:kpds5=135:kpds6=200:kpds7=0:
> > > > > > > > > TR=10:P1=0:P2=0:TimeU=1:atmos
> > > > > > > > > col:anl:NAve=0
> > > > > > > > >
5:518736:d=15062500:CWAT:kpds5=76:kpds6=200:kpds7=0:TR=
> > > > > > > > > 10:P1=0:P2=0:TimeU=1:atmos
> > > > > > > > > col:anl:NAve=0
> > > > > > > > >
6:648420:d=15062500:var137:kpds5=137:kpds6=200:kpds7=0:
> > > > > > > > > TR=10:P1=0:P2=0:TimeU=1:atmos
> > > > > > > > > col:anl:NAve=0
> > > > > > > > >
> > > > > > > > > and I am running in wcoss and wgrib points to
> > > > > > > > > which wgrib
> > > > > > > > > alias wgrib='/nwprod/util/exec/wgrib'
> > > > > > > > > /nwprod/util/exec/wgrib
> > > > > > > > >
> > > > > > > > > So, I used "AOD" in the first record to plot like :
> > > > > > > > >
> > > > > > > > > plot_data_plane out.grb aod.ps 'name="AOD";
level="L0";'
> > -v 6
> > > > > > > > > DEBUG 1: Opening data file: out.grb
> > > > > > > > > DEBUG 4:
Met2dDataFileFactory::new_met_2d_data_file() ->
> > > created
> > > > > new
> > > > > > > > > Met2dDataFile object of type "FileType_Gb1".
> > > > > > > > > DEBUG 4: VarInfoFactory::new_var_info() -> created
new
> > VarInfo
> > > > > object
> > > > > > > of
> > > > > > > > > type "FileType_Gb1".
> > > > > > > > > DEBUG 3: MetGrib1DataFile::data_plane_array() ->
Found
> range
> > > > match
> > > > > > for
> > > > > > > > > VarInfo "AOD/L0" in GRIB record 1 of GRIB file
"out.grb".
> > > > > > > > > DEBUG 3: MetGrib1DataFile::data_plane_array() ->
Found 1
> > GRIB
> > > > > > records
> > > > > > > > > matching VarInfo "AOD/L0" in GRIB file "out.grb".
> > > > > > > > > DEBUG 1: Creating postscript file: aod.ps
> > > > > > > > > DEBUG 4: draw_map() -> Plotting map data file:
> > > > > > > > > /awc/noscrub/Julie.Prestopnik/
> met/6.0/share/met/map/country_
> > > data
> > > > > > > > > DEBUG 4: draw_map() -> Plotting map data file:
> > > > > > > > >
/awc/noscrub/Julie.Prestopnik/met/6.0/share/met/map/usa_
> > > > state_data
> > > > > > > > > DEBUG 4: draw_map() -> Plotting map data file:
> > > > > > > > >
/awc/noscrub/Julie.Prestopnik/met/6.0/share/met/map/major_
> > > > > lakes_data
> > > > > > > > >
> > > > > > > > > I do not think I have $MET_GRIB_TABLES environment
variable
> > > set.
> > > > > > > > >
> > > > > > > > > Partha
> > > > > > > > >
> > > > > > > > > On Wed, Jan 17, 2018 at 1:33 PM, John Halley Gotway
via RT
> <
> > > > > > > > > met_help at ucar.edu> wrote:
> > > > > > > > >
> > > > > > > > > > Partha,
> > > > > > > > > >
> > > > > > > > > > I'm confused, when I run wgrib on the data you
sent, I
> > don't
> > > > see
> > > > > > any
> > > > > > > > > > mention of AOD. But I do get a warning message
about an
> > > > > undefined
> > > > > > > > > > parameter table:
> > > > > > > > > >
> > > > > > > > > > [johnhg at number5]% wgrib out.grb
> > > > > > > > > >
> > > > > > > > > > Undefined parameter table (center 7-0 table 141),
using
> > > > NCEP-opn
> > > > > > > > > >
1:0:d=15062500:MSLMA:kpds5=129:kpds6=200:kpds7=0:TR=10:
> > > > > > > > > > P1=0:P2=0:TimeU=1:atmos
> > > > > > > > > > col:anl:NAve=0
> > > > > > > > > > 2:129684:d=15062500:KX:kpds5=
> 133:kpds6=200:kpds7=0:TR=10:
> > > > > > > > > > P1=0:P2=0:TimeU=1:atmos
> > > > > > > > > > col:anl:NAve=0
> > > > > > > > > > 3:259368:d=15062500:SX:kpds5=
> 134:kpds6=200:kpds7=0:TR=10:
> > > > > > > > > > P1=0:P2=0:TimeU=1:atmos
> > > > > > > > > > col:anl:NAve=0
> > > > > > > > > >
4:389052:d=15062500:MCONV:kpds5=135:kpds6=200:kpds7=0:
> > > > > > > > > > TR=10:P1=0:P2=0:TimeU=1:atmos
> > > > > > > > > > col:anl:NAve=0
> > > > > > > > > >
5:518736:d=15062500:CWAT:kpds5=76:kpds6=200:kpds7=0:TR=
> > > > > > > > > > 10:P1=0:P2=0:TimeU=1:atmos
> > > > > > > > > > col:anl:NAve=0
> > > > > > > > > >
6:648420:d=15062500:TSLSA:kpds5=137:kpds6=200:kpds7=0:
> > > > > > > > > > TR=10:P1=0:P2=0:TimeU=1:atmos
> > > > > > > > > > col:anl:NAve=0
> > > > > > > > > >
> > > > > > > > > > All of these 6 records appear to have different
GRIB
> codes.
> > > > > Based
> > > > > > on
> > > > > > > > the
> > > > > > > > > > question you asked, I was expecting to find
multiple
> > records
> > > > with
> > > > > > the
> > > > > > > > > > *same* GRIB code but were different in some other
way
> (i.e.
> > > > > > different
> > > > > > > > > level
> > > > > > > > > > type or different statistical processing type).
> > > > > > > > > >
> > > > > > > > > > You say these are all AOD, but the plot_data_plane
> command
> > > you
> > > > > sent
> > > > > > > > does
> > > > > > > > > > *not* work for me:
> > > > > > > > > >
> > > > > > > > > > plot_data_plane out.grb aod.ps 'name="AOD";
> level="A0";'
> > > > > > > > > > ERROR : plot_data_plane -> trouble getting field
> > > "name="AOD";
> > > > > > > > > level="A0";"
> > > > > > > > > > from file "out.grb"
> > > > > > > > > >
> > > > > > > > > > Please check to see if you have the
$MET_GRIB_TABLES
> > > > environment
> > > > > > > > variable
> > > > > > > > > > set. If so, please send me the table(s) to which
it
> > points.
> > > > > > > > > >
> > > > > > > > > > Thanks,
> > > > > > > > > > John
> > > > > > > > > >
> > > > > > > > > >
> > > > > > > > > > On Wed, Jan 17, 2018 at 10:35 AM, Partha
Bhattacharjee -
> > NOAA
> > > > > > > Affiliate
> > > > > > > > > via
> > > > > > > > > > RT <met_help at ucar.edu> wrote:
> > > > > > > > > >
> > > > > > > > > > >
> > > > > > > > > > > Wed Jan 17 10:35:05 2018: Request 83618 was
acted upon.
> > > > > > > > > > > Transaction: Ticket created by
> > > partha.bhattacharjee at noaa.gov
> > > > > > > > > > > Queue: met_help
> > > > > > > > > > > Subject: Multiple records of AOD in a grib
file
> > > > > > > > > > > Owner: Nobody
> > > > > > > > > > > Requestors: partha.bhattacharjee at noaa.gov
> > > > > > > > > > > Status: new
> > > > > > > > > > > Ticket <URL: https://rt.rap.ucar.edu/rt/
> > > > > > > > Ticket/Display.html?id=83618
> > > > > > > > > >
> > > > > > > > > > >
> > > > > > > > > > >
> > > > > > > > > > > Hi
> > > > > > > > > > >
> > > > > > > > > > > I have an input GRIB1 file with multiple records
of
> > > > AOD,before
> > > > > I
> > > > > > > used
> > > > > > > > > > > copygb to interpolate that into (-90:90;-
180:180} grid
> to
> > > use
> > > > > for
> > > > > > > > > > > Grid-Stat. I am wondering how to read multiple
AOD
> > records
> > > in
> > > > > > that
> > > > > > > > > file,
> > > > > > > > > > is
> > > > > > > > > > > there any Grib table specific number or template
I can
> > use
> > > ?
> > > > > > > > > > >
> > > > > > > > > > > MET able to read the interpolated file record
when I
> use
> > > this
> > > > > > > > command :
> > > > > > > > > > > plot_data_plane out.grb aod.ps 'name="AOD";
> > level="A0";'
> > > > and
> > > > > > > > > generated
> > > > > > > > > > > the
> > > > > > > > > > > attached figure.
> > > > > > > > > > >
> > > > > > > > > > > Attached are original and interpolated grib1
files.
> > > > > > > > > > >
> > > > > > > > > > > Thanks.
> > > > > > > > > > > Partha
> > > > > > > > > > >
> > > > > > > > > > >
> > > > > > > > > >
> > > > > > > > > >
> > > > > > > > >
> > > > > > > > >
> > > > > > > > > --
> > > > > > > > > Partha S. Bhattacharjee, Ph.D.
> > > > > > > > > I.M. Systems Group
> > > > > > > > > Environmental Modeling Centre
> > > > > > > > > NCWCP
> > > > > > > > > 5830 University Research Court.
> > > > > > > > > College Park, MD 20740
> > > > > > > > >
> > > > > > > > >
> > > > > > > >
> > > > > > > >
> > > > > > >
> > > > > > >
> > > > > > > --
> > > > > > > Partha S. Bhattacharjee, Ph.D.
> > > > > > > I.M. Systems Group
> > > > > > > Environmental Modeling Centre
> > > > > > > NCWCP
> > > > > > > 5830 University Research Court.
> > > > > > > College Park, MD 20740
> > > > > > >
> > > > > > >
> > > > > >
> > > > > >
> > > > >
> > > > >
> > > > > --
> > > > > Partha S. Bhattacharjee, Ph.D.
> > > > > I.M. Systems Group
> > > > > Environmental Modeling Centre
> > > > > NCWCP
> > > > > 5830 University Research Court.
> > > > > College Park, MD 20740
> > > > >
> > > > >
> > > >
> > > >
> > >
> > >
> > > --
> > > Partha S. Bhattacharjee, Ph.D.
> > > I.M. Systems Group
> > > Environmental Modeling Centre
> > > NCWCP
> > > 5830 University Research Court.
> > > College Park, MD 20740
> > >
> > >
> >
> >
>
>
> --
> Partha S. Bhattacharjee, Ph.D.
> I.M. Systems Group
> Environmental Modeling Centre
> NCWCP
> 5830 University Research Court.
> College Park, MD 20740
>
>
------------------------------------------------
More information about the Met_help
mailing list