[Met_help] [rt.rap.ucar.edu #58094] History for grib_ptv
Paul Oldenburg via RT
met_help at ucar.edu
Thu Sep 6 13:23:42 MDT 2012
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
I update to METv4 + patches 20120801. I am trying to use a different table
version 129 in order to use the variable 180 L1
in the file METv4.0/data/config/README indicates the use of GRIB1_ptv and
GRIB1_rec .
in the config file is the following setup:
...
{
name = "180";
level = [ "L1" ];
cat_thresh = [ >=40 , >= 60, >= 110 , >= 120 ];
GRIB1_ptv = [129];
}
];
but I get the following
./jal_point_stat.sh
*** Running Point-Stat on Jalisco WRF data ***
DEBUG 1: Default Config File:
/Users/pronos/Documents/WRF/METv4.0/data/config/PointStatConfig_default
DEBUG 1: User Config File: config/PointStatConfig_O3
ERROR :
ERROR : VarInfoGrib::set_dict() - unrecognized GRIB1 field abbreviation
'180' for table version 2
ERROR :
when using the following
{
GRIB1_rec = [180];
GRIB1_ptv = [129];
level = [ "L1" ];
cat_thresh = [ >=40 , >= 60, >= 110 , >= 120 ];
}
I got
*** Running Point-Stat on Jalisco WRF data ***
DEBUG 1: Default Config File:
/Users/pronos/Documents/WRF/METv4.0/data/config/PointStatConfig_default
DEBUG 1: User Config File: config/PointStatConfig_O3
ERROR :
ERROR : VarInfoGrib::set_dict() - either name or GRIB1_ptv and GRIB1_rec
must be specified in field information
ERROR :
How can I setup those variables?
--
Dr. Agustin Garcia-Reynoso
Atmosfera <http://www.atmosfera.unam.mx/fqa>
*Centro de Ciencias de la Atmosfera, *UNAM
Tel (52) 55 *56224073* Fax (52) 55 56160789
Circuito de la Investigacion Cientifica s/n Mexico DF 04510
----------------------------------------------------------------
Complete Ticket History
----------------------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #58094] grib_ptv
From: Paul Oldenburg
Time: Tue Sep 04 08:36:23 2012
Agustin,
I looked up the field you mentioned, table 129 variable 180. In the
standard NCEP GRIB code tables, which is what MET
uses by default, this field corresponds to Ozone concentration in PPB
abbreviated as "OZCON". See the web site:
http://www.nco.ncep.noaa.gov/pmb/docs/on388/table2.html
The information on this website is stored in MET in the table
$MET_BASE/data/table_files/nceptab_flat.txt. If you plan
to use different GRIB code tables, you should copy and rename this
file and modify the original version to suit your
purposes. Then, MET will look for records in your GRIB files based on
the GRIB code abbreviations you specify in your
version of the file.
If you are interested in verifying ozone concentration, then use the
following settings in the point_stat configuration
file:
{
name = "OZCON";
GRIB1_ptv = 129;
level = [ "L1" ];
cat_thresh = [ >=40 , >=60, >=110 , >=120 ];
}
Please let me know if you have any questions or further troubles.
Thanks,
Paul
On 09/01/2012 08:21 AM, Jose Agustín García Reynoso via RT wrote:
>
> Sat Sep 01 08:21:20 2012: Request 58094 was acted upon.
> Transaction: Ticket created by agustin at atmosfera.unam.mx
> Queue: met_help
> Subject: grib_ptv
> Owner: Nobody
> Requestors: agustin at atmosfera.unam.mx
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=58094 >
>
>
> I update to METv4 + patches 20120801. I am trying to use a different
table
> version 129 in order to use the variable 180 L1
> in the file METv4.0/data/config/README indicates the use of
GRIB1_ptv and
> GRIB1_rec .
>
> in the config file is the following setup:
> ...
> {
> name = "180";
> level = [ "L1" ];
> cat_thresh = [ >=40 , >= 60, >= 110 , >= 120 ];
> GRIB1_ptv = [129];
> }
> ];
>
> but I get the following
>
> ./jal_point_stat.sh
>
> *** Running Point-Stat on Jalisco WRF data ***
> DEBUG 1: Default Config File:
>
/Users/pronos/Documents/WRF/METv4.0/data/config/PointStatConfig_default
> DEBUG 1: User Config File: config/PointStatConfig_O3
> ERROR :
> ERROR : VarInfoGrib::set_dict() - unrecognized GRIB1 field
abbreviation
> '180' for table version 2
> ERROR :
>
> when using the following
>
> {
> GRIB1_rec = [180];
> GRIB1_ptv = [129];
> level = [ "L1" ];
> cat_thresh = [ >=40 , >= 60, >= 110 , >= 120 ];
> }
>
> I got
>
> *** Running Point-Stat on Jalisco WRF data ***
> DEBUG 1: Default Config File:
>
/Users/pronos/Documents/WRF/METv4.0/data/config/PointStatConfig_default
> DEBUG 1: User Config File: config/PointStatConfig_O3
> ERROR :
> ERROR : VarInfoGrib::set_dict() - either name or GRIB1_ptv and
GRIB1_rec
> must be specified in field information
> ERROR :
>
> How can I setup those variables?
>
> --
> Dr. Agustin Garcia-Reynoso
> Atmosfera <http://www.atmosfera.unam.mx/fqa>
> *Centro de Ciencias de la Atmosfera, *UNAM
> Tel (52) 55 *56224073* Fax (52) 55 56160789
> Circuito de la Investigacion Cientifica s/n Mexico DF 04510
>
------------------------------------------------
Subject: grib_ptv
From: Jose Agustín García Reynoso
Time: Tue Sep 04 15:21:26 2012
Paul,
I did wat you suggest,
I tried with
GRIB1_ptv = 129;
GRIB1_ptv = [129];
and GRIB1_ptv = "129";
but still having the table 2 as reference table.
DEBUG 1: Default Config File:
/Users/pronos/Documents/WRF/METv4.0/data/config/PointStatConfig_default
DEBUG 1: User Config File:
/Users/pronos/Documents/WRF/METv4.0/scripts/config/PointStatConfig_O3
ERROR :
ERROR : VarInfoGrib::set_dict() - unrecognized GRIB1 field
abbreviation
'OZCON' for table version 2
ERROR :
the log File using -v 6 is attached
thanks
2012/9/4 Paul Oldenburg via RT <met_help at ucar.edu>
> Agustin,
>
> I looked up the field you mentioned, table 129 variable 180. In the
> standard NCEP GRIB code tables, which is what MET
> uses by default, this field corresponds to Ozone concentration in
PPB
> abbreviated as "OZCON". See the web site:
>
> http://www.nco.ncep.noaa.gov/pmb/docs/on388/table2.html
>
> The information on this website is stored in MET in the table
> $MET_BASE/data/table_files/nceptab_flat.txt. If you plan
> to use different GRIB code tables, you should copy and rename this
file
> and modify the original version to suit your
> purposes. Then, MET will look for records in your GRIB files based
on the
> GRIB code abbreviations you specify in your
> version of the file.
>
> If you are interested in verifying ozone concentration, then use the
> following settings in the point_stat configuration
> file:
>
> {
> name = "OZCON";
> GRIB1_ptv = 129;
> level = [ "L1" ];
> cat_thresh = [ >=40 , >=60, >=110 , >=120 ];
> }
>
> Please let me know if you have any questions or further troubles.
>
> Thanks,
>
> Paul
>
>
> On 09/01/2012 08:21 AM, Jose Agustín García Reynoso via RT wrote:
> >
> > Sat Sep 01 08:21:20 2012: Request 58094 was acted upon.
> > Transaction: Ticket created by agustin at atmosfera.unam.mx
> > Queue: met_help
> > Subject: grib_ptv
> > Owner: Nobody
> > Requestors: agustin at atmosfera.unam.mx
> > Status: new
> > Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=58094 >
> >
> >
> > I update to METv4 + patches 20120801. I am trying to use a
different
> table
> > version 129 in order to use the variable 180 L1
> > in the file METv4.0/data/config/README indicates the use of
GRIB1_ptv
> and
> > GRIB1_rec .
> >
> > in the config file is the following setup:
> > ...
> > {
> > name = "180";
> > level = [ "L1" ];
> > cat_thresh = [ >=40 , >= 60, >= 110 , >= 120 ];
> > GRIB1_ptv = [129];
> > }
> > ];
> >
> > but I get the following
> >
> > ./jal_point_stat.sh
> >
> > *** Running Point-Stat on Jalisco WRF data ***
> > DEBUG 1: Default Config File:
> >
/Users/pronos/Documents/WRF/METv4.0/data/config/PointStatConfig_default
> > DEBUG 1: User Config File: config/PointStatConfig_O3
> > ERROR :
> > ERROR : VarInfoGrib::set_dict() - unrecognized GRIB1 field
abbreviation
> > '180' for table version 2
> > ERROR :
> >
> > when using the following
> >
> > {
> > GRIB1_rec = [180];
> > GRIB1_ptv = [129];
> > level = [ "L1" ];
> > cat_thresh = [ >=40 , >= 60, >= 110 , >= 120 ];
> > }
> >
> > I got
> >
> > *** Running Point-Stat on Jalisco WRF data ***
> > DEBUG 1: Default Config File:
> >
/Users/pronos/Documents/WRF/METv4.0/data/config/PointStatConfig_default
> > DEBUG 1: User Config File: config/PointStatConfig_O3
> > ERROR :
> > ERROR : VarInfoGrib::set_dict() - either name or GRIB1_ptv and
GRIB1_rec
> > must be specified in field information
> > ERROR :
> >
> > How can I setup those variables?
> >
> > --
> > Dr. Agustin Garcia-Reynoso
> > Atmosfera <http://www.atmosfera.unam.mx/fqa>
> > *Centro de Ciencias de la Atmosfera, *UNAM
> > Tel (52) 55 *56224073* Fax (52) 55 56160789
> > Circuito de la Investigacion Cientifica s/n Mexico DF 04510
> >
>
>
>
--
Dr. Jose Agustin Garcia Reynoso
Atmosfera <http://www.atmosfera.unam.mx/fqa>
*Centro de Ciencias de la Atmosfera, *UNAM
Tel (52) 55 *56224073* Fax (52) 55 56160789
Circuito de la Investigacion Cientifica s/n Mexico DF 04510
------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #58094] grib_ptv
From: Paul Oldenburg
Time: Tue Sep 04 15:24:23 2012
Agustin,
Can you please send me the data file that you are trying to verify?
That might be the easiest way for me to track this
problem. It is possible that there is a bug in this functionality,
since it is relatively new. If your data file is
larger then about 2MB, please upload it to our FTP site using the
instructions here:
http://www.dtcenter.org/met/users/support/met_help.php#ftp
If you have trouble, please let me know.
Thanks,
Paul
On 09/04/2012 03:21 PM, Jose Agustín García Reynoso via RT wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=58094 >
>
> Paul,
>
> I did wat you suggest,
> I tried with
> GRIB1_ptv = 129;
> GRIB1_ptv = [129];
> and GRIB1_ptv = "129";
>
> but still having the table 2 as reference table.
>
>
> DEBUG 1: Default Config File:
>
/Users/pronos/Documents/WRF/METv4.0/data/config/PointStatConfig_default
> DEBUG 1: User Config File:
>
/Users/pronos/Documents/WRF/METv4.0/scripts/config/PointStatConfig_O3
> ERROR :
> ERROR : VarInfoGrib::set_dict() - unrecognized GRIB1 field
abbreviation
> 'OZCON' for table version 2
> ERROR :
>
> the log File using -v 6 is attached
>
> thanks
>
> 2012/9/4 Paul Oldenburg via RT <met_help at ucar.edu>
>
>> Agustin,
>>
>> I looked up the field you mentioned, table 129 variable 180. In
the
>> standard NCEP GRIB code tables, which is what MET
>> uses by default, this field corresponds to Ozone concentration in
PPB
>> abbreviated as "OZCON". See the web site:
>>
>> http://www.nco.ncep.noaa.gov/pmb/docs/on388/table2.html
>>
>> The information on this website is stored in MET in the table
>> $MET_BASE/data/table_files/nceptab_flat.txt. If you plan
>> to use different GRIB code tables, you should copy and rename this
file
>> and modify the original version to suit your
>> purposes. Then, MET will look for records in your GRIB files based
on the
>> GRIB code abbreviations you specify in your
>> version of the file.
>>
>> If you are interested in verifying ozone concentration, then use
the
>> following settings in the point_stat configuration
>> file:
>>
>> {
>> name = "OZCON";
>> GRIB1_ptv = 129;
>> level = [ "L1" ];
>> cat_thresh = [ >=40 , >=60, >=110 , >=120 ];
>> }
>>
>> Please let me know if you have any questions or further troubles.
>>
>> Thanks,
>>
>> Paul
>>
>>
>> On 09/01/2012 08:21 AM, Jose Agustín García Reynoso via RT wrote:
>>>
>>> Sat Sep 01 08:21:20 2012: Request 58094 was acted upon.
>>> Transaction: Ticket created by agustin at atmosfera.unam.mx
>>> Queue: met_help
>>> Subject: grib_ptv
>>> Owner: Nobody
>>> Requestors: agustin at atmosfera.unam.mx
>>> Status: new
>>> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=58094 >
>>>
>>>
>>> I update to METv4 + patches 20120801. I am trying to use a
different
>> table
>>> version 129 in order to use the variable 180 L1
>>> in the file METv4.0/data/config/README indicates the use of
GRIB1_ptv
>> and
>>> GRIB1_rec .
>>>
>>> in the config file is the following setup:
>>> ...
>>> {
>>> name = "180";
>>> level = [ "L1" ];
>>> cat_thresh = [ >=40 , >= 60, >= 110 , >= 120 ];
>>> GRIB1_ptv = [129];
>>> }
>>> ];
>>>
>>> but I get the following
>>>
>>> ./jal_point_stat.sh
>>>
>>> *** Running Point-Stat on Jalisco WRF data ***
>>> DEBUG 1: Default Config File:
>>>
/Users/pronos/Documents/WRF/METv4.0/data/config/PointStatConfig_default
>>> DEBUG 1: User Config File: config/PointStatConfig_O3
>>> ERROR :
>>> ERROR : VarInfoGrib::set_dict() - unrecognized GRIB1 field
abbreviation
>>> '180' for table version 2
>>> ERROR :
>>>
>>> when using the following
>>>
>>> {
>>> GRIB1_rec = [180];
>>> GRIB1_ptv = [129];
>>> level = [ "L1" ];
>>> cat_thresh = [ >=40 , >= 60, >= 110 , >= 120 ];
>>> }
>>>
>>> I got
>>>
>>> *** Running Point-Stat on Jalisco WRF data ***
>>> DEBUG 1: Default Config File:
>>>
/Users/pronos/Documents/WRF/METv4.0/data/config/PointStatConfig_default
>>> DEBUG 1: User Config File: config/PointStatConfig_O3
>>> ERROR :
>>> ERROR : VarInfoGrib::set_dict() - either name or GRIB1_ptv and
GRIB1_rec
>>> must be specified in field information
>>> ERROR :
>>>
>>> How can I setup those variables?
>>>
>>> --
>>> Dr. Agustin Garcia-Reynoso
>>> Atmosfera <http://www.atmosfera.unam.mx/fqa>
>>> *Centro de Ciencias de la Atmosfera, *UNAM
>>> Tel (52) 55 *56224073* Fax (52) 55 56160789
>>> Circuito de la Investigacion Cientifica s/n Mexico DF 04510
>>>
>>
>>
>>
>
>
------------------------------------------------
Subject: grib_ptv
From: Jose Agustín García Reynoso
Time: Tue Sep 04 19:28:37 2012
Paul,
I uploaded the following files:
PointStatConfig_O3_org
WRFPRS_d01.000
WRFPRS_d01.001
rama_2010.txt
The forecast are WRFPRS files
The observation stations are rama_2001.txt in ascii you have to
convert in
nc using
bin/ascii2nc rama_2010.txt out/ascii2nc/simat_20100100.nc -v 1
for the analysis:
bin/point_stat WRFPRS_d01.000 out/ascii2nc/simat_20100100.nc
PointStatConfig_O3_org
-v 1
or
bin/point_stat WRFPRS_d01.001 out/ascii2nc/simat_20100100.nc
PointStatConfig_O3_org
-v 1
Also we want to analyze the following compounds (NO, NO2, CO, SO2)
fcst_field[] = [ "141/L1","142/L1", "148/L1", "232/L1" ];
fcst_thresh[] = [ "ge105.0 ge210.0","ge105.0 ge210.0", "ge550.0
ge1100.0","ge65 ge130" ];
for grib_ptv = 141
Thanks
2012/9/4 Paul Oldenburg via RT <met_help at ucar.edu>
> Agustin,
>
> Can you please send me the data file that you are trying to verify?
That
> might be the easiest way for me to track this
> problem. It is possible that there is a bug in this functionality,
since
> it is relatively new. If your data file is
> larger then about 2MB, please upload it to our FTP site using the
> instructions here:
>
> http://www.dtcenter.org/met/users/support/met_help.php#ftp
>
> If you have trouble, please let me know.
>
> Thanks,
>
> Paul
>
>
> On 09/04/2012 03:21 PM, Jose Agustín García Reynoso via RT wrote:
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=58094 >
> >
> > Paul,
> >
> > I did wat you suggest,
> > I tried with
> > GRIB1_ptv = 129;
> > GRIB1_ptv = [129];
> > and GRIB1_ptv = "129";
> >
> > but still having the table 2 as reference table.
> >
> >
> > DEBUG 1: Default Config File:
> >
/Users/pronos/Documents/WRF/METv4.0/data/config/PointStatConfig_default
> > DEBUG 1: User Config File:
> >
/Users/pronos/Documents/WRF/METv4.0/scripts/config/PointStatConfig_O3
> > ERROR :
> > ERROR : VarInfoGrib::set_dict() - unrecognized GRIB1 field
abbreviation
> > 'OZCON' for table version 2
> > ERROR :
> >
> > the log File using -v 6 is attached
> >
> > thanks
> >
> > 2012/9/4 Paul Oldenburg via RT <met_help at ucar.edu>
> >
> >> Agustin,
> >>
> >> I looked up the field you mentioned, table 129 variable 180. In
the
> >> standard NCEP GRIB code tables, which is what MET
> >> uses by default, this field corresponds to Ozone concentration in
PPB
> >> abbreviated as "OZCON". See the web site:
> >>
> >> http://www.nco.ncep.noaa.gov/pmb/docs/on388/table2.html
> >>
> >> The information on this website is stored in MET in the table
> >> $MET_BASE/data/table_files/nceptab_flat.txt. If you plan
> >> to use different GRIB code tables, you should copy and rename
this file
> >> and modify the original version to suit your
> >> purposes. Then, MET will look for records in your GRIB files
based on
> the
> >> GRIB code abbreviations you specify in your
> >> version of the file.
> >>
> >> If you are interested in verifying ozone concentration, then use
the
> >> following settings in the point_stat configuration
> >> file:
> >>
> >> {
> >> name = "OZCON";
> >> GRIB1_ptv = 129;
> >> level = [ "L1" ];
> >> cat_thresh = [ >=40 , >=60, >=110 , >=120 ];
> >> }
> >>
> >> Please let me know if you have any questions or further troubles.
> >>
> >> Thanks,
> >>
> >> Paul
> >>
> >>
> >> On 09/01/2012 08:21 AM, Jose Agustín García Reynoso via RT wrote:
> >>>
> >>> Sat Sep 01 08:21:20 2012: Request 58094 was acted upon.
> >>> Transaction: Ticket created by agustin at atmosfera.unam.mx
> >>> Queue: met_help
> >>> Subject: grib_ptv
> >>> Owner: Nobody
> >>> Requestors: agustin at atmosfera.unam.mx
> >>> Status: new
> >>> Ticket <URL:
> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=58094 >
> >>>
> >>>
> >>> I update to METv4 + patches 20120801. I am trying to use a
different
> >> table
> >>> version 129 in order to use the variable 180 L1
> >>> in the file METv4.0/data/config/README indicates the use of
GRIB1_ptv
> >> and
> >>> GRIB1_rec .
> >>>
> >>> in the config file is the following setup:
> >>> ...
> >>> {
> >>> name = "180";
> >>> level = [ "L1" ];
> >>> cat_thresh = [ >=40 , >= 60, >= 110 , >= 120 ];
> >>> GRIB1_ptv = [129];
> >>> }
> >>> ];
> >>>
> >>> but I get the following
> >>>
> >>> ./jal_point_stat.sh
> >>>
> >>> *** Running Point-Stat on Jalisco WRF data ***
> >>> DEBUG 1: Default Config File:
> >>>
/Users/pronos/Documents/WRF/METv4.0/data/config/PointStatConfig_default
> >>> DEBUG 1: User Config File: config/PointStatConfig_O3
> >>> ERROR :
> >>> ERROR : VarInfoGrib::set_dict() - unrecognized GRIB1 field
> abbreviation
> >>> '180' for table version 2
> >>> ERROR :
> >>>
> >>> when using the following
> >>>
> >>> {
> >>> GRIB1_rec = [180];
> >>> GRIB1_ptv = [129];
> >>> level = [ "L1" ];
> >>> cat_thresh = [ >=40 , >= 60, >= 110 , >= 120 ];
> >>> }
> >>>
> >>> I got
> >>>
> >>> *** Running Point-Stat on Jalisco WRF data ***
> >>> DEBUG 1: Default Config File:
> >>>
/Users/pronos/Documents/WRF/METv4.0/data/config/PointStatConfig_default
> >>> DEBUG 1: User Config File: config/PointStatConfig_O3
> >>> ERROR :
> >>> ERROR : VarInfoGrib::set_dict() - either name or GRIB1_ptv and
> GRIB1_rec
> >>> must be specified in field information
> >>> ERROR :
> >>>
> >>> How can I setup those variables?
> >>>
> >>> --
> >>> Dr. Agustin Garcia-Reynoso
> >>> Atmosfera <http://www.atmosfera.unam.mx/fqa>
> >>> *Centro de Ciencias de la Atmosfera, *UNAM
> >>> Tel (52) 55 *56224073* Fax (52) 55 56160789
> >>> Circuito de la Investigacion Cientifica s/n Mexico DF 04510
> >>>
> >>
> >>
> >>
> >
> >
>
>
>
--
Dr. Jose Agustin Garcia Reynoso
Atmosfera <http://www.atmosfera.unam.mx/fqa>
*Centro de Ciencias de la Atmosfera, *UNAM
Tel (52) 55 *56224073* Fax (52) 55 56160789
Circuito de la Investigacion Cientifica s/n Mexico DF 04510
------------------------------------------------
Subject: grib_ptv
From: Paul Oldenburg
Time: Wed Sep 05 09:14:31 2012
Agustin,
I tested the data that you sent and was able to verify several matched
pairs for each of the variables that you
mentioned. The GRIB1_ptv functionality did work as expected, it just
needed to be set up so that each variable is
listed with the correct parameter table version number. I found the
correct parameter table version by searching the
nceptab_flat.txt file that is included in the MET source distribution:
bash$ egrep 'OZCON|"NO2?"|"CO"|"SO"'
$MET_BASE/data/table_files/nceptab_flat.txt
180 129 "OZCON" "Ozone concentration" "ppb"
141 141 "NO" "Nitrogen Oxide" "ppbV"
142 141 "NO2" "Nitrogen Dioxide" "ppbV"
148 141 "CO" "Carbon Monoxide" "ppbV"
The first column is the GRIB code and the second is the parameter
table version number. You can see in the attached
configuration file that these values are set for each parameter. I
ran point_stat with the attached configuration file
using the following command:
bash$ $MET_BASE/bin/point_stat \
WRFPRS_d01.000 \
rama_2010.nc \
PointStatConfig_O3_org \
-outdir . -v 3 > ps_test.log 2>&1 &
The output log is also attached. Note that I changed the boot.n_rep
setting from 1000 to 0 for testing. This feature
is not very practical in the situation where there are relatively few
matched pairs, as there are in this case. I also
modified the output_flag to write out only a single stat file, for
brevity.
Please let me know if you have any questions about this information.
Paul
On 09/04/2012 07:28 PM, Jose Agustín García Reynoso via RT wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=58094 >
>
> Paul,
>
> I uploaded the following files:
> PointStatConfig_O3_org
> WRFPRS_d01.000
> WRFPRS_d01.001
> rama_2010.txt
>
> The forecast are WRFPRS files
> The observation stations are rama_2001.txt in ascii you have to
convert in
> nc using
>
> bin/ascii2nc rama_2010.txt out/ascii2nc/simat_20100100.nc -v 1
>
> for the analysis:
>
> bin/point_stat WRFPRS_d01.000 out/ascii2nc/simat_20100100.nc
> PointStatConfig_O3_org
> -v 1
>
> or
>
> bin/point_stat WRFPRS_d01.001 out/ascii2nc/simat_20100100.nc
> PointStatConfig_O3_org
> -v 1
>
> Also we want to analyze the following compounds (NO, NO2, CO, SO2)
>
> fcst_field[] = [ "141/L1","142/L1", "148/L1", "232/L1" ];
> fcst_thresh[] = [ "ge105.0 ge210.0","ge105.0 ge210.0", "ge550.0
> ge1100.0","ge65 ge130" ];
>
> for grib_ptv = 141
>
> Thanks
>
>
> 2012/9/4 Paul Oldenburg via RT <met_help at ucar.edu>
>
>> Agustin,
>>
>> Can you please send me the data file that you are trying to verify?
That
>> might be the easiest way for me to track this
>> problem. It is possible that there is a bug in this functionality,
since
>> it is relatively new. If your data file is
>> larger then about 2MB, please upload it to our FTP site using the
>> instructions here:
>>
>> http://www.dtcenter.org/met/users/support/met_help.php#ftp
>>
>> If you have trouble, please let me know.
>>
>> Thanks,
>>
>> Paul
>>
>>
>> On 09/04/2012 03:21 PM, Jose Agustín García Reynoso via RT wrote:
>>>
>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=58094 >
>>>
>>> Paul,
>>>
>>> I did wat you suggest,
>>> I tried with
>>> GRIB1_ptv = 129;
>>> GRIB1_ptv = [129];
>>> and GRIB1_ptv = "129";
>>>
>>> but still having the table 2 as reference table.
>>>
>>>
>>> DEBUG 1: Default Config File:
>>>
/Users/pronos/Documents/WRF/METv4.0/data/config/PointStatConfig_default
>>> DEBUG 1: User Config File:
>>>
/Users/pronos/Documents/WRF/METv4.0/scripts/config/PointStatConfig_O3
>>> ERROR :
>>> ERROR : VarInfoGrib::set_dict() - unrecognized GRIB1 field
abbreviation
>>> 'OZCON' for table version 2
>>> ERROR :
>>>
>>> the log File using -v 6 is attached
>>>
>>> thanks
>>>
>>> 2012/9/4 Paul Oldenburg via RT <met_help at ucar.edu>
>>>
>>>> Agustin,
>>>>
>>>> I looked up the field you mentioned, table 129 variable 180. In
the
>>>> standard NCEP GRIB code tables, which is what MET
>>>> uses by default, this field corresponds to Ozone concentration in
PPB
>>>> abbreviated as "OZCON". See the web site:
>>>>
>>>> http://www.nco.ncep.noaa.gov/pmb/docs/on388/table2.html
>>>>
>>>> The information on this website is stored in MET in the table
>>>> $MET_BASE/data/table_files/nceptab_flat.txt. If you plan
>>>> to use different GRIB code tables, you should copy and rename
this file
>>>> and modify the original version to suit your
>>>> purposes. Then, MET will look for records in your GRIB files
based on
>> the
>>>> GRIB code abbreviations you specify in your
>>>> version of the file.
>>>>
>>>> If you are interested in verifying ozone concentration, then use
the
>>>> following settings in the point_stat configuration
>>>> file:
>>>>
>>>> {
>>>> name = "OZCON";
>>>> GRIB1_ptv = 129;
>>>> level = [ "L1" ];
>>>> cat_thresh = [ >=40 , >=60, >=110 , >=120 ];
>>>> }
>>>>
>>>> Please let me know if you have any questions or further troubles.
>>>>
>>>> Thanks,
>>>>
>>>> Paul
>>>>
>>>>
>>>> On 09/01/2012 08:21 AM, Jose Agustín García Reynoso via RT wrote:
>>>>>
>>>>> Sat Sep 01 08:21:20 2012: Request 58094 was acted upon.
>>>>> Transaction: Ticket created by agustin at atmosfera.unam.mx
>>>>> Queue: met_help
>>>>> Subject: grib_ptv
>>>>> Owner: Nobody
>>>>> Requestors: agustin at atmosfera.unam.mx
>>>>> Status: new
>>>>> Ticket <URL:
>> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=58094 >
>>>>>
>>>>>
>>>>> I update to METv4 + patches 20120801. I am trying to use a
different
>>>> table
>>>>> version 129 in order to use the variable 180 L1
>>>>> in the file METv4.0/data/config/README indicates the use of
GRIB1_ptv
>>>> and
>>>>> GRIB1_rec .
>>>>>
>>>>> in the config file is the following setup:
>>>>> ...
>>>>> {
>>>>> name = "180";
>>>>> level = [ "L1" ];
>>>>> cat_thresh = [ >=40 , >= 60, >= 110 , >= 120 ];
>>>>> GRIB1_ptv = [129];
>>>>> }
>>>>> ];
>>>>>
>>>>> but I get the following
>>>>>
>>>>> ./jal_point_stat.sh
>>>>>
>>>>> *** Running Point-Stat on Jalisco WRF data ***
>>>>> DEBUG 1: Default Config File:
>>>>>
/Users/pronos/Documents/WRF/METv4.0/data/config/PointStatConfig_default
>>>>> DEBUG 1: User Config File: config/PointStatConfig_O3
>>>>> ERROR :
>>>>> ERROR : VarInfoGrib::set_dict() - unrecognized GRIB1 field
>> abbreviation
>>>>> '180' for table version 2
>>>>> ERROR :
>>>>>
>>>>> when using the following
>>>>>
>>>>> {
>>>>> GRIB1_rec = [180];
>>>>> GRIB1_ptv = [129];
>>>>> level = [ "L1" ];
>>>>> cat_thresh = [ >=40 , >= 60, >= 110 , >= 120 ];
>>>>> }
>>>>>
>>>>> I got
>>>>>
>>>>> *** Running Point-Stat on Jalisco WRF data ***
>>>>> DEBUG 1: Default Config File:
>>>>>
/Users/pronos/Documents/WRF/METv4.0/data/config/PointStatConfig_default
>>>>> DEBUG 1: User Config File: config/PointStatConfig_O3
>>>>> ERROR :
>>>>> ERROR : VarInfoGrib::set_dict() - either name or GRIB1_ptv and
>> GRIB1_rec
>>>>> must be specified in field information
>>>>> ERROR :
>>>>>
>>>>> How can I setup those variables?
>>>>>
>>>>> --
>>>>> Dr. Agustin Garcia-Reynoso
>>>>> Atmosfera <http://www.atmosfera.unam.mx/fqa>
>>>>> *Centro de Ciencias de la Atmosfera, *UNAM
>>>>> Tel (52) 55 *56224073* Fax (52) 55 56160789
>>>>> Circuito de la Investigacion Cientifica s/n Mexico DF 04510
>>>>>
>>>>
>>>>
>>>>
>>>
>>>
>>
>>
>>
>
>
------------------------------------------------
Subject: grib_ptv
From: Paul Oldenburg
Time: Wed Sep 05 09:14:31 2012
////////////////////////////////////////////////////////////////////////////////
//
// Point-Stat configuration file.
//
// For additional information, see the MET_BASE/data/config/README
file.
//
////////////////////////////////////////////////////////////////////////////////
//
// Output model name to be written
//
model = "WRF";
////////////////////////////////////////////////////////////////////////////////
//
// Forecast and observation fields to be verified
//
fcst = {
wind_thresh = [ NA ];
message_type = [ "ADPSFC" ];
field = [
{
name = "TMP";
level = [ "Z2" ];
cat_thresh = [ >283, >288, >293 ];
GRIB1_ptv = 2;
},
{
name = "RH";
level = [ "Z2" ];
cat_thresh = [ >=10 ];
GRIB1_ptv = 2;
},
{
name = "UGRD";
level = [ "Z10" ];
cat_thresh = [ >=5 ];
GRIB1_ptv = 2;
},
{
name = "VGRD";
level = [ "Z10" ];
cat_thresh = [ >=5 ];
GRIB1_ptv = 2;
},
{
name = "OZCON";
level = [ "L1" ];
cat_thresh = [ >=50, >=110, >=120 ];
GRIB1_ptv = 129;
},
{
name = "NO";
level = [ "L1" ];
cat_thresh = [ >=105, >=210 ];
GRIB1_ptv = 141;
},
{
name = "NO2";
level = [ "L1" ];
cat_thresh = [ >=105, >=210 ];
GRIB1_ptv = 141;
},
{
name = "CO";
level = [ "L1" ];
cat_thresh = [ >=550, >=1100 ];
GRIB1_ptv = 141;
},
{
name = "SO2";
level = [ "L1" ];
cat_thresh = [ >=65, >=130 ];
GRIB1_ptv = 141;
}
];
};
obs = fcst;
////////////////////////////////////////////////////////////////////////////////
//
// Point observation time window
//
obs_window = {
beg = -2700;
end = 2700;
}
////////////////////////////////////////////////////////////////////////////////
//
// Verification masking regions
//
mask = {
grid = [ ];
poly = [ "MET_BASE/data/poly/MEX.poly" ];
sid = "";
};
////////////////////////////////////////////////////////////////////////////////
//
// Confidence interval settings
//
ci_alpha = [ 0.05 ];
boot = {
interval = PCTILE;
rep_prop = 1.0;
n_rep = 0;
rng = "mt19937";
seed = "1";
};
////////////////////////////////////////////////////////////////////////////////
//
// Interpolation methods
//
interp = {
vld_thresh = 1.0;
type = [
{
method = UW_MEAN;
width = 2;
}
];
};
////////////////////////////////////////////////////////////////////////////////
//
// Statistical output types
//
output_flag = {
fho = STAT;
ctc = STAT;
cts = STAT;
mctc = NONE;
mcts = NONE;
cnt = STAT;
sl1l2 = STAT;
sal1l2 = NONE;
vl1l2 = STAT;
val1l2 = NONE;
pct = NONE;
pstd = NONE;
pjc = NONE;
prc = NONE;
mpr = NONE;
};
////////////////////////////////////////////////////////////////////////////////
duplicate_flag = NONE;
rank_corr_flag = TRUE;
tmp_dir = "/tmp";
output_prefix = "jal_o3";
version = "V4.0";
////////////////////////////////////////////////////////////////////////////////
------------------------------------------------
Subject: grib_ptv
From: Paul Oldenburg
Time: Wed Sep 05 09:14:31 2012
DEBUG 1: Default Config File:
/d1/pgoldenb/opt/MET_builds/METv4.0/data/config/PointStatConfig_default
DEBUG 1: User Config File: PointStatConfig_O3_org
GSL_RNG_TYPE=mt19937
GSL_RNG_SEED=1
DEBUG 1: Forecast File: WRFPRS_d01.000
DEBUG 1: Climatology File: none
DEBUG 1: Observation File: rama_2010.nc
DEBUG 2:
DEBUG 2:
--------------------------------------------------------------------------------
DEBUG 2:
DEBUG 2: Reading data for TMP/Z2.
DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found range match for
VarInfo "TMP/Z2" in GRIB record 250 of GRIB file "WRFPRS_d01.000".
DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 1 GRIB records
matching VarInfo "TMP/Z2" in GRIB file "WRFPRS_d01.000".
DEBUG 2: For TMP/Z2 found 1 forecast levels and 0 climatology levels.
DEBUG 2:
DEBUG 2:
--------------------------------------------------------------------------------
DEBUG 2:
DEBUG 2: Reading data for RH/Z2.
DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found range match for
VarInfo "RH/Z2" in GRIB record 253 of GRIB file "WRFPRS_d01.000".
DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 1 GRIB records
matching VarInfo "RH/Z2" in GRIB file "WRFPRS_d01.000".
DEBUG 2: For RH/Z2 found 1 forecast levels and 0 climatology levels.
DEBUG 2:
DEBUG 2:
--------------------------------------------------------------------------------
DEBUG 2:
DEBUG 2: Reading data for UGRD/Z10.
DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found range match for
VarInfo "UGRD/Z10" in GRIB record 254 of GRIB file "WRFPRS_d01.000".
DEBUG 3: MetGrib1DataFile::rotate_winds() -> Have U-wind record,
reading V-wind record.
DEBUG 3: MetGrib1DataFile::data_plane_scalar() -> Found exact match
for VarInfo "UGRD/Z10" in GRIB record 255 of GRIB file
"WRFPRS_d01.000".
DEBUG 3: Rotating U and V wind components from grid-relative to earth-
relative.
DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 1 GRIB records
matching VarInfo "UGRD/Z10" in GRIB file "WRFPRS_d01.000".
DEBUG 2: For UGRD/Z10 found 1 forecast levels and 0 climatology
levels.
DEBUG 2:
DEBUG 2:
--------------------------------------------------------------------------------
DEBUG 2:
DEBUG 2: Reading data for VGRD/Z10.
DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found range match for
VarInfo "VGRD/Z10" in GRIB record 255 of GRIB file "WRFPRS_d01.000".
DEBUG 3: MetGrib1DataFile::rotate_winds() -> Have V-wind record,
reading U-wind record.
DEBUG 3: MetGrib1DataFile::data_plane_scalar() -> Found exact match
for VarInfo "VGRD/Z10" in GRIB record 254 of GRIB file
"WRFPRS_d01.000".
DEBUG 3: Rotating U and V wind components from grid-relative to earth-
relative.
DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 1 GRIB records
matching VarInfo "VGRD/Z10" in GRIB file "WRFPRS_d01.000".
DEBUG 2: For VGRD/Z10 found 1 forecast levels and 0 climatology
levels.
DEBUG 2:
DEBUG 2:
--------------------------------------------------------------------------------
DEBUG 2:
DEBUG 2: Reading data for OZCON/L1.
DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found range match for
VarInfo "OZCON/L1" in GRIB record 16 of GRIB file "WRFPRS_d01.000".
DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 1 GRIB records
matching VarInfo "OZCON/L1" in GRIB file "WRFPRS_d01.000".
DEBUG 2: For OZCON/L1 found 1 forecast levels and 0 climatology
levels.
DEBUG 2:
DEBUG 2:
--------------------------------------------------------------------------------
DEBUG 2:
DEBUG 2: Reading data for NO/L1.
DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found range match for
VarInfo "NO/L1" in GRIB record 11 of GRIB file "WRFPRS_d01.000".
DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 1 GRIB records
matching VarInfo "NO/L1" in GRIB file "WRFPRS_d01.000".
DEBUG 2: For NO/L1 found 1 forecast levels and 0 climatology levels.
DEBUG 2:
DEBUG 2:
--------------------------------------------------------------------------------
DEBUG 2:
DEBUG 2: Reading data for NO2/L1.
DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found range match for
VarInfo "NO2/L1" in GRIB record 12 of GRIB file "WRFPRS_d01.000".
DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 1 GRIB records
matching VarInfo "NO2/L1" in GRIB file "WRFPRS_d01.000".
DEBUG 2: For NO2/L1 found 1 forecast levels and 0 climatology levels.
DEBUG 2:
DEBUG 2:
--------------------------------------------------------------------------------
DEBUG 2:
DEBUG 2: Reading data for CO/L1.
DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found range match for
VarInfo "CO/L1" in GRIB record 14 of GRIB file "WRFPRS_d01.000".
DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 1 GRIB records
matching VarInfo "CO/L1" in GRIB file "WRFPRS_d01.000".
DEBUG 2: For CO/L1 found 1 forecast levels and 0 climatology levels.
DEBUG 2:
DEBUG 2:
--------------------------------------------------------------------------------
DEBUG 2:
DEBUG 2: Reading data for SO2/L1.
DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found range match for
VarInfo "SO2/L1" in GRIB record 13 of GRIB file "WRFPRS_d01.000".
DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 1 GRIB records
matching VarInfo "SO2/L1" in GRIB file "WRFPRS_d01.000".
DEBUG 2: For SO2/L1 found 1 forecast levels and 0 climatology levels.
DEBUG 2:
DEBUG 2:
--------------------------------------------------------------------------------
DEBUG 2:
DEBUG 2: Searching 2680560 observations from 297840 messages.
DEBUG 2:
DEBUG 2:
--------------------------------------------------------------------------------
DEBUG 2:
DEBUG 2: Processing TMP/Z2 versus TMP/Z2, for observation type ADPSFC,
over region MEX, for interpolation method UW_MEAN(4), using 14 pairs.
DEBUG 3: Number of matched pairs = 14
DEBUG 3: Observations processed = 2680560
DEBUG 3: Rejected: GRIB code = 2382720
DEBUG 3: Rejected: valid time = 297806
DEBUG 3: Rejected: bad obs value = 20
DEBUG 3: Rejected: off the grid = 0
DEBUG 3: Rejected: level mismatch = 0
DEBUG 3: Rejected: message type = 0
DEBUG 3: Rejected: masking region = 0
DEBUG 3: Rejected: bad fcst value = 0
DEBUG 2: Computing Categorical Statistics.
DEBUG 2: Computing Continuous Statistics.
DEBUG 2: Computing Scalar Partial Sums.
DEBUG 2:
DEBUG 2:
--------------------------------------------------------------------------------
DEBUG 2:
DEBUG 2: Processing RH/Z2 versus RH/Z2, for observation type ADPSFC,
over region MEX, for interpolation method UW_MEAN(4), using 14 pairs.
DEBUG 3: Number of matched pairs = 14
DEBUG 3: Observations processed = 2680560
DEBUG 3: Rejected: GRIB code = 2382720
DEBUG 3: Rejected: valid time = 297806
DEBUG 3: Rejected: bad obs value = 20
DEBUG 3: Rejected: off the grid = 0
DEBUG 3: Rejected: level mismatch = 0
DEBUG 3: Rejected: message type = 0
DEBUG 3: Rejected: masking region = 0
DEBUG 3: Rejected: bad fcst value = 0
DEBUG 2: Computing Categorical Statistics.
DEBUG 2: Computing Continuous Statistics.
DEBUG 2: Computing Scalar Partial Sums.
DEBUG 2:
DEBUG 2:
--------------------------------------------------------------------------------
DEBUG 2:
DEBUG 2: Processing UGRD/Z10 versus UGRD/Z10, for observation type
ADPSFC, over region MEX, for interpolation method UW_MEAN(4), using 10
pairs.
DEBUG 3: Number of matched pairs = 10
DEBUG 3: Observations processed = 2680560
DEBUG 3: Rejected: GRIB code = 2382720
DEBUG 3: Rejected: valid time = 297806
DEBUG 3: Rejected: bad obs value = 24
DEBUG 3: Rejected: off the grid = 0
DEBUG 3: Rejected: level mismatch = 0
DEBUG 3: Rejected: message type = 0
DEBUG 3: Rejected: masking region = 0
DEBUG 3: Rejected: bad fcst value = 0
DEBUG 2: Computing Categorical Statistics.
DEBUG 2: Computing Continuous Statistics.
DEBUG 2:
DEBUG 2:
--------------------------------------------------------------------------------
DEBUG 2:
DEBUG 2: Processing VGRD/Z10 versus VGRD/Z10, for observation type
ADPSFC, over region MEX, for interpolation method UW_MEAN(4), using 10
pairs.
DEBUG 3: Number of matched pairs = 10
DEBUG 3: Observations processed = 2680560
DEBUG 3: Rejected: GRIB code = 2382720
DEBUG 3: Rejected: valid time = 297806
DEBUG 3: Rejected: bad obs value = 24
DEBUG 3: Rejected: off the grid = 0
DEBUG 3: Rejected: level mismatch = 0
DEBUG 3: Rejected: message type = 0
DEBUG 3: Rejected: masking region = 0
DEBUG 3: Rejected: bad fcst value = 0
DEBUG 2: Computing Categorical Statistics.
DEBUG 2: Computing Continuous Statistics.
DEBUG 2: Computing Vector Partial Sums, for forecast wind speed NA,
for observation wind speed NA, using 10 pairs.
DEBUG 2:
DEBUG 2:
--------------------------------------------------------------------------------
DEBUG 2:
DEBUG 2: Processing OZCON/L1 versus OZCON/L1, for observation type
ADPSFC, over region MEX, for interpolation method UW_MEAN(4), using 19
pairs.
DEBUG 3: Number of matched pairs = 19
DEBUG 3: Observations processed = 2680560
DEBUG 3: Rejected: GRIB code = 2382720
DEBUG 3: Rejected: valid time = 297806
DEBUG 3: Rejected: bad obs value = 15
DEBUG 3: Rejected: off the grid = 0
DEBUG 3: Rejected: level mismatch = 0
DEBUG 3: Rejected: message type = 0
DEBUG 3: Rejected: masking region = 0
DEBUG 3: Rejected: bad fcst value = 0
DEBUG 2: Computing Categorical Statistics.
DEBUG 2: Computing Continuous Statistics.
DEBUG 2: Computing Scalar Partial Sums.
DEBUG 2:
DEBUG 2:
--------------------------------------------------------------------------------
DEBUG 2:
DEBUG 2: Processing NO/L1 versus NO/L1, for observation type ADPSFC,
over region MEX, for interpolation method UW_MEAN(4), using 14 pairs.
DEBUG 3: Number of matched pairs = 14
DEBUG 3: Observations processed = 2680560
DEBUG 3: Rejected: GRIB code = 2382720
DEBUG 3: Rejected: valid time = 297806
DEBUG 3: Rejected: bad obs value = 20
DEBUG 3: Rejected: off the grid = 0
DEBUG 3: Rejected: level mismatch = 0
DEBUG 3: Rejected: message type = 0
DEBUG 3: Rejected: masking region = 0
DEBUG 3: Rejected: bad fcst value = 0
DEBUG 2: Computing Categorical Statistics.
DEBUG 2: Computing Continuous Statistics.
DEBUG 2: Computing Scalar Partial Sums.
DEBUG 2:
DEBUG 2:
--------------------------------------------------------------------------------
DEBUG 2:
DEBUG 2: Processing NO2/L1 versus NO2/L1, for observation type ADPSFC,
over region MEX, for interpolation method UW_MEAN(4), using 14 pairs.
DEBUG 3: Number of matched pairs = 14
DEBUG 3: Observations processed = 2680560
DEBUG 3: Rejected: GRIB code = 2382720
DEBUG 3: Rejected: valid time = 297806
DEBUG 3: Rejected: bad obs value = 20
DEBUG 3: Rejected: off the grid = 0
DEBUG 3: Rejected: level mismatch = 0
DEBUG 3: Rejected: message type = 0
DEBUG 3: Rejected: masking region = 0
DEBUG 3: Rejected: bad fcst value = 0
DEBUG 2: Computing Categorical Statistics.
DEBUG 2: Computing Continuous Statistics.
DEBUG 2: Computing Scalar Partial Sums.
DEBUG 2:
DEBUG 2:
--------------------------------------------------------------------------------
DEBUG 2:
DEBUG 2: Processing CO/L1 versus CO/L1, for observation type ADPSFC,
over region MEX, for interpolation method UW_MEAN(4), using 14 pairs.
DEBUG 3: Number of matched pairs = 14
DEBUG 3: Observations processed = 2680560
DEBUG 3: Rejected: GRIB code = 2382720
DEBUG 3: Rejected: valid time = 297806
DEBUG 3: Rejected: bad obs value = 20
DEBUG 3: Rejected: off the grid = 0
DEBUG 3: Rejected: level mismatch = 0
DEBUG 3: Rejected: message type = 0
DEBUG 3: Rejected: masking region = 0
DEBUG 3: Rejected: bad fcst value = 0
DEBUG 2: Computing Categorical Statistics.
DEBUG 2: Computing Continuous Statistics.
DEBUG 2: Computing Scalar Partial Sums.
DEBUG 2:
DEBUG 2:
--------------------------------------------------------------------------------
DEBUG 2:
DEBUG 2: Processing SO2/L1 versus SO2/L1, for observation type ADPSFC,
over region MEX, for interpolation method UW_MEAN(4), using 22 pairs.
DEBUG 3: Number of matched pairs = 22
DEBUG 3: Observations processed = 2680560
DEBUG 3: Rejected: GRIB code = 2382720
DEBUG 3: Rejected: valid time = 297806
DEBUG 3: Rejected: bad obs value = 12
DEBUG 3: Rejected: off the grid = 0
DEBUG 3: Rejected: level mismatch = 0
DEBUG 3: Rejected: message type = 0
DEBUG 3: Rejected: masking region = 0
DEBUG 3: Rejected: bad fcst value = 0
DEBUG 2: Computing Categorical Statistics.
DEBUG 2: Computing Continuous Statistics.
DEBUG 2: Computing Scalar Partial Sums.
DEBUG 2:
DEBUG 2:
--------------------------------------------------------------------------------
DEBUG 2:
DEBUG 1: Output file:
./point_stat_jal_o3_000000L_20100502_120000V.stat
------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #58094] grib_ptv
From: Jose Agustín García Reynoso
Time: Thu Sep 06 13:15:29 2012
Paul,
Now is working as expected .
May be will be very useful to include an example in the README file
in the METv4.0/data/config directory.
Thanks
Jose Agustín García Reynoso http://www.atmosfera.unam.mx/fqa
Centro de Ciencias de la Atmosfera http://www.atmosfera.unam.mx
52-55-56224073
Universidad Nacional Autonoma de Mexico 52-55-56224396 (fax)
Circuito Exterior s/n Mexico DF 04510
El 05/09/2012, a las 10:14, Paul Oldenburg via RT <met_help at ucar.edu>
escribió:
> Agustin,
>
> I tested the data that you sent and was able to verify several
matched pairs for each of the variables that you
> mentioned. The GRIB1_ptv functionality did work as expected, it
just needed to be set up so that each variable is
> listed with the correct parameter table version number. I found the
correct parameter table version by searching the
> nceptab_flat.txt file that is included in the MET source
distribution:
>
> bash$ egrep 'OZCON|"NO2?"|"CO"|"SO"'
$MET_BASE/data/table_files/nceptab_flat.txt
> 180 129 "OZCON" "Ozone concentration" "ppb"
> 141 141 "NO" "Nitrogen Oxide" "ppbV"
> 142 141 "NO2" "Nitrogen Dioxide" "ppbV"
> 148 141 "CO" "Carbon Monoxide" "ppbV"
>
> The first column is the GRIB code and the second is the parameter
table version number. You can see in the attached
> configuration file that these values are set for each parameter. I
ran point_stat with the attached configuration file
> using the following command:
>
> bash$ $MET_BASE/bin/point_stat \
> WRFPRS_d01.000 \
> rama_2010.nc \
> PointStatConfig_O3_org \
> -outdir . -v 3 > ps_test.log 2>&1 &
>
> The output log is also attached. Note that I changed the boot.n_rep
setting from 1000 to 0 for testing. This feature
> is not very practical in the situation where there are relatively
few matched pairs, as there are in this case. I also
> modified the output_flag to write out only a single stat file, for
brevity.
>
> Please let me know if you have any questions about this information.
>
> Paul
>
>
> On 09/04/2012 07:28 PM, Jose Agustín García Reynoso via RT wrote:
>>
>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=58094 >
>>
>> Paul,
>>
>> I uploaded the following files:
>> PointStatConfig_O3_org
>> WRFPRS_d01.000
>> WRFPRS_d01.001
>> rama_2010.txt
>>
>> The forecast are WRFPRS files
>> The observation stations are rama_2001.txt in ascii you have to
convert in
>> nc using
>>
>> bin/ascii2nc rama_2010.txt out/ascii2nc/simat_20100100.nc -v 1
>>
>> for the analysis:
>>
>> bin/point_stat WRFPRS_d01.000 out/ascii2nc/simat_20100100.nc
>> PointStatConfig_O3_org
>> -v 1
>>
>> or
>>
>> bin/point_stat WRFPRS_d01.001 out/ascii2nc/simat_20100100.nc
>> PointStatConfig_O3_org
>> -v 1
>>
>> Also we want to analyze the following compounds (NO, NO2, CO, SO2)
>>
>> fcst_field[] = [ "141/L1","142/L1", "148/L1", "232/L1" ];
>> fcst_thresh[] = [ "ge105.0 ge210.0","ge105.0 ge210.0", "ge550.0
>> ge1100.0","ge65 ge130" ];
>>
>> for grib_ptv = 141
>>
>> Thanks
>>
>>
>> 2012/9/4 Paul Oldenburg via RT <met_help at ucar.edu>
>>
>>> Agustin,
>>>
>>> Can you please send me the data file that you are trying to
verify? That
>>> might be the easiest way for me to track this
>>> problem. It is possible that there is a bug in this
functionality, since
>>> it is relatively new. If your data file is
>>> larger then about 2MB, please upload it to our FTP site using the
>>> instructions here:
>>>
>>> http://www.dtcenter.org/met/users/support/met_help.php#ftp
>>>
>>> If you have trouble, please let me know.
>>>
>>> Thanks,
>>>
>>> Paul
>>>
>>>
>>> On 09/04/2012 03:21 PM, Jose Agustín García Reynoso via RT wrote:
>>>>
>>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=58094 >
>>>>
>>>> Paul,
>>>>
>>>> I did wat you suggest,
>>>> I tried with
>>>> GRIB1_ptv = 129;
>>>> GRIB1_ptv = [129];
>>>> and GRIB1_ptv = "129";
>>>>
>>>> but still having the table 2 as reference table.
>>>>
>>>>
>>>> DEBUG 1: Default Config File:
>>>>
/Users/pronos/Documents/WRF/METv4.0/data/config/PointStatConfig_default
>>>> DEBUG 1: User Config File:
>>>>
/Users/pronos/Documents/WRF/METv4.0/scripts/config/PointStatConfig_O3
>>>> ERROR :
>>>> ERROR : VarInfoGrib::set_dict() - unrecognized GRIB1 field
abbreviation
>>>> 'OZCON' for table version 2
>>>> ERROR :
>>>>
>>>> the log File using -v 6 is attached
>>>>
>>>> thanks
>>>>
>>>> 2012/9/4 Paul Oldenburg via RT <met_help at ucar.edu>
>>>>
>>>>> Agustin,
>>>>>
>>>>> I looked up the field you mentioned, table 129 variable 180. In
the
>>>>> standard NCEP GRIB code tables, which is what MET
>>>>> uses by default, this field corresponds to Ozone concentration
in PPB
>>>>> abbreviated as "OZCON". See the web site:
>>>>>
>>>>> http://www.nco.ncep.noaa.gov/pmb/docs/on388/table2.html
>>>>>
>>>>> The information on this website is stored in MET in the table
>>>>> $MET_BASE/data/table_files/nceptab_flat.txt. If you plan
>>>>> to use different GRIB code tables, you should copy and rename
this file
>>>>> and modify the original version to suit your
>>>>> purposes. Then, MET will look for records in your GRIB files
based on
>>> the
>>>>> GRIB code abbreviations you specify in your
>>>>> version of the file.
>>>>>
>>>>> If you are interested in verifying ozone concentration, then use
the
>>>>> following settings in the point_stat configuration
>>>>> file:
>>>>>
>>>>> {
>>>>> name = "OZCON";
>>>>> GRIB1_ptv = 129;
>>>>> level = [ "L1" ];
>>>>> cat_thresh = [ >=40 , >=60, >=110 , >=120 ];
>>>>> }
>>>>>
>>>>> Please let me know if you have any questions or further
troubles.
>>>>>
>>>>> Thanks,
>>>>>
>>>>> Paul
>>>>>
>>>>>
>>>>> On 09/01/2012 08:21 AM, Jose Agustín García Reynoso via RT
wrote:
>>>>>>
>>>>>> Sat Sep 01 08:21:20 2012: Request 58094 was acted upon.
>>>>>> Transaction: Ticket created by agustin at atmosfera.unam.mx
>>>>>> Queue: met_help
>>>>>> Subject: grib_ptv
>>>>>> Owner: Nobody
>>>>>> Requestors: agustin at atmosfera.unam.mx
>>>>>> Status: new
>>>>>> Ticket <URL:
>>> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=58094 >
>>>>>>
>>>>>>
>>>>>> I update to METv4 + patches 20120801. I am trying to use a
different
>>>>> table
>>>>>> version 129 in order to use the variable 180 L1
>>>>>> in the file METv4.0/data/config/README indicates the use of
GRIB1_ptv
>>>>> and
>>>>>> GRIB1_rec .
>>>>>>
>>>>>> in the config file is the following setup:
>>>>>> ...
>>>>>> {
>>>>>> name = "180";
>>>>>> level = [ "L1" ];
>>>>>> cat_thresh = [ >=40 , >= 60, >= 110 , >= 120 ];
>>>>>> GRIB1_ptv = [129];
>>>>>> }
>>>>>> ];
>>>>>>
>>>>>> but I get the following
>>>>>>
>>>>>> ./jal_point_stat.sh
>>>>>>
>>>>>> *** Running Point-Stat on Jalisco WRF data ***
>>>>>> DEBUG 1: Default Config File:
>>>>>>
/Users/pronos/Documents/WRF/METv4.0/data/config/PointStatConfig_default
>>>>>> DEBUG 1: User Config File: config/PointStatConfig_O3
>>>>>> ERROR :
>>>>>> ERROR : VarInfoGrib::set_dict() - unrecognized GRIB1 field
>>> abbreviation
>>>>>> '180' for table version 2
>>>>>> ERROR :
>>>>>>
>>>>>> when using the following
>>>>>>
>>>>>> {
>>>>>> GRIB1_rec = [180];
>>>>>> GRIB1_ptv = [129];
>>>>>> level = [ "L1" ];
>>>>>> cat_thresh = [ >=40 , >= 60, >= 110 , >= 120 ];
>>>>>> }
>>>>>>
>>>>>> I got
>>>>>>
>>>>>> *** Running Point-Stat on Jalisco WRF data ***
>>>>>> DEBUG 1: Default Config File:
>>>>>>
/Users/pronos/Documents/WRF/METv4.0/data/config/PointStatConfig_default
>>>>>> DEBUG 1: User Config File: config/PointStatConfig_O3
>>>>>> ERROR :
>>>>>> ERROR : VarInfoGrib::set_dict() - either name or GRIB1_ptv and
>>> GRIB1_rec
>>>>>> must be specified in field information
>>>>>> ERROR :
>>>>>>
>>>>>> How can I setup those variables?
>>>>>>
>>>>>> --
>>>>>> Dr. Agustin Garcia-Reynoso
>>>>>> Atmosfera <http://www.atmosfera.unam.mx/fqa>
>>>>>> *Centro de Ciencias de la Atmosfera, *UNAM
>>>>>> Tel (52) 55 *56224073* Fax (52) 55 56160789
>>>>>> Circuito de la Investigacion Cientifica s/n Mexico DF 04510
>>>>>>
>>>>>
>>>>>
>>>>>
>>>>
>>>>
>>>
>>>
>>>
>>
>>
>
>
> <PointStatConfig_O3_org.txt><ps_test.log>
------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #58094] grib_ptv
From: Paul Oldenburg
Time: Thu Sep 06 13:23:19 2012
Agustin,
I'm glad to hear that point_stat is working for you. There are a
couple simple examples of how to use MET tools
included in the $MET_BASE/scripts directory, and there are also more
in-depth online tutorials available on the MET
website, if you are interested:
http://www.dtcenter.org/met/users/support/online_tutorial/index.php
http://www.dtcenter.org/met/users/support/tutorial.php
Thanks for your suggestion.
Paul
On 09/06/2012 01:15 PM, Jose Agustín García Reynoso via RT wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=58094 >
>
> Paul,
>
> Now is working as expected .
> May be will be very useful to include an example in the README
file in the METv4.0/data/config directory.
>
> Thanks
>
> Jose Agustín García Reynoso http://www.atmosfera.unam.mx/fqa
> Centro de Ciencias de la Atmosfera http://www.atmosfera.unam.mx
> 52-55-56224073
> Universidad Nacional Autonoma de Mexico 52-55-56224396 (fax)
> Circuito Exterior s/n Mexico DF 04510
>
>
>
>
>
>
> El 05/09/2012, a las 10:14, Paul Oldenburg via RT
<met_help at ucar.edu> escribió:
>
>> Agustin,
>>
>> I tested the data that you sent and was able to verify several
matched pairs for each of the variables that you
>> mentioned. The GRIB1_ptv functionality did work as expected, it
just needed to be set up so that each variable is
>> listed with the correct parameter table version number. I found
the correct parameter table version by searching the
>> nceptab_flat.txt file that is included in the MET source
distribution:
>>
>> bash$ egrep 'OZCON|"NO2?"|"CO"|"SO"'
$MET_BASE/data/table_files/nceptab_flat.txt
>> 180 129 "OZCON" "Ozone concentration" "ppb"
>> 141 141 "NO" "Nitrogen Oxide" "ppbV"
>> 142 141 "NO2" "Nitrogen Dioxide" "ppbV"
>> 148 141 "CO" "Carbon Monoxide" "ppbV"
>>
>> The first column is the GRIB code and the second is the parameter
table version number. You can see in the attached
>> configuration file that these values are set for each parameter. I
ran point_stat with the attached configuration file
>> using the following command:
>>
>> bash$ $MET_BASE/bin/point_stat \
>> WRFPRS_d01.000 \
>> rama_2010.nc \
>> PointStatConfig_O3_org \
>> -outdir . -v 3 > ps_test.log 2>&1 &
>>
>> The output log is also attached. Note that I changed the
boot.n_rep setting from 1000 to 0 for testing. This feature
>> is not very practical in the situation where there are relatively
few matched pairs, as there are in this case. I also
>> modified the output_flag to write out only a single stat file, for
brevity.
>>
>> Please let me know if you have any questions about this
information.
>>
>> Paul
>>
>>
>> On 09/04/2012 07:28 PM, Jose Agustín García Reynoso via RT wrote:
>>>
>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=58094 >
>>>
>>> Paul,
>>>
>>> I uploaded the following files:
>>> PointStatConfig_O3_org
>>> WRFPRS_d01.000
>>> WRFPRS_d01.001
>>> rama_2010.txt
>>>
>>> The forecast are WRFPRS files
>>> The observation stations are rama_2001.txt in ascii you have to
convert in
>>> nc using
>>>
>>> bin/ascii2nc rama_2010.txt out/ascii2nc/simat_20100100.nc -v 1
>>>
>>> for the analysis:
>>>
>>> bin/point_stat WRFPRS_d01.000 out/ascii2nc/simat_20100100.nc
>>> PointStatConfig_O3_org
>>> -v 1
>>>
>>> or
>>>
>>> bin/point_stat WRFPRS_d01.001 out/ascii2nc/simat_20100100.nc
>>> PointStatConfig_O3_org
>>> -v 1
>>>
>>> Also we want to analyze the following compounds (NO, NO2, CO, SO2)
>>>
>>> fcst_field[] = [ "141/L1","142/L1", "148/L1", "232/L1" ];
>>> fcst_thresh[] = [ "ge105.0 ge210.0","ge105.0 ge210.0", "ge550.0
>>> ge1100.0","ge65 ge130" ];
>>>
>>> for grib_ptv = 141
>>>
>>> Thanks
>>>
>>>
>>> 2012/9/4 Paul Oldenburg via RT <met_help at ucar.edu>
>>>
>>>> Agustin,
>>>>
>>>> Can you please send me the data file that you are trying to
verify? That
>>>> might be the easiest way for me to track this
>>>> problem. It is possible that there is a bug in this
functionality, since
>>>> it is relatively new. If your data file is
>>>> larger then about 2MB, please upload it to our FTP site using the
>>>> instructions here:
>>>>
>>>> http://www.dtcenter.org/met/users/support/met_help.php#ftp
>>>>
>>>> If you have trouble, please let me know.
>>>>
>>>> Thanks,
>>>>
>>>> Paul
>>>>
>>>>
>>>> On 09/04/2012 03:21 PM, Jose Agustín García Reynoso via RT wrote:
>>>>>
>>>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=58094 >
>>>>>
>>>>> Paul,
>>>>>
>>>>> I did wat you suggest,
>>>>> I tried with
>>>>> GRIB1_ptv = 129;
>>>>> GRIB1_ptv = [129];
>>>>> and GRIB1_ptv = "129";
>>>>>
>>>>> but still having the table 2 as reference table.
>>>>>
>>>>>
>>>>> DEBUG 1: Default Config File:
>>>>>
/Users/pronos/Documents/WRF/METv4.0/data/config/PointStatConfig_default
>>>>> DEBUG 1: User Config File:
>>>>>
/Users/pronos/Documents/WRF/METv4.0/scripts/config/PointStatConfig_O3
>>>>> ERROR :
>>>>> ERROR : VarInfoGrib::set_dict() - unrecognized GRIB1 field
abbreviation
>>>>> 'OZCON' for table version 2
>>>>> ERROR :
>>>>>
>>>>> the log File using -v 6 is attached
>>>>>
>>>>> thanks
>>>>>
>>>>> 2012/9/4 Paul Oldenburg via RT <met_help at ucar.edu>
>>>>>
>>>>>> Agustin,
>>>>>>
>>>>>> I looked up the field you mentioned, table 129 variable 180.
In the
>>>>>> standard NCEP GRIB code tables, which is what MET
>>>>>> uses by default, this field corresponds to Ozone concentration
in PPB
>>>>>> abbreviated as "OZCON". See the web site:
>>>>>>
>>>>>> http://www.nco.ncep.noaa.gov/pmb/docs/on388/table2.html
>>>>>>
>>>>>> The information on this website is stored in MET in the table
>>>>>> $MET_BASE/data/table_files/nceptab_flat.txt. If you plan
>>>>>> to use different GRIB code tables, you should copy and rename
this file
>>>>>> and modify the original version to suit your
>>>>>> purposes. Then, MET will look for records in your GRIB files
based on
>>>> the
>>>>>> GRIB code abbreviations you specify in your
>>>>>> version of the file.
>>>>>>
>>>>>> If you are interested in verifying ozone concentration, then
use the
>>>>>> following settings in the point_stat configuration
>>>>>> file:
>>>>>>
>>>>>> {
>>>>>> name = "OZCON";
>>>>>> GRIB1_ptv = 129;
>>>>>> level = [ "L1" ];
>>>>>> cat_thresh = [ >=40 , >=60, >=110 , >=120 ];
>>>>>> }
>>>>>>
>>>>>> Please let me know if you have any questions or further
troubles.
>>>>>>
>>>>>> Thanks,
>>>>>>
>>>>>> Paul
>>>>>>
>>>>>>
>>>>>> On 09/01/2012 08:21 AM, Jose Agustín García Reynoso via RT
wrote:
>>>>>>>
>>>>>>> Sat Sep 01 08:21:20 2012: Request 58094 was acted upon.
>>>>>>> Transaction: Ticket created by agustin at atmosfera.unam.mx
>>>>>>> Queue: met_help
>>>>>>> Subject: grib_ptv
>>>>>>> Owner: Nobody
>>>>>>> Requestors: agustin at atmosfera.unam.mx
>>>>>>> Status: new
>>>>>>> Ticket <URL:
>>>> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=58094 >
>>>>>>>
>>>>>>>
>>>>>>> I update to METv4 + patches 20120801. I am trying to use a
different
>>>>>> table
>>>>>>> version 129 in order to use the variable 180 L1
>>>>>>> in the file METv4.0/data/config/README indicates the use of
GRIB1_ptv
>>>>>> and
>>>>>>> GRIB1_rec .
>>>>>>>
>>>>>>> in the config file is the following setup:
>>>>>>> ...
>>>>>>> {
>>>>>>> name = "180";
>>>>>>> level = [ "L1" ];
>>>>>>> cat_thresh = [ >=40 , >= 60, >= 110 , >= 120 ];
>>>>>>> GRIB1_ptv = [129];
>>>>>>> }
>>>>>>> ];
>>>>>>>
>>>>>>> but I get the following
>>>>>>>
>>>>>>> ./jal_point_stat.sh
>>>>>>>
>>>>>>> *** Running Point-Stat on Jalisco WRF data ***
>>>>>>> DEBUG 1: Default Config File:
>>>>>>>
/Users/pronos/Documents/WRF/METv4.0/data/config/PointStatConfig_default
>>>>>>> DEBUG 1: User Config File: config/PointStatConfig_O3
>>>>>>> ERROR :
>>>>>>> ERROR : VarInfoGrib::set_dict() - unrecognized GRIB1 field
>>>> abbreviation
>>>>>>> '180' for table version 2
>>>>>>> ERROR :
>>>>>>>
>>>>>>> when using the following
>>>>>>>
>>>>>>> {
>>>>>>> GRIB1_rec = [180];
>>>>>>> GRIB1_ptv = [129];
>>>>>>> level = [ "L1" ];
>>>>>>> cat_thresh = [ >=40 , >= 60, >= 110 , >= 120 ];
>>>>>>> }
>>>>>>>
>>>>>>> I got
>>>>>>>
>>>>>>> *** Running Point-Stat on Jalisco WRF data ***
>>>>>>> DEBUG 1: Default Config File:
>>>>>>>
/Users/pronos/Documents/WRF/METv4.0/data/config/PointStatConfig_default
>>>>>>> DEBUG 1: User Config File: config/PointStatConfig_O3
>>>>>>> ERROR :
>>>>>>> ERROR : VarInfoGrib::set_dict() - either name or GRIB1_ptv
and
>>>> GRIB1_rec
>>>>>>> must be specified in field information
>>>>>>> ERROR :
>>>>>>>
>>>>>>> How can I setup those variables?
>>>>>>>
>>>>>>> --
>>>>>>> Dr. Agustin Garcia-Reynoso
>>>>>>> Atmosfera <http://www.atmosfera.unam.mx/fqa>
>>>>>>> *Centro de Ciencias de la Atmosfera, *UNAM
>>>>>>> Tel (52) 55 *56224073* Fax (52) 55 56160789
>>>>>>> Circuito de la Investigacion Cientifica s/n Mexico DF 04510
>>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>
>>>>>
>>>>
>>>>
>>>>
>>>
>>>
>>
>>
>> <PointStatConfig_O3_org.txt><ps_test.log>
>
------------------------------------------------
More information about the Met_help
mailing list