[Met_help] [rt.rap.ucar.edu #77365] History for Settings for config file for GriB2 data
John Halley Gotway via RT
met_help at ucar.edu
Fri Aug 12 11:45:00 MDT 2016
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
Hi,
I'm trying to get point_stat to verify a GriB2 model file but I'm having trouble figuring out the right config file settings. I've never tried using MET on GriB2 data before.
I think I have the "level" set wrong in the "fcst" and "obs" definitions. I know the parm codes of the fields and when I run MET it agrees with the lookup table as to which variable I am seeking.
Specifically I am trying to verify surface and skin temperature fields against surface temperature observations that I have converted using "ascii2nc". The obs file worked great for a similar test using GriB1 data so I believe I don't need to make changes to those config file settings for the obs. It's the model I think I have incorrect.
Later there's a mask setting to use the land/water flag field from the model file to subset the verification into land only and water only. Again, I did that for a GriB1 case which worked.
The GRIB2_parm numbers all match the lookup table, and the output when I run MET suggests it is looking for the correct fields that I want. It's just not finding anything matching in the model file for "L0". See below for the output from MET.
I've placed my GriB2 model and NetCDF observation files, plus my config file, in "sittel_data" on your FTP server under "met_help".
Thanks,
Matt
DEBUG 1: Default Config File: /home/qcteam/met-5.1/share/met/config/PointStatConfig_default
DEBUG 1: User Config File: /home/qcteam/met-5.1/cfg/sfctmp/point_stat_GALWEM_OBS.cfg
DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() -> created new Met2dDataFile object of type "FileType_Gb2".
DEBUG 4: Switching the GRIB2 radius of the earth value of 6371.23 km to 6371.2 km for internal consistency.
DEBUG 4:
DEBUG 4: Latitude/Longitude Grid Data:
DEBUG 4: lat_ll: -89.9219
DEBUG 4: lon_ll: -0.117187
DEBUG 4: delta_lat: 0.15625
DEBUG 4: delta_lon: 0.234375
DEBUG 4: Nlat: 1152
DEBUG 4: Nlon: 1536
DEBUG 4:
DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo object of type "FileType_Gb2".
DEBUG 4: parse_conf_sid_exc() -> Station ID exclusion list contains 0 entries.
DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo object of type "FileType_Gb2".
DEBUG 4: parse_conf_sid_exc() -> Station ID exclusion list contains 0 entries.
GSL_RNG_TYPE=mt19937
GSL_RNG_SEED=1
DEBUG 1: Forecast File: /home/hildebre/sfctmp/um/PS.557WW_SC.U_DI.F_GP.GALWEM-GD_GR.C17KM_AR.GLOBAL_DD.20160713_CY.00_FH.012_DF.GR2
DEBUG 1: Observation File: /home/qcteam/met-5.1/temp/2016071312.sfc.nc
DEBUG 2:
DEBUG 2: --------------------------------------------------------------------------------
DEBUG 2:
DEBUG 2: Reading data for TMP/L0.
WARNING:
WARNING: MetGrib2DataFile::data_plane_array() - No matching records found for 'TMP/L0'
WARNING:
WARNING:
WARNING: process_fcst_climo_files() -> no fields matching TMP/L0 found in file: /home/hildebre/sfctmp/um/PS.557WW_SC.U_DI.F_GP.GALWEM-GD_GR.C17KM_AR.GLOBAL_DD.20160713_CY.00_FH.012_DF.GR2
WARNING:
DEBUG 2:
DEBUG 2: --------------------------------------------------------------------------------
DEBUG 2:
DEBUG 2: Reading data for SKINT/L0.
WARNING:
WARNING: MetGrib2DataFile::data_plane_array() - No matching records found for 'SKINT/L0'
WARNING:
WARNING:
WARNING: process_fcst_climo_files() -> no fields matching SKINT/L0 found in file: /home/hildebre/sfctmp/um/PS.557WW_SC.U_DI.F_GP.GALWEM-GD_GR.C17KM_AR.GLOBAL_DD.20160713_CY.00_FH.012_DF.GR2
WARNING:
DEBUG 2:
DEBUG 2: --------------------------------------------------------------------------------
DEBUG 2:
DEBUG 2: Searching 71616 observations from 11936 messages.
DEBUG 2:
DEBUG 2: --------------------------------------------------------------------------------
DEBUG 2:
DEBUG 1: Output file: /home/qcteam/met-5.1/stats/sfctmp/GALWEM_OBS/point_stat_-02-46-39L_19700101_000000V.stat
DEBUG 1: Output file: /home/qcteam/met-5.1/stats/sfctmp/GALWEM_OBS/point_stat_-02-46-39L_19700101_000000V_cnt.txt
// SIGNED //
Matthew C. Sittel
University Corporation for Atmospheric Research
16 WS/WXN, 557 WW, Offutt AFB, NE
matthew.sittel.ctr at us.af.mil
----------------------------------------------------------------
Complete Ticket History
----------------------------------------------------------------
Subject: Settings for config file for GriB2 data
From: John Halley Gotway
Time: Thu Jul 28 12:05:16 2016
Hi Matt,
Here's the wgrib2 output for that file:
1:0:d=2016071300:TMP:2 m above ground:12 hour fcst:
2:929176:d=2016071300:LAND:surface:12 hour fcst:
3:976581:d=2016071300:SKINT:surface:12 hour fcst:
Looks like you have 2-m temperature, a land/sea mask, and skin
temperature. When I'm having trouble reading data, I always use
plot_data_plane to test things out.
2-m Temperature is easy. Just use 'name="TMP"; level="Z2";'
plot_data_plane
PS.557WW_SC.U_DI.F_GP.GALWEM-
GD_GR.C17KM_AR.GLOBAL_DD.20160713_CY.00_FH.012_DF.GR2
tmp_z2.ps 'name="TMP"; level="Z2";'
However, LAND and SKINT aren't working very well. I ran it through
the
debugger and see that the problem is with the level value stored in
the
GRIB2 file. The debugger lists the level value as -2147483647.
One way to get around this is to just use the record number instead.
So
the following name/level combinations work:
- name="TMP"; level="R1";
- name="LAND"; level="R2";
- name="SKINT"; level="R3";
But that isn't very satisfying! I see in the code how we can fix this
behavior so that 'name="SKINT"; level="L0";' will work.
I'll commit that change for version 5.2 which will likely be released
next
week.
Are you OK using "R3" for SKINT for now and then you could use "L0"
once
5.2 is out?
Thanks,
John
On Thu, Jul 28, 2016 at 7:53 AM, SITTEL, MATTHEW C CTR USAF AFWA 16
WS/WXN
via RT <met_help at ucar.edu> wrote:
>
> Thu Jul 28 07:53:14 2016: Request 77365 was acted upon.
> Transaction: Ticket created by matthew.sittel.ctr at us.af.mil
> Queue: met_help
> Subject: Settings for config file for GriB2 data
> Owner: Nobody
> Requestors: matthew.sittel.ctr at us.af.mil
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=77365 >
>
>
> Hi,
>
> I'm trying to get point_stat to verify a GriB2 model file but I'm
having
> trouble figuring out the right config file settings. I've never
tried
> using MET on GriB2 data before.
>
> I think I have the "level" set wrong in the "fcst" and "obs"
definitions.
> I know the parm codes of the fields and when I run MET it agrees
with the
> lookup table as to which variable I am seeking.
>
> Specifically I am trying to verify surface and skin temperature
fields
> against surface temperature observations that I have converted using
> "ascii2nc". The obs file worked great for a similar test using GriB1
data
> so I believe I don't need to make changes to those config file
settings for
> the obs. It's the model I think I have incorrect.
>
> Later there's a mask setting to use the land/water flag field from
the
> model file to subset the verification into land only and water only.
> Again, I did that for a GriB1 case which worked.
>
> The GRIB2_parm numbers all match the lookup table, and the output
when I
> run MET suggests it is looking for the correct fields that I want.
It's
> just not finding anything matching in the model file for "L0". See
below
> for the output from MET.
>
> I've placed my GriB2 model and NetCDF observation files, plus my
config
> file, in "sittel_data" on your FTP server under "met_help".
>
> Thanks,
> Matt
>
> DEBUG 1: Default Config File:
> /home/qcteam/met-5.1/share/met/config/PointStatConfig_default
> DEBUG 1: User Config File:
> /home/qcteam/met-5.1/cfg/sfctmp/point_stat_GALWEM_OBS.cfg
> DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() -> created new
> Met2dDataFile object of type "FileType_Gb2".
> DEBUG 4: Switching the GRIB2 radius of the earth value of 6371.23 km
to
> 6371.2 km for internal consistency.
> DEBUG 4:
> DEBUG 4: Latitude/Longitude Grid Data:
> DEBUG 4: lat_ll: -89.9219
> DEBUG 4: lon_ll: -0.117187
> DEBUG 4: delta_lat: 0.15625
> DEBUG 4: delta_lon: 0.234375
> DEBUG 4: Nlat: 1152
> DEBUG 4: Nlon: 1536
> DEBUG 4:
> DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo
object of
> type "FileType_Gb2".
> DEBUG 4: parse_conf_sid_exc() -> Station ID exclusion list contains
0
> entries.
> DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo
object of
> type "FileType_Gb2".
> DEBUG 4: parse_conf_sid_exc() -> Station ID exclusion list contains
0
> entries.
> GSL_RNG_TYPE=mt19937
> GSL_RNG_SEED=1
> DEBUG 1: Forecast File:
> /home/hildebre/sfctmp/um/PS.557WW_SC.U_DI.F_GP.GALWEM-
GD_GR.C17KM_AR.GLOBAL_DD.20160713_CY.00_FH.012_DF.GR2
> DEBUG 1: Observation File: /home/qcteam/met-
5.1/temp/2016071312.sfc.nc
> DEBUG 2:
> DEBUG 2:
>
--------------------------------------------------------------------------------
> DEBUG 2:
> DEBUG 2: Reading data for TMP/L0.
> WARNING:
> WARNING: MetGrib2DataFile::data_plane_array() - No matching records
found
> for 'TMP/L0'
> WARNING:
> WARNING:
> WARNING: process_fcst_climo_files() -> no fields matching TMP/L0
found in
> file:
> /home/hildebre/sfctmp/um/PS.557WW_SC.U_DI.F_GP.GALWEM-
GD_GR.C17KM_AR.GLOBAL_DD.20160713_CY.00_FH.012_DF.GR2
> WARNING:
> DEBUG 2:
> DEBUG 2:
>
--------------------------------------------------------------------------------
> DEBUG 2:
> DEBUG 2: Reading data for SKINT/L0.
> WARNING:
> WARNING: MetGrib2DataFile::data_plane_array() - No matching records
found
> for 'SKINT/L0'
> WARNING:
> WARNING:
> WARNING: process_fcst_climo_files() -> no fields matching SKINT/L0
found
> in file:
> /home/hildebre/sfctmp/um/PS.557WW_SC.U_DI.F_GP.GALWEM-
GD_GR.C17KM_AR.GLOBAL_DD.20160713_CY.00_FH.012_DF.GR2
> WARNING:
> DEBUG 2:
> DEBUG 2:
>
--------------------------------------------------------------------------------
> DEBUG 2:
> DEBUG 2: Searching 71616 observations from 11936 messages.
> DEBUG 2:
> DEBUG 2:
>
--------------------------------------------------------------------------------
> DEBUG 2:
> DEBUG 1: Output file:
> /home/qcteam/met-5.1/stats/sfctmp/GALWEM_OBS/point_stat_-02-46-
39L_19700101_000000V.stat
> DEBUG 1: Output file:
> /home/qcteam/met-5.1/stats/sfctmp/GALWEM_OBS/point_stat_-02-46-
39L_19700101_000000V_cnt.txt
>
> // SIGNED //
> Matthew C. Sittel
> University Corporation for Atmospheric Research
> 16 WS/WXN, 557 WW, Offutt AFB, NE
> matthew.sittel.ctr at us.af.mil
>
>
>
>
------------------------------------------------
Subject: RE: [rt.rap.ucar.edu #77365] Settings for config file for GriB2 data
From: SITTEL, MATTHEW C CTR USAF AFWA 16 WS/WXN
Time: Thu Jul 28 13:19:04 2016
Hey John,
Thanks for the info. Hard coding the "R1" et al. fields will be okay
for this project, but obviously in production I'm hesitant to do that
since I can't be sure fields won't be added or deleted which might
change any field's relative position in the GriB output.
As long as there is a workaround I don't think you have to add that to
5.2 unless you want to. The GriB2 data I sent was made in-house and
probably doesn't have any adherence to guidelines for variable
definition in the header. I'm not sure. If that is an issue I'll be
sure to let the developers know. Anything specific I should have them
set for future reference?
Thanks,
Matt
-----Original Message-----
From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
Sent: Thursday, July 28, 2016 1:05 PM
To: SITTEL, MATTHEW C CTR USAF AFMC AFLCMC/HBAW-OL
<matthew.sittel.ctr at us.af.mil>
Cc: CRAIG, ROBERT J GS-12 USAF ACC 16 WS/WXN
<robert.craig.2 at us.af.mil>
Subject: Re: [rt.rap.ucar.edu #77365] Settings for config file for
GriB2 data
Hi Matt,
Here's the wgrib2 output for that file:
1:0:d=2016071300:TMP:2 m above ground:12 hour fcst:
2:929176:d=2016071300:LAND:surface:12 hour fcst:
3:976581:d=2016071300:SKINT:surface:12 hour fcst:
Looks like you have 2-m temperature, a land/sea mask, and skin
temperature. When I'm having trouble reading data, I always use
plot_data_plane to test things out.
2-m Temperature is easy. Just use 'name="TMP"; level="Z2";'
plot_data_plane
PS.557WW_SC.U_DI.F_GP.GALWEM-
GD_GR.C17KM_AR.GLOBAL_DD.20160713_CY.00_FH.012_DF.GR2
tmp_z2.ps 'name="TMP"; level="Z2";'
However, LAND and SKINT aren't working very well. I ran it through
the debugger and see that the problem is with the level value stored
in the
GRIB2 file. The debugger lists the level value as -2147483647.
One way to get around this is to just use the record number instead.
So the following name/level combinations work:
- name="TMP"; level="R1";
- name="LAND"; level="R2";
- name="SKINT"; level="R3";
But that isn't very satisfying! I see in the code how we can fix this
behavior so that 'name="SKINT"; level="L0";' will work.
I'll commit that change for version 5.2 which will likely be released
next week.
Are you OK using "R3" for SKINT for now and then you could use "L0"
once
5.2 is out?
Thanks,
John
On Thu, Jul 28, 2016 at 7:53 AM, SITTEL, MATTHEW C CTR USAF AFWA 16
WS/WXN via RT <met_help at ucar.edu> wrote:
>
> Thu Jul 28 07:53:14 2016: Request 77365 was acted upon.
> Transaction: Ticket created by matthew.sittel.ctr at us.af.mil
> Queue: met_help
> Subject: Settings for config file for GriB2 data
> Owner: Nobody
> Requestors: matthew.sittel.ctr at us.af.mil
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=77365
> >
>
>
> Hi,
>
> I'm trying to get point_stat to verify a GriB2 model file but I'm
> having trouble figuring out the right config file settings. I've
> never tried using MET on GriB2 data before.
>
> I think I have the "level" set wrong in the "fcst" and "obs"
definitions.
> I know the parm codes of the fields and when I run MET it agrees
with
> the lookup table as to which variable I am seeking.
>
> Specifically I am trying to verify surface and skin temperature
fields
> against surface temperature observations that I have converted using
> "ascii2nc". The obs file worked great for a similar test using GriB1
> data so I believe I don't need to make changes to those config file
> settings for the obs. It's the model I think I have incorrect.
>
> Later there's a mask setting to use the land/water flag field from
the
> model file to subset the verification into land only and water only.
> Again, I did that for a GriB1 case which worked.
>
> The GRIB2_parm numbers all match the lookup table, and the output
when
> I run MET suggests it is looking for the correct fields that I want.
> It's just not finding anything matching in the model file for "L0".
> See below for the output from MET.
>
> I've placed my GriB2 model and NetCDF observation files, plus my
> config file, in "sittel_data" on your FTP server under "met_help".
>
> Thanks,
> Matt
>
> DEBUG 1: Default Config File:
> /home/qcteam/met-5.1/share/met/config/PointStatConfig_default
> DEBUG 1: User Config File:
> /home/qcteam/met-5.1/cfg/sfctmp/point_stat_GALWEM_OBS.cfg
> DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() -> created new
> Met2dDataFile object of type "FileType_Gb2".
> DEBUG 4: Switching the GRIB2 radius of the earth value of 6371.23 km
> to
> 6371.2 km for internal consistency.
> DEBUG 4:
> DEBUG 4: Latitude/Longitude Grid Data:
> DEBUG 4: lat_ll: -89.9219
> DEBUG 4: lon_ll: -0.117187
> DEBUG 4: delta_lat: 0.15625
> DEBUG 4: delta_lon: 0.234375
> DEBUG 4: Nlat: 1152
> DEBUG 4: Nlon: 1536
> DEBUG 4:
> DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo
object
> of type "FileType_Gb2".
> DEBUG 4: parse_conf_sid_exc() -> Station ID exclusion list contains
0
> entries.
> DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo
object
> of type "FileType_Gb2".
> DEBUG 4: parse_conf_sid_exc() -> Station ID exclusion list contains
0
> entries.
> GSL_RNG_TYPE=mt19937
> GSL_RNG_SEED=1
> DEBUG 1: Forecast File:
> /home/hildebre/sfctmp/um/PS.557WW_SC.U_DI.F_GP.GALWEM-
GD_GR.C17KM_AR.G
> LOBAL_DD.20160713_CY.00_FH.012_DF.GR2
> DEBUG 1: Observation File: /home/qcteam/met-
5.1/temp/2016071312.sfc.nc
> DEBUG 2:
> DEBUG 2:
>
----------------------------------------------------------------------
> ----------
> DEBUG 2:
> DEBUG 2: Reading data for TMP/L0.
> WARNING:
> WARNING: MetGrib2DataFile::data_plane_array() - No matching records
> found for 'TMP/L0'
> WARNING:
> WARNING:
> WARNING: process_fcst_climo_files() -> no fields matching TMP/L0
found
> in
> file:
> /home/hildebre/sfctmp/um/PS.557WW_SC.U_DI.F_GP.GALWEM-
GD_GR.C17KM_AR.G
> LOBAL_DD.20160713_CY.00_FH.012_DF.GR2
> WARNING:
> DEBUG 2:
> DEBUG 2:
>
----------------------------------------------------------------------
> ----------
> DEBUG 2:
> DEBUG 2: Reading data for SKINT/L0.
> WARNING:
> WARNING: MetGrib2DataFile::data_plane_array() - No matching records
> found for 'SKINT/L0'
> WARNING:
> WARNING:
> WARNING: process_fcst_climo_files() -> no fields matching SKINT/L0
> found in file:
> /home/hildebre/sfctmp/um/PS.557WW_SC.U_DI.F_GP.GALWEM-
GD_GR.C17KM_AR.G
> LOBAL_DD.20160713_CY.00_FH.012_DF.GR2
> WARNING:
> DEBUG 2:
> DEBUG 2:
>
----------------------------------------------------------------------
> ----------
> DEBUG 2:
> DEBUG 2: Searching 71616 observations from 11936 messages.
> DEBUG 2:
> DEBUG 2:
>
----------------------------------------------------------------------
> ----------
> DEBUG 2:
> DEBUG 1: Output file:
> /home/qcteam/met-5.1/stats/sfctmp/GALWEM_OBS/point_stat_-02-46-
39L_197
> 00101_000000V.stat
> DEBUG 1: Output file:
> /home/qcteam/met-5.1/stats/sfctmp/GALWEM_OBS/point_stat_-02-46-
39L_197
> 00101_000000V_cnt.txt
>
> // SIGNED //
> Matthew C. Sittel
> University Corporation for Atmospheric Research
> 16 WS/WXN, 557 WW, Offutt AFB, NE
> matthew.sittel.ctr at us.af.mil
>
>
>
>
------------------------------------------------
Subject: Settings for config file for GriB2 data
From: John Halley Gotway
Time: Fri Jul 29 11:50:16 2016
Matt,
For LAND and SKINT, here's whats going on...
MET reads the level type value from the GRIB2 PDS. For both of these,
the
level type is set to a value of 1. That means "Ground or Water
Surface",
based on this table:
http://www.nco.ncep.noaa.gov/pmb/docs/grib2/grib2_table4-5.shtml
MET was proceeding to parse the actual level values from octets 25-28
and
31-34:
http://www.nco.ncep.noaa.gov/pmb/docs/grib2/grib2_temp4-0.shtml
In your data, the level values were set to that odd number
-2147483647.
The fix I added for MET version 5.2 is to check for a level type of 1
(Ground of Water Surface) and in that case set the level values to 0
instead of actually reading them from the PDS.
I don't really know what's "correct" here. Perhaps setting the level
type
to 1 should be sufficient and you shouldn't have to worry about
setting the
level values themselves? If you are going to set them, 0 seems like a
good
choice.
We haven't had this error with other GRIB2 datasets so presumably, the
level values have been set to 0. That would make your GRIB2 data more
consistent with other data producers.
John
On Thu, Jul 28, 2016 at 1:19 PM, SITTEL, MATTHEW C CTR USAF AFWA 16
WS/WXN
via RT <met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=77365 >
>
> Hey John,
>
> Thanks for the info. Hard coding the "R1" et al. fields will be
okay for
> this project, but obviously in production I'm hesitant to do that
since I
> can't be sure fields won't be added or deleted which might change
any
> field's relative position in the GriB output.
>
> As long as there is a workaround I don't think you have to add that
to 5.2
> unless you want to. The GriB2 data I sent was made in-house and
probably
> doesn't have any adherence to guidelines for variable definition in
the
> header. I'm not sure. If that is an issue I'll be sure to let the
> developers know. Anything specific I should have them set for
future
> reference?
>
> Thanks,
> Matt
>
> -----Original Message-----
> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> Sent: Thursday, July 28, 2016 1:05 PM
> To: SITTEL, MATTHEW C CTR USAF AFMC AFLCMC/HBAW-OL <
> matthew.sittel.ctr at us.af.mil>
> Cc: CRAIG, ROBERT J GS-12 USAF ACC 16 WS/WXN
<robert.craig.2 at us.af.mil>
> Subject: Re: [rt.rap.ucar.edu #77365] Settings for config file for
GriB2
> data
>
> Hi Matt,
>
> Here's the wgrib2 output for that file:
> 1:0:d=2016071300:TMP:2 m above ground:12 hour fcst:
> 2:929176:d=2016071300:LAND:surface:12 hour fcst:
> 3:976581:d=2016071300:SKINT:surface:12 hour fcst:
>
> Looks like you have 2-m temperature, a land/sea mask, and skin
> temperature. When I'm having trouble reading data, I always use
> plot_data_plane to test things out.
>
> 2-m Temperature is easy. Just use 'name="TMP"; level="Z2";'
>
> plot_data_plane
>
> PS.557WW_SC.U_DI.F_GP.GALWEM-
GD_GR.C17KM_AR.GLOBAL_DD.20160713_CY.00_FH.012_DF.GR2
> tmp_z2.ps 'name="TMP"; level="Z2";'
>
> However, LAND and SKINT aren't working very well. I ran it through
the
> debugger and see that the problem is with the level value stored in
the
> GRIB2 file. The debugger lists the level value as -2147483647.
>
> One way to get around this is to just use the record number instead.
So
> the following name/level combinations work:
> - name="TMP"; level="R1";
> - name="LAND"; level="R2";
> - name="SKINT"; level="R3";
>
> But that isn't very satisfying! I see in the code how we can fix
this
> behavior so that 'name="SKINT"; level="L0";' will work.
>
> I'll commit that change for version 5.2 which will likely be
released next
> week.
>
> Are you OK using "R3" for SKINT for now and then you could use "L0"
once
> 5.2 is out?
>
> Thanks,
> John
>
>
>
>
> On Thu, Jul 28, 2016 at 7:53 AM, SITTEL, MATTHEW C CTR USAF AFWA 16
WS/WXN
> via RT <met_help at ucar.edu> wrote:
>
> >
> > Thu Jul 28 07:53:14 2016: Request 77365 was acted upon.
> > Transaction: Ticket created by matthew.sittel.ctr at us.af.mil
> > Queue: met_help
> > Subject: Settings for config file for GriB2 data
> > Owner: Nobody
> > Requestors: matthew.sittel.ctr at us.af.mil
> > Status: new
> > Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=77365
> > >
> >
> >
> > Hi,
> >
> > I'm trying to get point_stat to verify a GriB2 model file but I'm
> > having trouble figuring out the right config file settings. I've
> > never tried using MET on GriB2 data before.
> >
> > I think I have the "level" set wrong in the "fcst" and "obs"
definitions.
> > I know the parm codes of the fields and when I run MET it agrees
with
> > the lookup table as to which variable I am seeking.
> >
> > Specifically I am trying to verify surface and skin temperature
fields
> > against surface temperature observations that I have converted
using
> > "ascii2nc". The obs file worked great for a similar test using
GriB1
> > data so I believe I don't need to make changes to those config
file
> > settings for the obs. It's the model I think I have incorrect.
> >
> > Later there's a mask setting to use the land/water flag field from
the
> > model file to subset the verification into land only and water
only.
> > Again, I did that for a GriB1 case which worked.
> >
> > The GRIB2_parm numbers all match the lookup table, and the output
when
> > I run MET suggests it is looking for the correct fields that I
want.
> > It's just not finding anything matching in the model file for
"L0".
> > See below for the output from MET.
> >
> > I've placed my GriB2 model and NetCDF observation files, plus my
> > config file, in "sittel_data" on your FTP server under "met_help".
> >
> > Thanks,
> > Matt
> >
> > DEBUG 1: Default Config File:
> > /home/qcteam/met-5.1/share/met/config/PointStatConfig_default
> > DEBUG 1: User Config File:
> > /home/qcteam/met-5.1/cfg/sfctmp/point_stat_GALWEM_OBS.cfg
> > DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() -> created
new
> > Met2dDataFile object of type "FileType_Gb2".
> > DEBUG 4: Switching the GRIB2 radius of the earth value of 6371.23
km
> > to
> > 6371.2 km for internal consistency.
> > DEBUG 4:
> > DEBUG 4: Latitude/Longitude Grid Data:
> > DEBUG 4: lat_ll: -89.9219
> > DEBUG 4: lon_ll: -0.117187
> > DEBUG 4: delta_lat: 0.15625
> > DEBUG 4: delta_lon: 0.234375
> > DEBUG 4: Nlat: 1152
> > DEBUG 4: Nlon: 1536
> > DEBUG 4:
> > DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo
object
> > of type "FileType_Gb2".
> > DEBUG 4: parse_conf_sid_exc() -> Station ID exclusion list
contains 0
> > entries.
> > DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo
object
> > of type "FileType_Gb2".
> > DEBUG 4: parse_conf_sid_exc() -> Station ID exclusion list
contains 0
> > entries.
> > GSL_RNG_TYPE=mt19937
> > GSL_RNG_SEED=1
> > DEBUG 1: Forecast File:
> > /home/hildebre/sfctmp/um/PS.557WW_SC.U_DI.F_GP.GALWEM-
GD_GR.C17KM_AR.G
> > LOBAL_DD.20160713_CY.00_FH.012_DF.GR2
> > DEBUG 1: Observation File: /home/qcteam/met-
5.1/temp/2016071312.sfc.nc
> > DEBUG 2:
> > DEBUG 2:
> >
----------------------------------------------------------------------
> > ----------
> > DEBUG 2:
> > DEBUG 2: Reading data for TMP/L0.
> > WARNING:
> > WARNING: MetGrib2DataFile::data_plane_array() - No matching
records
> > found for 'TMP/L0'
> > WARNING:
> > WARNING:
> > WARNING: process_fcst_climo_files() -> no fields matching TMP/L0
found
> > in
> > file:
> > /home/hildebre/sfctmp/um/PS.557WW_SC.U_DI.F_GP.GALWEM-
GD_GR.C17KM_AR.G
> > LOBAL_DD.20160713_CY.00_FH.012_DF.GR2
> > WARNING:
> > DEBUG 2:
> > DEBUG 2:
> >
----------------------------------------------------------------------
> > ----------
> > DEBUG 2:
> > DEBUG 2: Reading data for SKINT/L0.
> > WARNING:
> > WARNING: MetGrib2DataFile::data_plane_array() - No matching
records
> > found for 'SKINT/L0'
> > WARNING:
> > WARNING:
> > WARNING: process_fcst_climo_files() -> no fields matching SKINT/L0
> > found in file:
> > /home/hildebre/sfctmp/um/PS.557WW_SC.U_DI.F_GP.GALWEM-
GD_GR.C17KM_AR.G
> > LOBAL_DD.20160713_CY.00_FH.012_DF.GR2
> > WARNING:
> > DEBUG 2:
> > DEBUG 2:
> >
----------------------------------------------------------------------
> > ----------
> > DEBUG 2:
> > DEBUG 2: Searching 71616 observations from 11936 messages.
> > DEBUG 2:
> > DEBUG 2:
> >
----------------------------------------------------------------------
> > ----------
> > DEBUG 2:
> > DEBUG 1: Output file:
> > /home/qcteam/met-5.1/stats/sfctmp/GALWEM_OBS/point_stat_-02-46-
39L_197
> > 00101_000000V.stat
> > DEBUG 1: Output file:
> > /home/qcteam/met-5.1/stats/sfctmp/GALWEM_OBS/point_stat_-02-46-
39L_197
> > 00101_000000V_cnt.txt
> >
> > // SIGNED //
> > Matthew C. Sittel
> > University Corporation for Atmospheric Research
> > 16 WS/WXN, 557 WW, Offutt AFB, NE
> > matthew.sittel.ctr at us.af.mil
> >
> >
> >
> >
>
>
>
>
------------------------------------------------
More information about the Met_help
mailing list