[gradsusr] HRAP precipitation grid - UPGRADED 2018 version
Kevin Levey
klevey at customweather.com
Tue Oct 23 13:31:37 EDT 2018
Hi
thanks, but unfortunately, that did not help. I also raised that even thought file supposedly NetCDF, when I opened the file in panoply, it shows it being an HDF5 file, so I changed the DTYPEo
DTYPE hdf5_grid
But get the same result as mentioned below.
This is the output in panoply of this file:
File "nws_precip_1day_conus.nc"
File type: Hierarchical Data Format, version 5
netcdf file:/Users/Stormchaser/Downloads/nws_precip_1day_conus.nc {
dimensions:
y = 881;
x = 1121;
variables:
float percent_of_normal(y=881, x=1121);
:_FillValue = -3.4028235E38f; // float
: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";
:grid_mapping = "crs";
:name = "percent_of_normal";
:coordinates = "x y";
:long_name = "Percent of normal";
:units = "Percent";
:missing_value = -3.4028234663852886E38; // double
: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]]";
:_ChunkSizes = 881, 1121; // int
double y(y=881);
:units = "Meter";
:long_name = "y coordinate of projection";
:standard_name = "projection_y_coordinate";
:_ChunkSizes = 881; // int
float departure_from_normal(y=881, x=1121);
:_FillValue = -3.4028235E38f; // float
: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";
:grid_mapping = "crs";
:name = "departure_from_normal";
:coordinates = "x y";
:long_name = "Departure from normal";
:units = "Inches";
:missing_value = -3.4028234663852886E38; // double
: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]]";
:_ChunkSizes = 881, 1121; // int
float normal(y=881, x=1121);
:_FillValue = -3.4028235E38f; // float
: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";
:grid_mapping = "crs";
:name = "normal";
:coordinates = "x y";
:long_name = "PRISM normals";
:units = "Inches";
:missing_value = -3.4028234663852886E38; // double
: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]]";
:_ChunkSizes = 881, 1121; // int
float observation(y=881, x=1121);
:_FillValue = -3.4028235E38f; // float
: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";
:grid_mapping = "crs";
:name = "observation";
:coordinates = "x y";
:long_name = "Observation";
:units = "Inches";
:missing_value = -3.4028234663852886E38; // double
: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]]";
:_ChunkSizes = 881, 1121; // int
int crs;
:straight_vertical_longitude_from_pole = -105.0; // double
:false_easting = 0.0; // double
:standard_parallel = 60.0; // double
:false_northing = 0.0; // double
:grid_mapping_name = "polar_stereographic";
:latitude_of_projection_origin = 90.0; // double
:earth_radius = 6371200.0; // double
: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_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]]";
: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=1121);
:units = "Meter";
:long_name = "x coordinate of projection";
:standard_name = "projection_x_coordinate";
:_ChunkSizes = 1121; // int
// global attributes:
:creation_time = "20181022T195235";
:data_time = "20181022";
:Conventions = "CF-1.6";
}
> On Oct 22, 2018, at 11:15 PM, yaqiang.wang at gmail.com wrote:
>
> You can try
>
>> PDEF 1121 881 nps 399 1599 -105.0 4.763 <tel:399%201599%20-105.0%204.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 <mailto:yaqiang.wang at gmail.com>
>
> http://www.meteothink.org <http://meteothink.org/>
> https://www.researchgate.net/profile/Yaqiang_Wang5 <https://www.researchgate.net/profile/Yaqiang_Wang5>
> **************************************************
>
> From: Kevin Levey <mailto:klevey at customweather.com>
> Date: 2018-10-23 03:38
> To: GrADS Users Forum <mailto:gradsusr at gradsusr.org>; Jennifer M Adams <mailto:jadams21 at gmu.edu>
> 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 <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 <tel:388%201591> -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 <mailto: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 <tel:388%201591> -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 <tel:401%201601> -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 <mailto:gradsusr at gradsusr.org>
>> http://gradsusr.org/mailman/listinfo/gradsusr <http://gradsusr.org/mailman/listinfo/gradsusr>
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org <mailto:gradsusr at gradsusr.org>
> http://gradsusr.org/mailman/listinfo/gradsusr <http://gradsusr.org/mailman/listinfo/gradsusr>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://gradsusr.org/pipermail/gradsusr/attachments/20181023/84ed2d5c/attachment-0001.html>
More information about the gradsusr
mailing list