[Met_help] [rt.rap.ucar.edu #58480] History for METv4 pcp_combine and pointstat
Paul Oldenburg via RT
met_help at ucar.edu
Wed Sep 26 12:35:37 MDT 2012
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
I am trying to run pointstat using gridded netcdf files generated by
pcp_combine and obs files from ascii2nc in METv4.0. I am unable to get
it to read the netcdf file properly. I am trying to modify the config
file two different ways to get this to work (shown below). Pointstat
either thinks I am giving it a grib file or I get a syntax error. I have
also listed the netcdf header of the file at the end. What am I doing
wrong here?
Command: point_stat wrfpcp_d01_2012092200_036.nc OBSASC_2012092200.nc
config/PointStatConfig -outdir ./ -v 2
Here is the 1st way I have the configured the config file:
fcst = {
wind_thresh = [ NA ];
message_type = [ "ADPSFC"];
field = [
{
name = "APCP_01";
level = [ "(*,*)" ];
cat_thresh = [ >0.0 ];
}
];
};
obs = fcst;
ERROR message: VarInfoGrib::set_dict() - unrecognized GRIB1 field
abbreviation 'APCP_01' for table version 2
(it thinks i am trying to pass it a grib file???)
=====================================================================
The 2nd way I am trying [as shown in the MET users guide] is this way: (
I get a syntax error here):
fcst_field [] = [ "APCP_01(*,*)" ];
obs = fcst;
ERROR message:
ERROR : yyerror() -> syntax error in file
"config/PointStatConfig_wrf_east12_d01_precip"
ERROR :
ERROR : line = 32
ERROR :
ERROR : column = 12
ERROR :
ERROR : text = "["
ERROR :
ERROR :
ERROR : fcst_field [] = [ "APCP_01(0,*,*)" ];
ncdump:
netcdf wrfpcp_d01_2012092200_036 {
dimensions:
lat = 151 ;
lon = 164 ;
variables:
float lat(lat, lon) ;
lat:long_name = "latitude" ;
lat:units = "degrees_north" ;
lat:standard_name = "latitude" ;
float lon(lat, lon) ;
lon:long_name = "longitude" ;
lon:units = "degrees_east" ;
lon:standard_name = "longitude" ;
float APCP_01(lat, lon) ;
APCP_01:name = "APCP_01" ;
APCP_01:long_name = "Total precipitation" ;
APCP_01:level = "A1" ;
APCP_01:units = "kg/m^2" ;
APCP_01:_FillValue = -9999.f ;
APCP_01:init_time = "20120922_000000" ;
APCP_01:init_time_ut = 1348272000 ;
APCP_01:valid_time = "20120923_120000" ;
APCP_01:valid_time_ut = 1348401600 ;
APCP_01:accum_time = "010000" ;
APCP_01:accum_time_sec = 3600 ;
// global attributes:
:FileOrigins = "File wrfpcp_d01_2012092200_036.nc generated
20120925_153004 UTC by the MET pcp_combine tool" ;
:MET_version = "V4.0" ;
:MET_tool = "pcp_combine" ;
:RunCommand = "Subtraction: wrfout_d01_2012092200_036.grb with
accumulation of 000000 minus wrfout_d01_2012092200_035.grb with
accumulation of 000000." ;
:Projection = "Lambert Conformal" ;
:scale_lat_1 = "60.000000" ;
:scale_lat_2 = "30.000000" ;
:lat_pin = "26.356000" ;
:lon_pin = "-89.670000" ;
:x_pin = "0.000000" ;
:y_pin = "0.000000" ;
:lon_orient = "-80.000000" ;
:d_km = "12.000000" ;
:r_km = "6371.200000" ;
:nx = "164" ;
:ny = "151 grid_points" ;
--
Aaron Sims
State Climate Office of North Carolina
North Carolina State University
Email: apsims at ncsu.edu
Phone: (919)513-2101
----------------------------------------------------------------
Complete Ticket History
----------------------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #58480] METv4 pcp_combine and pointstat
From: Paul Oldenburg
Time: Wed Sep 26 12:00:52 2012
Aaron,
I think what might be going on here is that you need to specify the
information for the observation field as if it was
in a GRIB file. I know this sounds confusing, but it is a bit of
legacy behavior. Try replacing this:
obs = fcst;
with this:
obs = {
wind_thresh = [ NA ];
message_type = [ "ADPSFC"];
field = [
{
name = "APCP";
level = [ "A1" ];
cat_thresh = [ >0.0 ];
}
];
};
If you still have trouble, could you please upload your data and
config file to our FTP site using the instructions here:
http://www.dtcenter.org/met/users/support/met_help.php#ftp
Then, I can take a closer look.
Thanks,
Paul
On 09/26/2012 11:54 AM, Aaron Sims via RT wrote:
>
> Wed Sep 26 11:54:26 2012: Request 58480 was acted upon.
> Transaction: Ticket created by apsims at ncsu.edu
> Queue: met_help
> Subject: METv4 pcp_combine and pointstat
> Owner: Nobody
> Requestors: apsims at ncsu.edu
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=58480 >
>
>
> I am trying to run pointstat using gridded netcdf files generated by
> pcp_combine and obs files from ascii2nc in METv4.0. I am unable to
get
> it to read the netcdf file properly. I am trying to modify the
config
> file two different ways to get this to work (shown below).
Pointstat
> either thinks I am giving it a grib file or I get a syntax error. I
have
> also listed the netcdf header of the file at the end. What am I
doing
> wrong here?
>
> Command: point_stat wrfpcp_d01_2012092200_036.nc
OBSASC_2012092200.nc
> config/PointStatConfig -outdir ./ -v 2
>
> Here is the 1st way I have the configured the config file:
> fcst = {
> wind_thresh = [ NA ];
> message_type = [ "ADPSFC"];
>
> field = [
> {
> name = "APCP_01";
> level = [ "(*,*)" ];
> cat_thresh = [ >0.0 ];
> }
> ];
>
> };
> obs = fcst;
>
> ERROR message: VarInfoGrib::set_dict() - unrecognized GRIB1 field
> abbreviation 'APCP_01' for table version 2
> (it thinks i am trying to pass it a grib file???)
>
=====================================================================
>
> The 2nd way I am trying [as shown in the MET users guide] is this
way: (
> I get a syntax error here):
> fcst_field [] = [ "APCP_01(*,*)" ];
> obs = fcst;
>
> ERROR message:
> ERROR : yyerror() -> syntax error in file
> "config/PointStatConfig_wrf_east12_d01_precip"
> ERROR :
> ERROR : line = 32
> ERROR :
> ERROR : column = 12
> ERROR :
> ERROR : text = "["
> ERROR :
> ERROR :
> ERROR : fcst_field [] = [ "APCP_01(0,*,*)" ];
>
>
>
> ncdump:
> netcdf wrfpcp_d01_2012092200_036 {
> dimensions:
> lat = 151 ;
> lon = 164 ;
> variables:
> float lat(lat, lon) ;
> lat:long_name = "latitude" ;
> lat:units = "degrees_north" ;
> lat:standard_name = "latitude" ;
> float lon(lat, lon) ;
> lon:long_name = "longitude" ;
> lon:units = "degrees_east" ;
> lon:standard_name = "longitude" ;
> float APCP_01(lat, lon) ;
> APCP_01:name = "APCP_01" ;
> APCP_01:long_name = "Total precipitation" ;
> APCP_01:level = "A1" ;
> APCP_01:units = "kg/m^2" ;
> APCP_01:_FillValue = -9999.f ;
> APCP_01:init_time = "20120922_000000" ;
> APCP_01:init_time_ut = 1348272000 ;
> APCP_01:valid_time = "20120923_120000" ;
> APCP_01:valid_time_ut = 1348401600 ;
> APCP_01:accum_time = "010000" ;
> APCP_01:accum_time_sec = 3600 ;
>
> // global attributes:
> :FileOrigins = "File wrfpcp_d01_2012092200_036.nc
generated
> 20120925_153004 UTC by the MET pcp_combine tool" ;
> :MET_version = "V4.0" ;
> :MET_tool = "pcp_combine" ;
> :RunCommand = "Subtraction: wrfout_d01_2012092200_036.grb
with
> accumulation of 000000 minus wrfout_d01_2012092200_035.grb with
> accumulation of 000000." ;
> :Projection = "Lambert Conformal" ;
> :scale_lat_1 = "60.000000" ;
> :scale_lat_2 = "30.000000" ;
> :lat_pin = "26.356000" ;
> :lon_pin = "-89.670000" ;
> :x_pin = "0.000000" ;
> :y_pin = "0.000000" ;
> :lon_orient = "-80.000000" ;
> :d_km = "12.000000" ;
> :r_km = "6371.200000" ;
> :nx = "164" ;
> :ny = "151 grid_points" ;
>
>
------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #58480] METv4 pcp_combine and pointstat
From: Aaron Sims
Time: Wed Sep 26 12:33:11 2012
Paul,
Thanks! That was it.
Aaron
On 09/26/2012 02:00 PM, Paul Oldenburg via RT wrote:
> Aaron,
>
> I think what might be going on here is that you need to specify the
information for the observation field as if it was
> in a GRIB file. I know this sounds confusing, but it is a bit of
legacy behavior. Try replacing this:
>
> obs = fcst;
>
> with this:
>
> obs = {
> wind_thresh = [ NA ];
> message_type = [ "ADPSFC"];
>
> field = [
> {
> name = "APCP";
> level = [ "A1" ];
> cat_thresh = [>0.0 ];
> }
> ];
>
> };
>
> If you still have trouble, could you please upload your data and
config file to our FTP site using the instructions here:
>
> http://www.dtcenter.org/met/users/support/met_help.php#ftp
>
> Then, I can take a closer look.
>
> Thanks,
>
> Paul
>
>
> On 09/26/2012 11:54 AM, Aaron Sims via RT wrote:
>> Wed Sep 26 11:54:26 2012: Request 58480 was acted upon.
>> Transaction: Ticket created by apsims at ncsu.edu
>> Queue: met_help
>> Subject: METv4 pcp_combine and pointstat
>> Owner: Nobody
>> Requestors: apsims at ncsu.edu
>> Status: new
>> Ticket<URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=58480>
>>
>>
>> I am trying to run pointstat using gridded netcdf files generated
by
>> pcp_combine and obs files from ascii2nc in METv4.0. I am unable to
get
>> it to read the netcdf file properly. I am trying to modify the
config
>> file two different ways to get this to work (shown below).
Pointstat
>> either thinks I am giving it a grib file or I get a syntax error. I
have
>> also listed the netcdf header of the file at the end. What am I
doing
>> wrong here?
>>
>> Command: point_stat wrfpcp_d01_2012092200_036.nc
OBSASC_2012092200.nc
>> config/PointStatConfig -outdir ./ -v 2
>>
>> Here is the 1st way I have the configured the config file:
>> fcst = {
>> wind_thresh = [ NA ];
>> message_type = [ "ADPSFC"];
>>
>> field = [
>> {
>> name = "APCP_01";
>> level = [ "(*,*)" ];
>> cat_thresh = [>0.0 ];
>> }
>> ];
>>
>> };
>> obs = fcst;
>>
>> ERROR message: VarInfoGrib::set_dict() - unrecognized GRIB1 field
>> abbreviation 'APCP_01' for table version 2
>> (it thinks i am trying to pass it a grib file???)
>>
=====================================================================
>>
>> The 2nd way I am trying [as shown in the MET users guide] is this
way: (
>> I get a syntax error here):
>> fcst_field [] = [ "APCP_01(*,*)" ];
>> obs = fcst;
>>
>> ERROR message:
>> ERROR : yyerror() -> syntax error in file
>> "config/PointStatConfig_wrf_east12_d01_precip"
>> ERROR :
>> ERROR : line = 32
>> ERROR :
>> ERROR : column = 12
>> ERROR :
>> ERROR : text = "["
>> ERROR :
>> ERROR :
>> ERROR : fcst_field [] = [ "APCP_01(0,*,*)" ];
>>
>>
>>
>> ncdump:
>> netcdf wrfpcp_d01_2012092200_036 {
>> dimensions:
>> lat = 151 ;
>> lon = 164 ;
>> variables:
>> float lat(lat, lon) ;
>> lat:long_name = "latitude" ;
>> lat:units = "degrees_north" ;
>> lat:standard_name = "latitude" ;
>> float lon(lat, lon) ;
>> lon:long_name = "longitude" ;
>> lon:units = "degrees_east" ;
>> lon:standard_name = "longitude" ;
>> float APCP_01(lat, lon) ;
>> APCP_01:name = "APCP_01" ;
>> APCP_01:long_name = "Total precipitation" ;
>> APCP_01:level = "A1" ;
>> APCP_01:units = "kg/m^2" ;
>> APCP_01:_FillValue = -9999.f ;
>> APCP_01:init_time = "20120922_000000" ;
>> APCP_01:init_time_ut = 1348272000 ;
>> APCP_01:valid_time = "20120923_120000" ;
>> APCP_01:valid_time_ut = 1348401600 ;
>> APCP_01:accum_time = "010000" ;
>> APCP_01:accum_time_sec = 3600 ;
>>
>> // global attributes:
>> :FileOrigins = "File wrfpcp_d01_2012092200_036.nc
generated
>> 20120925_153004 UTC by the MET pcp_combine tool" ;
>> :MET_version = "V4.0" ;
>> :MET_tool = "pcp_combine" ;
>> :RunCommand = "Subtraction:
wrfout_d01_2012092200_036.grb with
>> accumulation of 000000 minus wrfout_d01_2012092200_035.grb with
>> accumulation of 000000." ;
>> :Projection = "Lambert Conformal" ;
>> :scale_lat_1 = "60.000000" ;
>> :scale_lat_2 = "30.000000" ;
>> :lat_pin = "26.356000" ;
>> :lon_pin = "-89.670000" ;
>> :x_pin = "0.000000" ;
>> :y_pin = "0.000000" ;
>> :lon_orient = "-80.000000" ;
>> :d_km = "12.000000" ;
>> :r_km = "6371.200000" ;
>> :nx = "164" ;
>> :ny = "151 grid_points" ;
>>
>>
>
--
Aaron Sims
State Climate Office of North Carolina
North Carolina State University
Email: apsims at ncsu.edu
Phone: (919)513-2101
------------------------------------------------
More information about the Met_help
mailing list