[Met_help] [rt.rap.ucar.edu #48451] History for RE: problem reading APCP in grid_stat
RAL HelpDesk {for John Halley Gotway}
met_help at ucar.edu
Tue Jul 19 16:02:59 MDT 2011
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
Hello met_help,
FYI, I was able to fix my problem. I realized that newer versions of MET/pcp_combine have appended the precip accumulation onto the name of the variable in the netcdf output (i.e. APCP becomes APCP_01 for 1-hour accumulated precip). I have adjusted my script to construct the APCP variable more generically based on the accumulation interval, such that:
fcst_field[] = [ "APCP_01/A01" ];
obs_field[] = [ "APCP_01/A01" ];
The above entry now works in grid_stat.
Hopefully, you didn't work too much on a reply, since I figured it out myself!
Jonathan
From: Case, Jonathan (MSFC-VP61)[ENSCO INC]
Sent: Tuesday, July 19, 2011 2:30 PM
To: 'met_help at ucar.edu'
Subject: problem reading APCP in grid_stat
Hello met_help,
I am upgrading to the current MET v3 software, updating my scripts in the process.
I believe I have everything set to go for grid_stat, except that the program is giving me an error as follows:
discover08:/discover/nobackup/jlcase/MET/GRIDSTAT > grid_stat WRFPRECIP_1h OBSPRECIP_1h /home/jlcase/MET/CURRENT/data/config/GridStatConfig_APCP_CONUS -outdir /discover/nobackup/jlcase/MET/GRIDSTAT/SPORT_rt_NCEP/2010071700 -v 2
GSL_RNG_TYPE=mt19937
GSL_RNG_SEED=752362830
Forecast File: WRFPRECIP_1h
Observation File: OBSPRECIP_1h
Configuration File: /home/jlcase/MET/CURRENT/data/config/GridStatConfig_APCP_CONUS
***WARNING***: process_scores() -> APCP(*,*) not found in file: WRFPRECIP_1h
----------------------
In the config file, I have the following set for a 1-hour precip verification:
fcst_field[] = [ "APCP/A1" ];
obs_field[] = [ "APCP/A1" ];
Finally, my script re-projects the WRF precip grids onto the Stage IV grid prior to subtracting via pcp_combine.
I am able to successfully obtain a fcst and obs netcdf file for the valid 1-hour precip window using copygb. However, it appears that grid_stat doesn't recognize the APCP grid in the netcdf file. Here is a "netcdf -h" of the fcst netcdf file that grid_stat is complaining about, just to confirm that the data are there and named properly:
discover08:/discover/nobackup/jlcase/MET/GRIDSTAT > ncdump -h WRFPRECIP_1h
netcdf WRFPRECIP_1h {
dimensions:
lat = 881 ;
lon = 1121 ;
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" ;
APCP_01:long_name = "Total precipitation" ;
APCP_01:level = "A1" ;
APCP_01:units = "kg/m^2" ;
APCP_01:grib_code = 61 ;
APCP_01:_FillValue = -9999.f ;
APCP_01:init_time = "20100717_000000" ;
APCP_01:init_time_ut = 1279324800 ;
APCP_01:valid_time = "20100717_120000" ;
APCP_01:valid_time_ut = 1279368000 ;
APCP_01:accum_time = "010000" ;
APCP_01:accum_time_sec = 3600 ;
// global attributes:
:FileOrigins = "File WRFPRECIP_1h generated 20110719_190238 UTC on host discover08 by the MET pcp_combine tool" ;
:MET_version = "V3.0" ;
:MET_tool = "pcp_combine" ;
:RunCommand = "Subtraction: WRFPRS012.tm00.d01.grb.tostiv with accumulation of 120000 minus WRFPRS011.tm00.d01.grb.tostiv with accumulation of 110000." ;
:Projection = "Polar Stereographic" ;
:hemisphere = "N" ;
:scale_lat = "60.000000 degrees_north" ;
:lat_pin = "23.117000" ;
:lon_pin = "-119.023000" ;
:x_pin = "0.000000" ;
:y_pin = "0.000000" ;
:lon_orient = "-105.000000" ;
:d_km = "4.763000 km" ;
:r_km = "6367.470000 km" ;
:nx = "1121" ;
:ny = "881" ;
}
Any ideas?
Jonathan
--------------------------------------------------------------------------------------------------
Jonathan Case, ENSCO Inc.
NASA Short-term Prediction Research and Transition Center (aka SPoRT Center)
320 Sparkman Drive, Room 3062
Huntsville, AL 35805
Voice: 256.961.7504
Fax: 256.961.7788
Emails: Jonathan.Case-1 at nasa.gov / case.jonathan at ensco.com
--------------------------------------------------------------------------------------------------
"Whether the weather is cold, or whether the weather is hot, we'll weather
the weather whether we like it or not!"
----------------------------------------------------------------
Complete Ticket History
----------------------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #48451] RE: problem reading APCP in grid_stat
From: John Halley Gotway
Time: Tue Jul 19 16:02:44 2011
Jon,
Great - glad you were able to figure it out!
John
On 07/19/2011 02:33 PM, RAL HelpDesk {for Case, Jonathan[ENSCO INC]}
wrote:
>
> Tue Jul 19 14:33:23 2011: Request 48451 was acted upon.
> Transaction: Ticket created by jonathan.case-1 at nasa.gov
> Queue: met_help
> Subject: RE: problem reading APCP in grid_stat
> Owner: Nobody
> Requestors: jonathan.case-1 at nasa.gov
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=48451 >
>
>
> Hello met_help,
>
> FYI, I was able to fix my problem. I realized that newer versions
of MET/pcp_combine have appended the precip accumulation onto the name
of the variable in the netcdf output (i.e. APCP becomes APCP_01 for 1-
hour accumulated precip). I have adjusted my script to construct the
APCP variable more generically based on the accumulation interval,
such that:
>
> fcst_field[] = [ "APCP_01/A01" ];
> obs_field[] = [ "APCP_01/A01" ];
>
> The above entry now works in grid_stat.
>
> Hopefully, you didn't work too much on a reply, since I figured it
out myself!
> Jonathan
>
> From: Case, Jonathan (MSFC-VP61)[ENSCO INC]
> Sent: Tuesday, July 19, 2011 2:30 PM
> To: 'met_help at ucar.edu'
> Subject: problem reading APCP in grid_stat
>
> Hello met_help,
>
> I am upgrading to the current MET v3 software, updating my scripts
in the process.
> I believe I have everything set to go for grid_stat, except that the
program is giving me an error as follows:
>
> discover08:/discover/nobackup/jlcase/MET/GRIDSTAT > grid_stat
WRFPRECIP_1h OBSPRECIP_1h
/home/jlcase/MET/CURRENT/data/config/GridStatConfig_APCP_CONUS -outdir
/discover/nobackup/jlcase/MET/GRIDSTAT/SPORT_rt_NCEP/2010071700 -v 2
> GSL_RNG_TYPE=mt19937
> GSL_RNG_SEED=752362830
> Forecast File: WRFPRECIP_1h
> Observation File: OBSPRECIP_1h
> Configuration File:
/home/jlcase/MET/CURRENT/data/config/GridStatConfig_APCP_CONUS
> ***WARNING***: process_scores() -> APCP(*,*) not found in file:
WRFPRECIP_1h
> ----------------------
>
> In the config file, I have the following set for a 1-hour precip
verification:
> fcst_field[] = [ "APCP/A1" ];
> obs_field[] = [ "APCP/A1" ];
>
>
> Finally, my script re-projects the WRF precip grids onto the Stage
IV grid prior to subtracting via pcp_combine.
> I am able to successfully obtain a fcst and obs netcdf file for the
valid 1-hour precip window using copygb. However, it appears that
grid_stat doesn't recognize the APCP grid in the netcdf file. Here is
a "netcdf -h" of the fcst netcdf file that grid_stat is complaining
about, just to confirm that the data are there and named properly:
> discover08:/discover/nobackup/jlcase/MET/GRIDSTAT > ncdump -h
WRFPRECIP_1h
> netcdf WRFPRECIP_1h {
> dimensions:
> lat = 881 ;
> lon = 1121 ;
> 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" ;
> APCP_01:long_name = "Total precipitation" ;
> APCP_01:level = "A1" ;
> APCP_01:units = "kg/m^2" ;
> APCP_01:grib_code = 61 ;
> APCP_01:_FillValue = -9999.f ;
> APCP_01:init_time = "20100717_000000" ;
> APCP_01:init_time_ut = 1279324800 ;
> APCP_01:valid_time = "20100717_120000" ;
> APCP_01:valid_time_ut = 1279368000 ;
> APCP_01:accum_time = "010000" ;
> APCP_01:accum_time_sec = 3600 ;
>
> // global attributes:
> :FileOrigins = "File WRFPRECIP_1h generated
20110719_190238 UTC on host discover08 by the MET pcp_combine tool" ;
> :MET_version = "V3.0" ;
> :MET_tool = "pcp_combine" ;
> :RunCommand = "Subtraction:
WRFPRS012.tm00.d01.grb.tostiv with accumulation of 120000 minus
WRFPRS011.tm00.d01.grb.tostiv with accumulation of 110000." ;
> :Projection = "Polar Stereographic" ;
> :hemisphere = "N" ;
> :scale_lat = "60.000000 degrees_north" ;
> :lat_pin = "23.117000" ;
> :lon_pin = "-119.023000" ;
> :x_pin = "0.000000" ;
> :y_pin = "0.000000" ;
> :lon_orient = "-105.000000" ;
> :d_km = "4.763000 km" ;
> :r_km = "6367.470000 km" ;
> :nx = "1121" ;
> :ny = "881" ;
> }
>
> Any ideas?
> Jonathan
>
>
--------------------------------------------------------------------------------------------------
> Jonathan Case, ENSCO Inc.
> NASA Short-term Prediction Research and Transition Center (aka SPoRT
Center)
> 320 Sparkman Drive, Room 3062
> Huntsville, AL 35805
> Voice: 256.961.7504
> Fax: 256.961.7788
> Emails: Jonathan.Case-1 at nasa.gov / case.jonathan at ensco.com
>
--------------------------------------------------------------------------------------------------
>
> "Whether the weather is cold, or whether the weather is hot, we'll
weather
> the weather whether we like it or not!"
>
------------------------------------------------
More information about the Met_help
mailing list