[gradsusr] HRAP precipitation grid - UPGRADED 2018 version
Kevin Levey
klevey at customweather.com
Mon Oct 22 15:38:22 EDT 2018
MON 22OCT18: 1230PDT
Hi Jennifer,
I was wondering if you could help me here. The HRAP data was updated on 27th JUN this year.
The CTL file that previously you helped Ryan with (below) no longer works.
According to information from the HRAP page:
https://water.weather.gov/precip/download.php <https://water.weather.gov/precip/download.php>
The new QPE NetCDFs generated from the NCEP Stage IV data are CF-compliant multi-variable NetCDFs. The variables they contain are:
Band 1 - Observation- Last 24 hours of QPE spanning 12Z to 12Z in inches
Band 2 - PRISM normals- PRISM normals in inches (see "Normal Precipitation" section on the About page) <https://water.weather.gov/precip/about.php>
Band 3 - Departure from normal- The departure from normal in inches
Band 4 - Percent of normal- The percent of normal
Here is output from today’s 1-day conus file:
ncdump -h /mnt/modeldata/domModel/HRAP/rainfall/2018/10/22/nws_precip_1day_conus.nc
netcdf nws_precip_1day_conus {
dimensions:
y = 881 ;
x = 1121 ;
variables:
int crs ;
crs:straight_vertical_longitude_from_pole = -105. ;
crs:false_easting = 0. ;
crs:standard_parallel = 60. ;
crs:false_northing = 0. ;
crs:grid_mapping_name = "polar_stereographic" ;
crs:latitude_of_projection_origin = 90. ;
crs:earth_radius = 6371200. ;
crs:proj4 = "+proj=stere +lat_0=90 +lat_ts=60 +lon_0=-105 +x_0=0 +y_0=0 +a=6371200 +b=6371200 +units=m +no_defs" ;
crs:crs_wkt = "PROJCS[\"NOAA_HRAP_Grid\", GEOGCS[\"GCS_NOAA_HRAP\", DATUM[\"D_NOAA_HRAP\", SPHEROID[\"Sphere\",6371200,0.0] ], PRIMEM[\"Greenwich\",0.0], UNIT[\"Degree\",0.0174532925199433] ], PROJECTION[\"Polar_Stereographic\"], PARAMETER[\"false_easting\",0.0], PARAMETER[\"false_northing\",0.0], PARAMETER[\"central_meridian\",-105.0], PARAMETER[\"latitude_of_origin\",60.0], UNIT[\"meter\",1.0]]" ;
crs:esri_pe_string = "PROJCS[\"NOAA_HRAP_Grid\", GEOGCS[\"GCS_NOAA_HRAP\", DATUM[\"D_NOAA_HRAP\", SPHEROID[\"Sphere\",6371200,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 observation(y, x) ;
observation:_FillValue = -3.402823e+38f ;
observation:proj4 = "+proj=stere +lat_0=90 +lat_ts=60 +lon_0=-105 +x_0=0 +y_0=0 +a=6371200 +b=6371200 +units=m +no_defs" ;
observation:grid_mapping = "crs" ;
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[\"NOAA_HRAP_Grid\", GEOGCS[\"GCS_NOAA_HRAP\", DATUM[\"D_NOAA_HRAP\", SPHEROID[\"Sphere\",6371200,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 +x_0=0 +y_0=0 +a=6371200 +b=6371200 +units=m +no_defs" ;
normal:grid_mapping = "crs" ;
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[\"NOAA_HRAP_Grid\", GEOGCS[\"GCS_NOAA_HRAP\", DATUM[\"D_NOAA_HRAP\", SPHEROID[\"Sphere\",6371200,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 +x_0=0 +y_0=0 +a=6371200 +b=6371200 +units=m +no_defs" ;
departure_from_normal:grid_mapping = "crs" ;
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[\"NOAA_HRAP_Grid\", GEOGCS[\"GCS_NOAA_HRAP\", DATUM[\"D_NOAA_HRAP\", SPHEROID[\"Sphere\",6371200,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 +x_0=0 +y_0=0 +a=6371200 +b=6371200 +units=m +no_defs" ;
percent_of_normal:grid_mapping = "crs" ;
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[\"NOAA_HRAP_Grid\", GEOGCS[\"GCS_NOAA_HRAP\", DATUM[\"D_NOAA_HRAP\", SPHEROID[\"Sphere\",6371200,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 = "20181022T185215" ;
:data_time = "20181022" ;
:Conventions = "CF-1.6" ;
}
I changed
>> amountofprecip=>pcp 1 y,x amountofprecip
to
observation=>precip 1 y,x Precip Observation
Grads runs and produces this output:
ga-> open /mnt/modeldata/domModel/HRAP/rainfall/grads_ctl/hrap_conus.ctl
Scanning description file: /mnt/modeldata/domModel/HRAP/rainfall/grads_ctl/hrap_conus.ctl
Data file /mnt/modeldata/domModel/HRAP/rainfall/2018/10/22/nws_precip_1day_conus.nc is open as file 1
LON set to -135 -60.05
LAT set to 20 54.95
LEV set to 1 1
Time values set: 2018:10:22:0 2018:10:22:0
E set to 1 1
Notice: Implied interpolation for file /mnt/modeldata/domModel/HRAP/rainfall/grads_ctl/hrap_conus.ctl
Interpolation will be performed on any data displayed from this file
ga-> q file 1
File 1 : Gridded 4km x 4km HRAP 24-Hour Observed Rainfall ending 8AM ET
Descriptor: /mnt/modeldata/domModel/HRAP/rainfall/grads_ctl/hrap_conus.ctl
Binary: /mnt/modeldata/domModel/HRAP/rainfall/2018/10/22/nws_precip_1day_conus.nc
Type = Gridded
Xsize = 1500 Ysize = 700 Zsize = 1 Tsize = 1 Esize = 1
Number of Variables = 1
precip 1 y,x Precip Observation
ga-> set gxout shaded
ga-> d precip
Notice: Automatic Grid Interpolation Taking Place
Contouring: -3.3e+38 to 0 interval 3e+37
but nothing is plotted - the map shows no data, only USA outline.
Here is my control file:
DSET /mnt/modeldata/domModel/HRAP/rainfall/2018/10/22/nws_precip_1day_conus.nc
DTYPE netcdf
TITLE Gridded 4km x 4km HRAP 24-Hour Observed Rainfall ending 8AM ET
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 00:00Z22oct2018 1dy
VARS 1
observation=>precip 1 y,x Precip Observation
ENDVARS
Any help here would be much appreciated.
Cheers
Kevin
> On Apr 18, 2014, at 10:43 AM, Jennifer Adams <jma at COLA.IGES.ORG> wrote:
>
> Hi, Ryan --
> I subtracted hrap.xor from 401 and hrap.yor from 1601 and now it seems to fit with this PDEF:
> pdef 1051 813 nps 388 1591 -105.000000 4.7625
> xdef 1500 linear -135 0.05
> ydef 700 linear 20 0.05
>
> Note the new XDEF and YDEF axes do not need to match the size of the original grid; I used 0.05 degrees as a grid increment which kind of matches the native grid size, and made the axes big enough to cover the whole domain:
>
> <Screen Shot 2014-04-18 at 1.42.01 PM.png>
>
>
> --Jennifer
>
> On Apr 9, 2014, at 5:14 PM, Maue, Ryan wrote:
>
>> Greetings,
>>
>> Does anyone have a control file handy for the netcdf files that are provided by NWS AHPS? My control file is a little bit off. The grid is the HRAP 4 km polar stereographic.
>>
>> Any help or guidance would be great!
>>
>> I've uploaded an example here as well as the ascii output from the conversion program:
>> http://models.weatherbell.com/news/nws_precip_conus_20140101.nc <http://models.weatherbell.com/news/nws_precip_conus_20140101.nc>
>> http://models.weatherbell.com/news/20140101.txt <http://models.weatherbell.com/news/20140101.txt>
>>
>> Location of the NWS Precip page:
>> http://water.weather.gov/precip/ <http://water.weather.gov/precip/>
>>
>> dset ^nws_precip_conus_20140101.nc
>> dtype netcdf
>> title amountofprecip
>> undef -1
>> pdef 1051 813 nps 401 1601 -105.000000 4.7625
>> xdef 1051 linear -131.402908 .06366019047619047619
>> ydef 813 linear 20.761557 .04260785837438423645
>> ydef 813 linear 20.749580 0.04262164576
>> zdef 1 linear 1 1
>> tdef 1 linear 00:00Z01jan2014 1dy
>> vars 1
>> amountofprecip=>pcp 1 y,x amountofprecip
>> endvars
>>
>> RYAN MAUE
>>
>>
>>
>>
>>
>>
>>
>>
>> _______________________________________________
>> gradsusr mailing list
>> gradsusr at gradsusr.org <mailto:gradsusr at gradsusr.org>
>> http://gradsusr.org/mailman/listinfo/gradsusr <http://gradsusr.org/mailman/listinfo/gradsusr>
>
> --
> Jennifer M. Adams
> Center for Ocean-Land-Atmosphere Studies (COLA)
> 111 Research Hall, Mail Stop 2B3
> George Mason University
> 4400 University Drive
> Fairfax, VA 22030
>
>
>
>
>
> _______________________________________________
> 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/20181022/fc99d1d5/attachment-0001.html>
More information about the gradsusr
mailing list