[Met_help] [rt.rap.ucar.edu #48446] History for problem reading APCP in grid_stat
RAL HelpDesk {for Paul Oldenburg}
met_help at ucar.edu
Fri Jul 22 13:53:04 MDT 2011
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
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
----------------------------------------------------------------
More information about the Met_help
mailing list