[gradsusr] Precipitation Output Files Changed from NOAA AHPS: Can't read in GrADS anoymore
James Spinks
jspin881 at hotmail.com
Sun Jul 30 00:08:56 EDT 2017
Hi James,
Thanks for the response. I was able to get the data to display, but had to add the xdef information because GrADS wouldn't open the file with out it. My only issue is that I do not believe the xdef and ydef information in the ctl file is accurate even though the file diplayed ok. Here is the modified ctl file:
DSET ^nws_precip_1day_20170728_conus.nc
dtype netcdf
options yrev
title amountofprecip
undef -3.40282346638529e+38
pdef 1051 813 nps 388 1591 -105.000000 4.7625
ydef 881 linear 20 0.05
xdef 1121 linear -125 0.05
zdef 1 linear 1 1
tdef 1 linear 00Z28Jul2017 1dy
vars 4
observation=>obs 1 y,x observation
normal=>pcp 1 y,x normal
departure_from_normal=>dpcp 1 y,x departure
percent_of_normal=>ppcp 1 y,x percent
endvars
This ctl displays ok, but I'm not certain if the xdef and ydef info is accurate. Can you or anyone else please confirm? Thanks
________________________________
From: James T. Potemra <jimp at hawaii.edu>
Sent: Friday, July 28, 2017 10:35 PM
To: James Spinks
Cc: GrADS Users Forum
Subject: Re: [gradsusr] Precipitation Output Files Changed from NOAA AHPS: Can't read in GrADS anoymore
Hi James,
I was able to use ncdump on this file. There is no variable called "amountofprecip", but instead there is "normal", "departure_from_normal", and "percent_of_normal". The complete listing is below. The following ctl works with this file:
DSET ^nws_precip_1day_20170728_conus.nc
dtype netcdf
options yrev
title amountofprecip
undef -3.40282346638529e+38
pdef 1051 813 nps 388 1591 -105.000000 4.7625
ydef 700 linear 20 0.05
zdef 1 linear 1 1
tdef 1 linear 00Z28Jul2017 1dy
vars 4
observation=>obs 1 y,x observation
normal=>pcp 1 y,x normal
departure_from_normal=>dpcp 1 y,x departure
percent_of_normal=>ppcp 1 y,x percent
endvars
Jim
netcdf nws_precip_1day_20170728_conus {
dimensions:
y = 881 ;
x = 1121 ;
polar_stereographic = UNLIMITED ; // (0 currently)
variables:
int polar_stereographic(polar_stereographic) ;
polar_stereographic:latitude_of_projection_origin = 90. ;
polar_stereographic:straight_vertical_longitude_from_pole = -105. ;
polar_stereographic:standard_parallel = 60. ;
polar_stereographic:false_northing = 0. ;
polar_stereographic:grid_mapping_name = "polar_stereographic" ;
polar_stereographic:false_easting = 0. ;
float observation(y, x) ;
observation:_FillValue = -3.402823e+38f ;
observation:proj4 = "+proj=stere +lat_0=90 +lat_ts=60 +lon_0=-105 +k=90 +x_0=0 +y_0=0 +a=6371200 +b=6371200 +units=m +no_defs" ;
observation:grid_mapping = "polar_stereographic" ;
observation:name = "observation" ;
observation:coordinates = "x y" ;
observation:long_name = "Observation" ;
observation:units = "Inches" ;
observation:missing_value = -3.40282346638529e+38 ;
observation:esri_pe_string = "PROJCS[\"Polar_Stereographic\",GEOGCS[\"GCS_Coordinate_System_imported_from_GRIB_file\",DATUM[\"D_unknown\",SPHEROID[\"Sphere\",6371200.0,0.0]],PRIMEM[\"Greenwich\",0.0],UNIT[\"Degree\",0.0174532925199433]],PROJECTION[\"Stereographic_North_Pole\"],PARAMETER[\"false_easting\",0.0],PARAMETER[\"false_northing\",0.0],PARAMETER[\"central_meridian\",-105.0],PARAMETER[\"standard_parallel_1\",60.0],UNIT[\"Meter\",1.0]]" ;
float normal(y, x) ;
normal:_FillValue = -3.402823e+38f ;
normal:proj4 = "+proj=stere +lat_0=90 +lat_ts=60 +lon_0=-105 +k=90 +x_0=0 +y_0=0 +a=6371200 +b=6371200 +units=m +no_defs" ;
normal:grid_mapping = "polar_stereographic" ;
normal:name = "normal" ;
normal:coordinates = "x y" ;
normal:long_name = "PRISM normals" ;
normal:units = "Inches" ;
normal:missing_value = -3.40282346638529e+38 ;
normal:esri_pe_string = "PROJCS[\"Polar_Stereographic\",GEOGCS[\"GCS_Coordinate_System_imported_from_GRIB_file\",DATUM[\"D_unknown\",SPHEROID[\"Sphere\",6371200.0,0.0]],PRIMEM[\"Greenwich\",0.0],UNIT[\"Degree\",0.0174532925199433]],PROJECTION[\"Stereographic_North_Pole\"],PARAMETER[\"false_easting\",0.0],PARAMETER[\"false_northing\",0.0],PARAMETER[\"central_meridian\",-105.0],PARAMETER[\"standard_parallel_1\",60.0],UNIT[\"Meter\",1.0]]" ;
float departure_from_normal(y, x) ;
departure_from_normal:_FillValue = -3.402823e+38f ;
departure_from_normal:proj4 = "+proj=stere +lat_0=90 +lat_ts=60 +lon_0=-105 +k=90 +x_0=0 +y_0=0 +a=6371200 +b=6371200 +units=m +no_defs" ;
departure_from_normal:grid_mapping = "polar_stereographic" ;
departure_from_normal:name = "departure_from_normal" ;
departure_from_normal:coordinates = "x y" ;
departure_from_normal:long_name = "Departure from normal" ;
departure_from_normal:units = "Inches" ;
departure_from_normal:missing_value = -3.40282346638529e+38 ;
departure_from_normal:esri_pe_string = "PROJCS[\"Polar_Stereographic\",GEOGCS[\"GCS_Coordinate_System_imported_from_GRIB_file\",DATUM[\"D_unknown\",SPHEROID[\"Sphere\",6371200.0,0.0]],PRIMEM[\"Greenwich\",0.0],UNIT[\"Degree\",0.0174532925199433]],PROJECTION[\"Stereographic_North_Pole\"],PARAMETER[\"false_easting\",0.0],PARAMETER[\"false_northing\",0.0],PARAMETER[\"central_meridian\",-105.0],PARAMETER[\"standard_parallel_1\",60.0],UNIT[\"Meter\",1.0]]" ;
float percent_of_normal(y, x) ;
percent_of_normal:_FillValue = -3.402823e+38f ;
percent_of_normal:proj4 = "+proj=stere +lat_0=90 +lat_ts=60 +lon_0=-105 +k=90 +x_0=0 +y_0=0 +a=6371200 +b=6371200 +units=m +no_defs" ;
percent_of_normal:grid_mapping = "polar_stereographic" ;
percent_of_normal:name = "percent_of_normal" ;
percent_of_normal:coordinates = "x y" ;
percent_of_normal:long_name = "Percent of normal" ;
percent_of_normal:units = "Percent" ;
percent_of_normal:missing_value = -3.40282346638529e+38 ;
percent_of_normal:esri_pe_string = "PROJCS[\"Polar_Stereographic\",GEOGCS[\"GCS_Coordinate_System_imported_from_GRIB_file\",DATUM[\"D_unknown\",SPHEROID[\"Sphere\",6371200.0,0.0]],PRIMEM[\"Greenwich\",0.0],UNIT[\"Degree\",0.0174532925199433]],PROJECTION[\"Stereographic_North_Pole\"],PARAMETER[\"false_easting\",0.0],PARAMETER[\"false_northing\",0.0],PARAMETER[\"central_meridian\",-105.0],PARAMETER[\"standard_parallel_1\",60.0],UNIT[\"Meter\",1.0]]" ;
double x(x) ;
x:units = "Meter" ;
x:long_name = "x coordinate of projection" ;
x:standard_name = "projection_x_coordinate" ;
double y(y) ;
y:units = "Meter" ;
y:long_name = "y coordinate of projection" ;
y:standard_name = "projection_y_coordinate" ;
// global attributes:
:creation_time = "20170729T001638" ;
:Conventions = "CF-1.6" ;
}
On 7/28/17 1:52 PM, James Spinks wrote:
Hi all,
NOAA AHPS changed their netcdf for format for the precipitation data as stated here: http://www.nws.noaa.gov/om/notification/scn17-32ahps_pcpnaaa.htm
will be available in the following formats:
- Climate and Forecast (CF)-Network Common Data Form (netCDF)
Precipitation products will no longer be available in the following
formats:
- netCDF not adhering to the netCDF CF metadata convention
I'm not exactly sure the difference of the two and how it effects grads. I previously had a .ctl file the would read the original netcdf file just fine:
DSET ^nws_precip_1day_20170728_conus.nc
dtype netcdf
title amountofprecip
undef -1
pdef 1051 813 nps 388 1591 -105.000000 4.7625
xdef 1500 linear -135 0.05
ydef 700 linear 20 0.05
zdef 1 linear 1 1
tdef 1 linear 00Z28Jul2017 1dy
vars 1
amountofprecip=>pcp 1 y,x amountofprecip
endvars
Now this ctl does not work nor does ncdump read the new file from http://water.weather.gov/precip/downloads/
Has anyone had any success opening the new updated precip data in GrADS from AHPS? If so, do you have a ctl file you could share?
Thanks,
James Spinks
_______________________________________________
gradsusr mailing list
gradsusr at gradsusr.org<mailto:gradsusr at gradsusr.org>
http://gradsusr.org/mailman/listinfo/gradsusr
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20170730/eaa71eb1/attachment.html
More information about the gradsusr
mailing list