[Met_help] [rt.rap.ucar.edu #75633] History for More Point Stat Questions
Julie Prestopnik via RT
met_help at ucar.edu
Thu Mar 24 13:45:19 MDT 2016
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
John or Julie, the point stat runs you were helping me with are working for our GALWEM (UM) model. Now I need to do the same thing for WRF (grib 1). In this case, I am verifying Visibility but it is called something different in the grib header (see the last entry below). I am comparing to the same obs file has we discussed before.
rec 57:6417046:date 2016030100 HEATX kpds5=224 kpds6=119 kpds7=9985 levels=(39,1) grid=255 0.99850 (ETA level) 1440min fcst:
HEATX= Heat index [K]
timerange 10 P1 5 P2 160 TimeU 0 nx 620 ny 210 GDS grid 0 num_in_ave 0 missing 0
center 57 subcenter 1 process 220 Table 133 scan: WE:SN winds(N/S)
latlon: lat -0.857000 to 34.843000 by 0.170000 nxny 130200
long -31.064000 to 80.480000 by 0.180000, (620 x 210) scan 64 mode 128 bdsgrid 1
min/max data 941.5 19999.9 num bits 16 BDS_Ref 9415 DecScale 1 BinScale 2
I have tried several different things. First I tried the following:
Fcst data
GRIB1_rec = 57;
GRIB1_lvl_typ=119;
GRIB1_lvl_typ=9985;
level = [ ];
cat_thresh =[<1609,<3218,<4828];
Obs data
name = "SOTYP"; since this is grib1 kpds5 name
level = ["Z0"];
cat_thresh =[<1609,<3218,<4828];
When I run this I get the following error:
['/h/WXQC/met-5.1/bin/point_stat', '/home/stronghs/dust_study/GCWRF/0301/GCWRF_20160301_CY12_FH27.GR1', '/h/data/global/WXQC/data/met/nc_obs/nc_2016030215', '/h/WXQC/met-5.1/data/config/PointStatConfig_dust_gcwrf', '-outdir', '/h/data/global/WXQC/data/met/data/']
DEBUG 1: Default Config File: /h/WXQC/met-5.1/share/met/config/PointStatConfig_default
DEBUG 1: User Config File: /h/WXQC/met-5.1/data/config/PointStatConfig_dust_gcwrf
ERROR :
ERROR : PointStatConfInfo::process_config() -> The number of verification tasks in "obs.field" must be non-zero and match the number in "fcst.field".
ERROR :
I also tried the following:
fcst = {
message_type = [ "" ];
field = [
{
name = "SOTYP";
//GRIB1_rec = 57;
//GRIB1_lvl_typ = 119;
//GRIB1_lvl_val1 = 9985;
level = [ "L9985" ]; or Z0, or Z2 or Z9985
cat_thresh = [<1609,<3218,<4828];
}
];
}
//obs=fcst
obs = {
message_type = [ "ADPSFC_MIN_030000" ];
sid_exc = [];
field = [
{
name = "SOTYP";
level = ["Z0"];
cat_thresh = [<1609,<3218,<4828];
}
I get the following error:
1: Forecast File: /home/stronghs/dust_study/GCWRF/0301/GCWRF_20160301_CY00_FH24.GR1
DEBUG 1: Observation File: /h/data/global/WXQC/data/met/nc_obs/nc_2016030200
DEBUG 2:
DEBUG 2: --------------------------------------------------------------------------------
DEBUG 2:
DEBUG 2: Reading data for SOTYP/L9985.
WARNING:
WARNING: process_fcst_climo_files() -> no fields matching SOTYP/L9985 found in file: /home/stronghs/dust_study/GCWRF/0301/GCWRF_20160301_CY00_FH24.GR1
WARNING:
DEBUG 2:
DEBUG 2: --------------------------------------------------------------------------------
DEBUG 2:
So, I don't know what else to try. I will push the config file and model file to your ftp site.
Thanks
Bob
----------------------------------------------------------------
Complete Ticket History
----------------------------------------------------------------
Subject: More Point Stat Questions
From: Julie Prestopnik
Time: Thu Mar 24 11:20:05 2016
Hi Bob. I'm sorry that you are having trouble with running
point_stat. I
made some changes to your PointStatConfig_dust_gcwrf configuration
file and
was able to get a successful run.
First, in the fcst seciton, I have:
fcst = {
message_type = [ "" ];
field = [
{
GRIB1_ptv = 133;
GRIB1_rec = 224;
level = [ "L9985" ];
cat_thresh = [<1609,<3218,<4828];
}
];
}
Then, in the obs section, I have:
obs = {
message_type = [ "ADPSFC_MIN_030000" ];
sid_exc = [];
field = [
{
GRIB1_ptv = 133;
GRIB1_rec = 211;
level = ["Z2"];
cat_thresh = [<1609,<3218,<4828];
}
];
}
Using those settings, it was able to find the fields. I had to use
"level = ["Z2"];" in the obs field instead of "Z0". In the
obs_arr in
the obs file, the height was listed as "2".
The fcst file we had was GCWRF_20160301_CY00_FH24.GR1 for 20160301,
but the
obs file we had was nc_2016030800 for 20160308, so I increased the
obs_window:
obs_window = {
beg = -7200000;
end = 7200000;
}
in order to get some results. Also, I didn't have your DST.poly file
so I
used the following instead:
mask = {
grid = [ "FULL" ];
poly = [];
sid = [];
}
Please go ahead and make the necessary changes and try it again and
let us
know how it goes.
Thanks,
Julie
On Wed, Mar 23, 2016 at 1:35 PM, robert.craig.2 at us.af.mil via RT <
met_help at ucar.edu> wrote:
>
> Wed Mar 23 13:35:01 2016: Request 75633 was acted upon.
> Transaction: Ticket created by robert.craig.2 at us.af.mil
> Queue: met_help
> Subject: More Point Stat Questions
> Owner: Nobody
> Requestors: robert.craig.2 at us.af.mil
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=75633 >
>
>
> John or Julie, the point stat runs you were helping me with are
working
> for our GALWEM (UM) model. Now I need to do the same thing for WRF
(grib
> 1). In this case, I am verifying Visibility but it is called
something
> different in the grib header (see the last entry below). I am
comparing to
> the same obs file has we discussed before.
>
> rec 57:6417046:date 2016030100 HEATX kpds5=224 kpds6=119 kpds7=9985
> levels=(39,1) grid=255 0.99850 (ETA level) 1440min fcst:
> HEATX= Heat index [K]
> timerange 10 P1 5 P2 160 TimeU 0 nx 620 ny 210 GDS grid 0
num_in_ave 0
> missing 0
> center 57 subcenter 1 process 220 Table 133 scan: WE:SN winds(N/S)
> latlon: lat -0.857000 to 34.843000 by 0.170000 nxny 130200
> long -31.064000 to 80.480000 by 0.180000, (620 x 210) scan
64
> mode 128 bdsgrid 1
> min/max data 941.5 19999.9 num bits 16 BDS_Ref 9415 DecScale 1
> BinScale 2
>
> I have tried several different things. First I tried the following:
>
> Fcst data
>
> GRIB1_rec = 57;
> GRIB1_lvl_typ=119;
> GRIB1_lvl_typ=9985;
> level = [ ];
> cat_thresh =[<1609,<3218,<4828];
>
> Obs data
>
> name = "SOTYP"; since this is grib1 kpds5 name
> level = ["Z0"];
> cat_thresh =[<1609,<3218,<4828];
>
> When I run this I get the following error:
>
> ['/h/WXQC/met-5.1/bin/point_stat',
> '/home/stronghs/dust_study/GCWRF/0301/GCWRF_20160301_CY12_FH27.GR1',
> '/h/data/global/WXQC/data/met/nc_obs/nc_2016030215',
> '/h/WXQC/met-5.1/data/config/PointStatConfig_dust_gcwrf', '-outdir',
> '/h/data/global/WXQC/data/met/data/']
> DEBUG 1: Default Config File:
> /h/WXQC/met-5.1/share/met/config/PointStatConfig_default
> DEBUG 1: User Config File:
> /h/WXQC/met-5.1/data/config/PointStatConfig_dust_gcwrf
> ERROR :
> ERROR : PointStatConfInfo::process_config() -> The number of
verification
> tasks in "obs.field" must be non-zero and match the number in
"fcst.field".
> ERROR :
>
> I also tried the following:
>
> fcst = {
> message_type = [ "" ];
> field = [
> {
> name = "SOTYP";
> //GRIB1_rec = 57;
> //GRIB1_lvl_typ = 119;
> //GRIB1_lvl_val1 = 9985;
> level = [ "L9985" ]; or Z0, or Z2 or Z9985
> cat_thresh = [<1609,<3218,<4828];
> }
>
>
> ];
>
> }
>
> //obs=fcst
>
> obs = {
> message_type = [ "ADPSFC_MIN_030000" ];
> sid_exc = [];
>
> field = [
> {
> name = "SOTYP";
> level = ["Z0"];
> cat_thresh = [<1609,<3218,<4828];
> }
>
>
>
> I get the following error:
>
> 1: Forecast File:
> /home/stronghs/dust_study/GCWRF/0301/GCWRF_20160301_CY00_FH24.GR1
> DEBUG 1: Observation File: /h/data/global/WXQC/data/met/nc_obs/nc_
> 2016030200
> DEBUG 2:
> DEBUG 2:
>
--------------------------------------------------------------------------------
> DEBUG 2:
> DEBUG 2: Reading data for SOTYP/L9985.
> WARNING:
> WARNING: process_fcst_climo_files() -> no fields matching
SOTYP/L9985
> found in file:
> /home/stronghs/dust_study/GCWRF/0301/GCWRF_20160301_CY00_FH24.GR1
> WARNING:
> DEBUG 2:
> DEBUG 2:
>
--------------------------------------------------------------------------------
> DEBUG 2:
>
> So, I don't know what else to try. I will push the config file and
model
> file to your ftp site.
>
> Thanks
> Bob
>
>
>
>
--
Julie Prestopnik
National Center for Atmospheric Research
Research Applications Laboratory
Phone: 303.497.8399
Email: jpresto at ucar.edu
------------------------------------------------
Subject: RE: [rt.rap.ucar.edu #75633] More Point Stat Questions
From: robert.craig.2 at us.af.mil
Time: Thu Mar 24 11:39:17 2016
Thanks Julie it worked. I see you guys got nailed with snow
yesterday, 2ft?
Are you picking table 133 since it has no defined parameters for the
record numbers being used? Also, GRID1_rec doesn’t refer to the
record number but the KPDS 5 value? I am trying to understand the
logic for next time. Am also trying to get the people here to fill
out the grib header correctly so MET doesn't have to be kludged to
read this non-standard parameters.
Again, thanks
Bob
-----Original Message-----
From: Julie Prestopnik via RT [mailto:met_help at ucar.edu]
Sent: Thursday, March 24, 2016 12:20 PM
To: CRAIG, ROBERT J GS-12 USAF ACC 16 WS/WXN
Subject: Re: [rt.rap.ucar.edu #75633] More Point Stat Questions
Hi Bob. I'm sorry that you are having trouble with running
point_stat. I made some changes to your PointStatConfig_dust_gcwrf
configuration file and was able to get a successful run.
First, in the fcst seciton, I have:
fcst = {
message_type = [ "" ];
field = [
{
GRIB1_ptv = 133;
GRIB1_rec = 224;
level = [ "L9985" ];
cat_thresh = [<1609,<3218,<4828];
}
];
}
Then, in the obs section, I have:
obs = {
message_type = [ "ADPSFC_MIN_030000" ];
sid_exc = [];
field = [
{
GRIB1_ptv = 133;
GRIB1_rec = 211;
level = ["Z2"];
cat_thresh = [<1609,<3218,<4828];
}
];
}
Using those settings, it was able to find the fields. I had to use
"level = ["Z2"];" in the obs field instead of "Z0". In the
obs_arr in
the obs file, the height was listed as "2".
The fcst file we had was GCWRF_20160301_CY00_FH24.GR1 for 20160301,
but the obs file we had was nc_2016030800 for 20160308, so I increased
the
obs_window:
obs_window = {
beg = -7200000;
end = 7200000;
}
in order to get some results. Also, I didn't have your DST.poly file
so I used the following instead:
mask = {
grid = [ "FULL" ];
poly = [];
sid = [];
}
Please go ahead and make the necessary changes and try it again and
let us know how it goes.
Thanks,
Julie
On Wed, Mar 23, 2016 at 1:35 PM, robert.craig.2 at us.af.mil via RT <
met_help at ucar.edu> wrote:
>
> Wed Mar 23 13:35:01 2016: Request 75633 was acted upon.
> Transaction: Ticket created by robert.craig.2 at us.af.mil
> Queue: met_help
> Subject: More Point Stat Questions
> Owner: Nobody
> Requestors: robert.craig.2 at us.af.mil
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=75633
> >
>
>
> John or Julie, the point stat runs you were helping me with are
> working for our GALWEM (UM) model. Now I need to do the same thing
for WRF (grib
> 1). In this case, I am verifying Visibility but it is called
something
> different in the grib header (see the last entry below). I am
> comparing to the same obs file has we discussed before.
>
> rec 57:6417046:date 2016030100 HEATX kpds5=224 kpds6=119 kpds7=9985
> levels=(39,1) grid=255 0.99850 (ETA level) 1440min fcst:
> HEATX= Heat index [K]
> timerange 10 P1 5 P2 160 TimeU 0 nx 620 ny 210 GDS grid 0
> num_in_ave 0 missing 0
> center 57 subcenter 1 process 220 Table 133 scan: WE:SN winds(N/S)
> latlon: lat -0.857000 to 34.843000 by 0.170000 nxny 130200
> long -31.064000 to 80.480000 by 0.180000, (620 x 210) scan
> 64 mode 128 bdsgrid 1
> min/max data 941.5 19999.9 num bits 16 BDS_Ref 9415 DecScale 1
> BinScale 2
>
> I have tried several different things. First I tried the following:
>
> Fcst data
>
> GRIB1_rec = 57;
> GRIB1_lvl_typ=119;
> GRIB1_lvl_typ=9985;
> level = [ ];
> cat_thresh =[<1609,<3218,<4828];
>
> Obs data
>
> name = "SOTYP"; since this is grib1 kpds5 name
> level = ["Z0"];
> cat_thresh =[<1609,<3218,<4828];
>
> When I run this I get the following error:
>
> ['/h/WXQC/met-5.1/bin/point_stat',
> '/home/stronghs/dust_study/GCWRF/0301/GCWRF_20160301_CY12_FH27.GR1',
> '/h/data/global/WXQC/data/met/nc_obs/nc_2016030215',
> '/h/WXQC/met-5.1/data/config/PointStatConfig_dust_gcwrf', '-outdir',
> '/h/data/global/WXQC/data/met/data/']
> DEBUG 1: Default Config File:
> /h/WXQC/met-5.1/share/met/config/PointStatConfig_default
> DEBUG 1: User Config File:
> /h/WXQC/met-5.1/data/config/PointStatConfig_dust_gcwrf
> ERROR :
> ERROR : PointStatConfInfo::process_config() -> The number of
> verification tasks in "obs.field" must be non-zero and match the
number in "fcst.field".
> ERROR :
>
> I also tried the following:
>
> fcst = {
> message_type = [ "" ];
> field = [
> {
> name = "SOTYP";
> //GRIB1_rec = 57;
> //GRIB1_lvl_typ = 119;
> //GRIB1_lvl_val1 = 9985;
> level = [ "L9985" ]; or Z0, or Z2 or Z9985
> cat_thresh = [<1609,<3218,<4828];
> }
>
>
> ];
>
> }
>
> //obs=fcst
>
> obs = {
> message_type = [ "ADPSFC_MIN_030000" ];
> sid_exc = [];
>
> field = [
> {
> name = "SOTYP";
> level = ["Z0"];
> cat_thresh = [<1609,<3218,<4828];
> }
>
>
>
> I get the following error:
>
> 1: Forecast File:
> /home/stronghs/dust_study/GCWRF/0301/GCWRF_20160301_CY00_FH24.GR1
> DEBUG 1: Observation File: /h/data/global/WXQC/data/met/nc_obs/nc_
> 2016030200
> DEBUG 2:
> DEBUG 2:
>
----------------------------------------------------------------------
> ----------
> DEBUG 2:
> DEBUG 2: Reading data for SOTYP/L9985.
> WARNING:
> WARNING: process_fcst_climo_files() -> no fields matching
SOTYP/L9985
> found in file:
> /home/stronghs/dust_study/GCWRF/0301/GCWRF_20160301_CY00_FH24.GR1
> WARNING:
> DEBUG 2:
> DEBUG 2:
>
----------------------------------------------------------------------
> ----------
> DEBUG 2:
>
> So, I don't know what else to try. I will push the config file and
> model file to your ftp site.
>
> Thanks
> Bob
>
>
>
>
--
Julie Prestopnik
National Center for Atmospheric Research Research Applications
Laboratory
Phone: 303.497.8399
Email: jpresto at ucar.edu
------------------------------------------------
Subject: More Point Stat Questions
From: Julie Prestopnik
Time: Thu Mar 24 13:45:00 2016
Hi Bob. I'm so glad it worked! We didn't quite get 2 ft, more like
1.25
ft, but still, it was quite a bit, especially for being close to 70
degrees
the day before. :)
To get 133, I first ran:
*wgrib GCWRF_20160301_CY00_FH24.GR1 | grep SOTYP*
Undefined parameter table (center 57-1 table 133), using NCEP-opn
57:6417046:d=16030100:SOTYP:kpds5=224:kpds6=119:kpds7=9985:TR=10:P1=5:P2=160:TimeU=0:0.99850
(ETA level):1440min fcst:NAve=0
Once I had the record number (57), I then ran:
*wgrib -d 57 -V GCWRF_20160301_CY00_FH24.GR1*
Undefined parameter table (center 57-1 table 133), using NCEP-opn
rec 57:6417046:date 2016030100 SOTYP kpds5=224 kpds6=119 kpds7=9985
levels=(39,1) grid=255 0.99850 (ETA level) 1440min fcst:
SOTYP=Soil type (Zobler) [0..9]
timerange 10 P1 5 P2 160 TimeU 0 nx 620 ny 210 GDS grid 0
num_in_ave 0
missing 0
center 57 subcenter 1 process 220 Table 133 scan: WE:SN winds(N/S)
latlon: lat -0.857000 to 34.843000 by 0.170000 nxny 130200
long -31.064000 to 80.480000 by 0.180000, (620 x 210) scan
64
mode 128 bdsgrid 1
min/max data 941.5 19999.9 num bits 16 BDS_Ref 9415 DecScale 1
BinScale 2
In that output, I got "Table 133".
You are correct that GRIB1_rec doesn't refer to the record number, but
rather to the grib code / kpds5 value. We need to update our
documentation
for that. Good catch.
Please let us know if you have any other questions or problems. I'll
go
ahead an close this ticket for now.
Thanks,
Julie
On Thu, Mar 24, 2016 at 11:39 AM, robert.craig.2 at us.af.mil via RT <
met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=75633 >
>
> Thanks Julie it worked. I see you guys got nailed with snow
yesterday,
> 2ft?
>
> Are you picking table 133 since it has no defined parameters for the
> record numbers being used? Also, GRID1_rec doesn’t refer to the
record
> number but the KPDS 5 value? I am trying to understand the logic
for next
> time. Am also trying to get the people here to fill out the grib
header
> correctly so MET doesn't have to be kludged to read this non-
standard
> parameters.
>
> Again, thanks
>
> Bob
>
> -----Original Message-----
> From: Julie Prestopnik via RT [mailto:met_help at ucar.edu]
> Sent: Thursday, March 24, 2016 12:20 PM
> To: CRAIG, ROBERT J GS-12 USAF ACC 16 WS/WXN
> Subject: Re: [rt.rap.ucar.edu #75633] More Point Stat Questions
>
> Hi Bob. I'm sorry that you are having trouble with running
point_stat. I
> made some changes to your PointStatConfig_dust_gcwrf configuration
file and
> was able to get a successful run.
>
> First, in the fcst seciton, I have:
>
> fcst = {
> message_type = [ "" ];
> field = [
> {
> GRIB1_ptv = 133;
> GRIB1_rec = 224;
> level = [ "L9985" ];
> cat_thresh = [<1609,<3218,<4828];
> }
> ];
>
> }
>
>
> Then, in the obs section, I have:
>
> obs = {
> message_type = [ "ADPSFC_MIN_030000" ];
> sid_exc = [];
>
> field = [
> {
> GRIB1_ptv = 133;
> GRIB1_rec = 211;
> level = ["Z2"];
> cat_thresh = [<1609,<3218,<4828];
> }
> ];
>
> }
>
> Using those settings, it was able to find the fields. I had to use
> "level = ["Z2"];" in the obs field instead of "Z0". In the
obs_arr in
> the obs file, the height was listed as "2".
>
> The fcst file we had was GCWRF_20160301_CY00_FH24.GR1 for 20160301,
but
> the obs file we had was nc_2016030800 for 20160308, so I increased
the
> obs_window:
>
> obs_window = {
> beg = -7200000;
> end = 7200000;
> }
>
> in order to get some results. Also, I didn't have your DST.poly
file so I
> used the following instead:
>
> mask = {
> grid = [ "FULL" ];
> poly = [];
> sid = [];
> }
>
>
> Please go ahead and make the necessary changes and try it again and
let us
> know how it goes.
>
> Thanks,
> Julie
>
>
> On Wed, Mar 23, 2016 at 1:35 PM, robert.craig.2 at us.af.mil via RT <
> met_help at ucar.edu> wrote:
>
> >
> > Wed Mar 23 13:35:01 2016: Request 75633 was acted upon.
> > Transaction: Ticket created by robert.craig.2 at us.af.mil
> > Queue: met_help
> > Subject: More Point Stat Questions
> > Owner: Nobody
> > Requestors: robert.craig.2 at us.af.mil
> > Status: new
> > Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=75633
> > >
> >
> >
> > John or Julie, the point stat runs you were helping me with are
> > working for our GALWEM (UM) model. Now I need to do the same
thing for
> WRF (grib
> > 1). In this case, I am verifying Visibility but it is called
something
> > different in the grib header (see the last entry below). I am
> > comparing to the same obs file has we discussed before.
> >
> > rec 57:6417046:date 2016030100 HEATX kpds5=224 kpds6=119
kpds7=9985
> > levels=(39,1) grid=255 0.99850 (ETA level) 1440min fcst:
> > HEATX= Heat index [K]
> > timerange 10 P1 5 P2 160 TimeU 0 nx 620 ny 210 GDS grid 0
> > num_in_ave 0 missing 0
> > center 57 subcenter 1 process 220 Table 133 scan: WE:SN
winds(N/S)
> > latlon: lat -0.857000 to 34.843000 by 0.170000 nxny 130200
> > long -31.064000 to 80.480000 by 0.180000, (620 x 210)
scan
> > 64 mode 128 bdsgrid 1
> > min/max data 941.5 19999.9 num bits 16 BDS_Ref 9415 DecScale
1
> > BinScale 2
> >
> > I have tried several different things. First I tried the
following:
> >
> > Fcst data
> >
> > GRIB1_rec = 57;
> > GRIB1_lvl_typ=119;
> > GRIB1_lvl_typ=9985;
> > level = [ ];
> > cat_thresh =[<1609,<3218,<4828];
> >
> > Obs data
> >
> > name = "SOTYP"; since this is grib1 kpds5 name
> > level = ["Z0"];
> > cat_thresh =[<1609,<3218,<4828];
> >
> > When I run this I get the following error:
> >
> > ['/h/WXQC/met-5.1/bin/point_stat',
> >
'/home/stronghs/dust_study/GCWRF/0301/GCWRF_20160301_CY12_FH27.GR1',
> > '/h/data/global/WXQC/data/met/nc_obs/nc_2016030215',
> > '/h/WXQC/met-5.1/data/config/PointStatConfig_dust_gcwrf', '-
outdir',
> > '/h/data/global/WXQC/data/met/data/']
> > DEBUG 1: Default Config File:
> > /h/WXQC/met-5.1/share/met/config/PointStatConfig_default
> > DEBUG 1: User Config File:
> > /h/WXQC/met-5.1/data/config/PointStatConfig_dust_gcwrf
> > ERROR :
> > ERROR : PointStatConfInfo::process_config() -> The number of
> > verification tasks in "obs.field" must be non-zero and match the
number
> in "fcst.field".
> > ERROR :
> >
> > I also tried the following:
> >
> > fcst = {
> > message_type = [ "" ];
> > field = [
> > {
> > name = "SOTYP";
> > //GRIB1_rec = 57;
> > //GRIB1_lvl_typ = 119;
> > //GRIB1_lvl_val1 = 9985;
> > level = [ "L9985" ]; or Z0, or Z2 or Z9985
> > cat_thresh = [<1609,<3218,<4828];
> > }
> >
> >
> > ];
> >
> > }
> >
> > //obs=fcst
> >
> > obs = {
> > message_type = [ "ADPSFC_MIN_030000" ];
> > sid_exc = [];
> >
> > field = [
> > {
> > name = "SOTYP";
> > level = ["Z0"];
> > cat_thresh = [<1609,<3218,<4828];
> > }
> >
> >
> >
> > I get the following error:
> >
> > 1: Forecast File:
> > /home/stronghs/dust_study/GCWRF/0301/GCWRF_20160301_CY00_FH24.GR1
> > DEBUG 1: Observation File: /h/data/global/WXQC/data/met/nc_obs/nc_
> > 2016030200
> > DEBUG 2:
> > DEBUG 2:
> >
----------------------------------------------------------------------
> > ----------
> > DEBUG 2:
> > DEBUG 2: Reading data for SOTYP/L9985.
> > WARNING:
> > WARNING: process_fcst_climo_files() -> no fields matching
SOTYP/L9985
> > found in file:
> > /home/stronghs/dust_study/GCWRF/0301/GCWRF_20160301_CY00_FH24.GR1
> > WARNING:
> > DEBUG 2:
> > DEBUG 2:
> >
----------------------------------------------------------------------
> > ----------
> > DEBUG 2:
> >
> > So, I don't know what else to try. I will push the config file
and
> > model file to your ftp site.
> >
> > Thanks
> > Bob
> >
> >
> >
> >
>
>
> --
> Julie Prestopnik
> National Center for Atmospheric Research Research Applications
Laboratory
> Phone: 303.497.8399
> Email: jpresto at ucar.edu
>
>
>
>
--
Julie Prestopnik
National Center for Atmospheric Research
Research Applications Laboratory
Phone: 303.497.8399
Email: jpresto at ucar.edu
------------------------------------------------
More information about the Met_help
mailing list