[gradsusr] Precipitation Output Files Changed from NOAA AHPS: Can't read in GrADS anoymore
James T. Potemra
jimp at hawaii.edu
Fri Jul 28 22:35:07 EDT 2017
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
>
>
> willbe 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/
> <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
> http://gradsusr.org/mailman/listinfo/gradsusr
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20170728/2a8673aa/attachment-0001.html
More information about the gradsusr
mailing list