<html><head><meta http-equiv="Content-Type" content="text/html; charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">MON 22OCT18: 1230PDT<div class=""><br class=""></div><div class="">Hi Jennifer,</div><div class=""><br class=""></div><div class="">I was wondering if you could help me here. The HRAP data was updated on 27th JUN this year. </div><div class=""><br class=""></div><div class="">The CTL file that previously you helped Ryan with (below) no longer works. </div><div class=""><br class=""></div><div class="">According to information from the HRAP page:</div><div class=""><br class=""></div><div class=""><a href="https://water.weather.gov/precip/download.php" class="">https://water.weather.gov/precip/download.php</a></div><div class=""><br class=""></div><div class=""><br class=""></div><div class=""><table xmlns="http://www.w3.org/1999/xhtml" cellspacing="0" cellpadding="0" dir="ltr" border="1" style="color: rgb(0, 0, 0); table-layout: fixed; font-size: 10pt; font-family: arial, sans, sans-serif; width: 0px; border-collapse: collapse; border: none;" class=""><colgroup class=""><col width="157" class=""></colgroup><tbody class=""><tr style="height: 21px;" class=""><td data-sheets-value="{"1":2,"2":"The new QPE NetCDFs generated from the NCEP Stage IV data are CF-compliant multi-variable NetCDFs. The variables they contain are:"}" style="border-width: 1px; border-style: solid; border-color: rgb(204, 204, 204) transparent rgb(204, 204, 204) rgb(204, 204, 204); overflow: visible; padding: 2px 0px; vertical-align: bottom; font-family: Arial;" class=""><div style="white-space: nowrap; overflow: hidden; position: relative; width: 894px; left: 3px;" class=""><div style="float: left;" class="">The new QPE NetCDFs generated from the NCEP Stage IV data are CF-compliant multi-variable NetCDFs. The variables they contain are:</div></div></td></tr><tr style="height: 21px;" class=""><td style="border: 1px solid rgb(204, 204, 204); overflow: hidden; padding: 2px 3px; vertical-align: bottom;" class=""></td></tr><tr style="height: 21px;" class=""><td data-sheets-value="{"1":2,"2":"Band 1 - Observation- Last 24 hours of QPE spanning 12Z to 12Z in inches"}" style="border-width: 1px; border-style: solid; border-color: rgb(204, 204, 204) transparent rgb(204, 204, 204) rgb(204, 204, 204); overflow: visible; padding: 2px 0px; vertical-align: bottom; font-family: Arial; font-style: italic;" class=""><div style="white-space: nowrap; overflow: hidden; position: relative; width: 635px; left: 3px;" class=""><div style="float: left;" class="">Band 1 - Observation- Last 24 hours of QPE spanning 12Z to 12Z in inches</div></div></td></tr><tr style="height: 21px;" class=""><td data-sheets-value="{"1":2,"2":"Band 2 - PRISM normals- PRISM normals in inches (see \"Normal Precipitation\" section on the About page)"}" style="border-width: 1px; border-style: solid; border-color: rgb(204, 204, 204) transparent rgb(204, 204, 204) rgb(204, 204, 204); overflow: visible; padding: 2px 0px; vertical-align: bottom; font-family: Arial; font-style: italic; text-decoration: underline;" class=""><div style="white-space: nowrap; overflow: hidden; position: relative; width: 635px; left: 3px;" class=""><div style="float: left;" class=""><a class="in-cell-link" href="https://water.weather.gov/precip/about.php" target="_blank">Band 2 - PRISM normals- PRISM normals in inches (see "Normal Precipitation" section on the About page)</a></div></div></td></tr><tr style="height: 21px;" class=""><td data-sheets-value="{"1":2,"2":"Band 3 - Departure from normal- The departure from normal in inches"}" style="border-width: 1px; border-style: solid; border-color: rgb(204, 204, 204) transparent rgb(204, 204, 204) rgb(204, 204, 204); overflow: visible; padding: 2px 0px; vertical-align: bottom; font-family: Arial; font-style: italic;" class=""><div style="white-space: nowrap; overflow: hidden; position: relative; width: 635px; left: 3px;" class=""><div style="float: left;" class="">Band 3 - Departure from normal- The departure from normal in inches</div></div></td></tr><tr style="height: 21px;" class=""><td data-sheets-value="{"1":2,"2":"Band 4 - Percent of normal- The percent of normal"}" style="border-width: 1px; border-style: solid; border-color: rgb(204, 204, 204) transparent rgb(204, 204, 204) rgb(204, 204, 204); overflow: visible; padding: 2px 0px; vertical-align: bottom; font-family: Arial; font-style: italic;" class=""><div style="white-space: nowrap; overflow: hidden; position: relative; width: 635px; left: 3px;" class=""><div style="float: left;" class="">Band 4 - Percent of normal- The percent of normal</div></div></td></tr></tbody></table><div class=""><br class=""></div></div><div class=""><br class=""></div><div class="">Here is output from today’s 1-day conus file:</div><div class=""><br class=""></div><div class=""><div class="">ncdump -h /mnt/modeldata/domModel/HRAP/rainfall/2018/10/22/nws_precip_1day_conus.nc</div><div class="">netcdf nws_precip_1day_conus {</div><div class="">dimensions:</div><div class="">        y = 881 ;</div><div class="">        x = 1121 ;</div><div class="">variables:</div><div class="">        int crs ;</div><div class="">                crs:straight_vertical_longitude_from_pole = -105. ;</div><div class="">                crs:false_easting = 0. ;</div><div class="">                crs:standard_parallel = 60. ;</div><div class="">                crs:false_northing = 0. ;</div><div class="">                crs:grid_mapping_name = "polar_stereographic" ;</div><div class="">                crs:latitude_of_projection_origin = 90. ;</div><div class="">                crs:earth_radius = 6371200. ;</div><div class="">                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" ;</div><div class="">                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]]" ;</div><div class="">                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]]" ;</div><div class="">        float observation(y, x) ;</div><div class="">                observation:_FillValue = -3.402823e+38f ;</div><div class="">                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" ;</div><div class="">                observation:grid_mapping = "crs" ;</div><div class="">                observation:name = "observation" ;</div><div class="">                observation:coordinates = "x y" ;</div><div class="">                observation:long_name = "Observation" ;</div><div class="">                observation:units = "Inches" ;</div><div class="">                observation:missing_value = -3.40282346638529e+38 ;</div><div class="">                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]]" ;</div><div class="">        float normal(y, x) ;</div><div class="">                normal:_FillValue = -3.402823e+38f ;</div><div class="">                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" ;</div><div class="">                normal:grid_mapping = "crs" ;</div><div class="">                normal:name = "normal" ;</div><div class="">                normal:coordinates = "x y" ;</div><div class="">                normal:long_name = "PRISM normals" ;</div><div class="">                normal:units = "Inches" ;</div><div class="">                normal:missing_value = -3.40282346638529e+38 ;</div><div class="">                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]]" ;</div><div class="">        float departure_from_normal(y, x) ;</div><div class="">                departure_from_normal:_FillValue = -3.402823e+38f ;</div><div class="">                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" ;</div><div class="">                departure_from_normal:grid_mapping = "crs" ;</div><div class="">                departure_from_normal:name = "departure_from_normal" ;</div><div class="">                departure_from_normal:coordinates = "x y" ;</div><div class="">                departure_from_normal:long_name = "Departure from normal" ;</div><div class="">                departure_from_normal:units = "Inches" ;</div><div class="">                departure_from_normal:missing_value = -3.40282346638529e+38 ;</div><div class="">                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]]" ;</div><div class="">        float percent_of_normal(y, x) ;</div><div class="">                percent_of_normal:_FillValue = -3.402823e+38f ;</div><div class="">                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" ;</div><div class="">                percent_of_normal:grid_mapping = "crs" ;</div><div class="">                percent_of_normal:name = "percent_of_normal" ;</div><div class="">                percent_of_normal:coordinates = "x y" ;</div><div class="">                percent_of_normal:long_name = "Percent of normal" ;</div><div class="">                percent_of_normal:units = "Percent" ;</div><div class="">                percent_of_normal:missing_value = -3.40282346638529e+38 ;</div><div class="">                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]]" ;</div><div class="">        double x(x) ;</div><div class="">                x:units = "Meter" ;</div><div class="">                x:long_name = "x coordinate of projection" ;</div><div class="">                x:standard_name = "projection_x_coordinate" ;</div><div class="">        double y(y) ;</div><div class="">                y:units = "Meter" ;</div><div class="">                y:long_name = "y coordinate of projection" ;</div><div class="">                y:standard_name = "projection_y_coordinate" ;</div><div class=""><br class=""></div><div class="">// global attributes:</div><div class="">                :creation_time = "20181022T185215" ;</div><div class="">                :data_time = "20181022" ;</div><div class="">                :Conventions = "CF-1.6" ;</div><div class="">}</div></div><div class=""><br class=""></div><div class=""><br class=""></div><div class="">I changed </div><div class=""><br class=""></div><div class=""><blockquote type="cite" class=""><div class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><div class=""><div class=""><div class=""><blockquote type="cite" class=""><div dir="ltr" class=""><div class="" style="font-size: 12pt; background-color: rgb(255, 255, 255); font-family: Calibri, Arial, Helvetica, sans-serif; position: static; z-index: auto;"><div class=""><div class=""><div class="">amountofprecip=>pcp  1 y,x amountofprecip</div></div></div></div></div></blockquote></div></div></div></div></blockquote><div class=""><br class=""></div>to </div><div class=""><br class=""></div><div class="">observation=>precip  1 y,x Precip Observation</div><div class=""><br class=""></div><div class="">Grads runs and produces this output:</div><div class=""><br class=""></div><div class=""><div class="">ga-> open /mnt/modeldata/domModel/HRAP/rainfall/grads_ctl/hrap_conus.ctl</div><div class="">Scanning description file:  /mnt/modeldata/domModel/HRAP/rainfall/grads_ctl/hrap_conus.ctl</div><div class="">Data file /mnt/modeldata/domModel/HRAP/rainfall/2018/10/22/nws_precip_1day_conus.nc is open as file 1</div><div class="">LON set to -135 -60.05 </div><div class="">LAT set to 20 54.95 </div><div class="">LEV set to 1 1 </div><div class="">Time values set: 2018:10:22:0 2018:10:22:0 </div><div class="">E set to 1 1 </div><div class="">Notice: Implied interpolation for file /mnt/modeldata/domModel/HRAP/rainfall/grads_ctl/hrap_conus.ctl</div><div class=""> Interpolation will be performed on any data displayed from this file</div><div class="">ga-> q file 1</div><div class="">File 1 : Gridded 4km x 4km HRAP 24-Hour Observed Rainfall ending 8AM ET</div><div class="">  Descriptor: /mnt/modeldata/domModel/HRAP/rainfall/grads_ctl/hrap_conus.ctl</div><div class="">  Binary: /mnt/modeldata/domModel/HRAP/rainfall/2018/10/22/nws_precip_1day_conus.nc</div><div class="">  Type = Gridded</div><div class="">  Xsize = 1500  Ysize = 700  Zsize = 1  Tsize = 1  Esize = 1</div><div class="">  Number of Variables = 1</div><div class="">     precip  1  y,x  Precip Observation</div><div class="">ga-> set gxout shaded</div><div class="">ga-> d precip</div><div class="">Notice:  Automatic Grid Interpolation Taking Place</div><div class="">Contouring: -3.3e+38 to 0 interval 3e+37 </div></div><div class=""><br class=""></div><div class="">but nothing is plotted - the map shows no data, only USA outline.</div><div class=""><br class=""></div><div class="">Here is my control file:</div><div class=""><br class=""></div><div class=""><div class="">DSET /mnt/modeldata/domModel/HRAP/rainfall/2018/10/22/nws_precip_1day_conus.nc</div><div class="">DTYPE netcdf </div><div class="">TITLE Gridded 4km x 4km HRAP 24-Hour Observed Rainfall ending 8AM ET</div><div class="">undef -1</div><div class="">PDEF 1051 813 nps 388 1591 -105.000000 4.7625</div><div class="">XDEF 1500 linear -135 0.05</div><div class="">YDEF 700 linear 20 0.05</div><div class="">ZDEF 1 linear 1 1</div><div class="">TDEF 1 linear 00:00Z22oct2018 1dy  </div><div class="">VARS 1</div><div class="">observation=>precip  1 y,x Precip Observation</div><div class="">ENDVARS</div><div class=""><br class=""></div><div class="">Any help here would be much appreciated. </div><div class=""><br class=""></div><div class="">Cheers</div><div class="">Kevin</div><div class=""><br class=""></div><div><br class=""><blockquote type="cite" class=""><div class="">On Apr 18, 2014, at 10:43 AM, Jennifer Adams <<a href="mailto:jma@COLA.IGES.ORG" class="">jma@COLA.IGES.ORG</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; " class="">Hi, Ryan -- <div class="">I subtracted hrap.xor from 401 and hrap.yor from 1601 and now it seems to fit with this PDEF:</div><div class=""><div class="">pdef 1051 813 nps 388 1591 -105.000000 4.7625</div><div class="">xdef 1500 linear -135 0.05</div><div class="">ydef 700 linear 20 0.05</div><div class=""><br class=""></div><div class="">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:</div><div class=""><br class=""></div><div class=""><span id="cid:72EE7728-60EE-4D3C-B81F-20ADFAC2C4C9@SME"><Screen Shot 2014-04-18 at 1.42.01 PM.png></span></div><div class=""><br class=""></div><div class=""><br class=""></div><div class="">--Jennifer</div><div class=""><br class=""></div><div class=""><div class=""><div class="">On Apr 9, 2014, at 5:14 PM, Maue, Ryan wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite" class=""><div dir="ltr" class=""><div style="font-size: 12pt; background-color: rgb(255, 255, 255); font-family: Calibri, Arial, Helvetica, sans-serif; position: static; z-index: auto;" class=""><div style="margin-top: 0px; margin-bottom: 0px; " class="">Greetings,<br class=""></div><div style="margin-top: 0px; margin-bottom: 0px; " class=""><br class=""></div><div style="margin-top: 0px; margin-bottom: 0px; " class="">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.<br class=""></div><div style="margin-top: 0px; margin-bottom: 0px; " class=""><br class=""></div><div style="margin-top: 0px; margin-bottom: 0px; " class="">Any help or guidance would be great!<br class=""></div><div style="margin-top: 0px; margin-bottom: 0px; " class=""><br class=""></div><div style="margin-top: 0px; margin-bottom: 0px; " class="">I've uploaded an example here as well as the ascii output from the conversion program:<br class=""></div><div style="margin-top: 0px; margin-bottom: 0px; " class=""><a href="http://models.weatherbell.com/news/nws_precip_conus_20140101.nc" id="lnk448130" class="">http://models.weatherbell.com/news/nws_precip_conus_20140101.nc</a><br class=""></div><div style="margin-top: 0px; margin-bottom: 0px; " class=""><a href="http://models.weatherbell.com/news/20140101.txt" id="lnk685508" class="">http://models.weatherbell.com/news/20140101.txt</a><br class=""></div><div style="margin-top: 0px; margin-bottom: 0px; " class=""><br class=""></div><div style="margin-top: 0px; margin-bottom: 0px; " class=""><span style="font-size: 12pt; " class="">Location of the NWS Precip page:</span><br class=""></div><div style="margin-top: 0px; margin-bottom: 0px; " class=""><a href="http://water.weather.gov/precip/" id="lnk839083" style="font-size: 12pt; " class="">http://water.weather.gov/precip/</a><br class=""></div><div style="margin-top: 0px; margin-bottom: 0px; " class="">​<br class=""></div><div class="">dset ^nws_precip_conus_20140101.nc</div><div class="">dtype netcdf</div><div class="">title amountofprecip</div><div class="">undef -1</div><div class="">pdef 1051 813 nps 401 1601 -105.000000 4.7625</div><div class=""><div class="">xdef 1051 linear -131.402908 .06366019047619047619</div><div class="">ydef 813 linear  20.761557  .04260785837438423645</div><div class=""><div class="">ydef 813 linear  20.749580    0.04262164576</div><div class="">zdef 1 linear 1 1</div><div class="">tdef 1 linear 00:00Z01jan2014 1dy</div><div class="">vars 1</div><div class="">amountofprecip=>pcp  1 y,x amountofprecip</div><div class="">endvars<br class=""></div><div class=""><br class=""></div><div class="">RYAN MAUE<br class=""></div><div class=""><br class=""><br class=""></div><div class=""><br class=""></div><br class=""><br class=""></div><br class=""><br class=""></div><div style="margin-top: 0px; margin-bottom: 0px; " class=""><br class=""></div></div>_______________________________________________<br class="">gradsusr mailing list<br class=""><a href="mailto:gradsusr@gradsusr.org" class="">gradsusr@gradsusr.org</a><br class=""><a href="http://gradsusr.org/mailman/listinfo/gradsusr" class="">http://gradsusr.org/mailman/listinfo/gradsusr</a><br class=""></div></blockquote></div><br class=""><div apple-content-edited="true" class="">
<span class="Apple-style-span" style="font-size: 12px;"><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; " class=""><span class="Apple-style-span" style="border-collapse: separate; border-spacing: 0px; font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-ligatures: normal; font-variant-caps: normal; font-variant-east-asian: normal; font-variant-position: normal; font-weight: normal; letter-spacing: normal; line-height: normal; -webkit-text-decorations-in-effect: none; text-indent: 0px; text-transform: none; orphans: 2; white-space: normal; widows: 2; word-spacing: 0px;"><span class="Apple-style-span" style="border-collapse: separate; border-spacing: 0px; font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-ligatures: normal; font-variant-caps: normal; font-variant-east-asian: normal; font-variant-position: normal; font-weight: normal; letter-spacing: normal; line-height: normal; -webkit-text-decorations-in-effect: none; text-indent: 0px; text-transform: none; orphans: 2; white-space: normal; widows: 2; word-spacing: 0px;"><div class="">--</div><div class="">Jennifer M. Adams<br class="">Center for Ocean-Land-Atmosphere Studies (COLA)<br class="">111 Research Hall, Mail Stop 2B3<br class="">George Mason University<br class="">4400 University Drive<br class="">Fairfax, VA 22030 <br class=""><br class=""></div><div class=""><br class="khtml-block-placeholder"></div><br class="Apple-interchange-newline"></span></span></div></span><br class="Apple-interchange-newline">
</div>
<br class=""></div></div></div>_______________________________________________<br class="">gradsusr mailing list<br class=""><a href="mailto:gradsusr@gradsusr.org" class="">gradsusr@gradsusr.org</a><br class="">http://gradsusr.org/mailman/listinfo/gradsusr<br class=""></div></blockquote></div><br class=""></div></body></html>