[gradsusr] HRAP precipitation grid - UPGRADED 2018 version
yaqiang.wang at gmail.com
yaqiang.wang at gmail.com
Tue Oct 23 02:15:07 EDT 2018
You can try
PDEF 1121 881 nps 399 1599 -105.0 4.763
Regards
Yaqiang
*************************************************
Dr. Yaqiang Wang
Chinese Academy of Meteorological Sciences (CAMS)
46, Zhong-Guan-Cun South Avenue
Beijing, 100081
China
yaqiang.wang at gmail.com
http://www.meteothink.org
https://www.researchgate.net/profile/Yaqiang_Wang5
**************************************************
From: Kevin Levey
Date: 2018-10-23 03:38
To: GrADS Users Forum; Jennifer M Adams
Subject: [gradsusr] HRAP precipitation grid - UPGRADED 2018 version
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
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)
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/20140101.txt
Location of the NWS Precip page:
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
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/20181023/7e2cd2fe/attachment-0001.html>
More information about the gradsusr
mailing list