<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
</head>
<body dir="ltr">
<div id="divtagdefaultwrapper" style="font-size:12pt;color:#000000;font-family:Calibri,Helvetica,sans-serif;" dir="ltr">
<p></p>
<div>
<p>Hi James,</p>
<p><br>
</p>
<p>Thanks for the response. I was able to get the data to display, but had to add the xdef information because GrADS wouldn't open the file with out it. My only issue is that I do not believe the xdef and ydef information in the ctl file is accurate even though
the file diplayed ok. Here is the modified ctl file:</p>
<p><br>
</p>
<div>DSET ^nws_precip_1day_20170728_conus.nc<br>
dtype netcdf<br>
options yrev<br>
title amountofprecip<br>
undef -3.40282346638529e+38<br>
pdef 1051 813 nps 388 1591 -105.000000 4.7625<br>
ydef 881 linear 20 0.05<br>
xdef 1121 linear -125 0.05<br>
zdef 1 linear 1 1<br>
tdef 1 linear 00Z28Jul2017 1dy<br>
vars 4<br>
observation=>obs 1 y,x observation<br>
normal=>pcp 1 y,x normal<br>
departure_from_normal=>dpcp 1 y,x departure<br>
percent_of_normal=>ppcp 1 y,x percent<br>
endvars<br>
</div>
<p><br>
</p>
This ctl displays ok, but I'm not certain if the xdef and ydef info is accurate. Can you or anyone else please confirm? Thanks</div>
<br>
<p></p>
<br>
<br>
<div style="color: rgb(0, 0, 0);">
<hr tabindex="-1" style="display:inline-block; width:98%">
<div id="divRplyFwdMsg" dir="ltr"><font style="font-size:11pt" color="#000000" face="Calibri, sans-serif"><b>From:</b> James T. Potemra <jimp@hawaii.edu><br>
<b>Sent:</b> Friday, July 28, 2017 10:35 PM<br>
<b>To:</b> James Spinks<br>
<b>Cc:</b> GrADS Users Forum<br>
<b>Subject:</b> Re: [gradsusr] Precipitation Output Files Changed from NOAA AHPS: Can't read in GrADS anoymore</font>
<div> </div>
</div>
<div>
<p>Hi James,</p>
<p><br>
</p>
<p>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:</p>
<p><br>
</p>
<p>DSET ^nws_precip_1day_20170728_conus.nc<br>
dtype netcdf<br>
options yrev<br>
title amountofprecip<br>
undef -3.40282346638529e+38</p>
<p>pdef 1051 813 nps 388 1591 -105.000000 4.7625<br>
ydef 700 linear 20 0.05<br>
zdef 1 linear 1 1<br>
tdef 1 linear 00Z28Jul2017 1dy<br>
vars 4<br>
observation=>obs 1 y,x observation<br>
normal=>pcp 1 y,x normal<br>
departure_from_normal=>dpcp 1 y,x departure<br>
percent_of_normal=>ppcp 1 y,x percent<br>
endvars<br>
<br>
</p>
<p><br>
</p>
<p>Jim</p>
<p></p>
<blockquote type="cite">netcdf nws_precip_1day_20170728_conus {<br>
dimensions:<br>
y = 881 ;<br>
x = 1121 ;<br>
polar_stereographic = UNLIMITED ; // (0 currently)<br>
variables:<br>
int polar_stereographic(polar_stereographic) ;<br>
polar_stereographic:latitude_of_projection_origin = 90. ;<br>
polar_stereographic:straight_vertical_longitude_from_pole = -105. ;<br>
polar_stereographic:standard_parallel = 60. ;<br>
polar_stereographic:false_northing = 0. ;<br>
polar_stereographic:grid_mapping_name = "polar_stereographic" ;<br>
polar_stereographic:false_easting = 0. ;<br>
float observation(y, x) ;<br>
observation:_FillValue = -3.402823e+38f ;<br>
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" ;<br>
observation:grid_mapping = "polar_stereographic" ;<br>
observation:name = "observation" ;<br>
observation:coordinates = "x y" ;<br>
observation:long_name = "Observation" ;<br>
observation:units = "Inches" ;<br>
observation:missing_value = -3.40282346638529e+38 ;<br>
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]]"
;<br>
float normal(y, x) ;<br>
normal:_FillValue = -3.402823e+38f ;<br>
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" ;<br>
normal:grid_mapping = "polar_stereographic" ;<br>
normal:name = "normal" ;<br>
normal:coordinates = "x y" ;<br>
normal:long_name = "PRISM normals" ;<br>
normal:units = "Inches" ;<br>
normal:missing_value = -3.40282346638529e+38 ;<br>
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]]"
;<br>
float departure_from_normal(y, x) ;<br>
departure_from_normal:_FillValue = -3.402823e+38f ;<br>
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" ;<br>
departure_from_normal:grid_mapping = "polar_stereographic" ;<br>
departure_from_normal:name = "departure_from_normal" ;<br>
departure_from_normal:coordinates = "x y" ;<br>
departure_from_normal:long_name = "Departure from normal" ;<br>
departure_from_normal:units = "Inches" ;<br>
departure_from_normal:missing_value = -3.40282346638529e+38 ;<br>
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]]"
;<br>
float percent_of_normal(y, x) ;<br>
percent_of_normal:_FillValue = -3.402823e+38f ;<br>
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" ;<br>
percent_of_normal:grid_mapping = "polar_stereographic" ;<br>
percent_of_normal:name = "percent_of_normal" ;<br>
percent_of_normal:coordinates = "x y" ;<br>
percent_of_normal:long_name = "Percent of normal" ;<br>
percent_of_normal:units = "Percent" ;<br>
percent_of_normal:missing_value = -3.40282346638529e+38 ;<br>
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]]"
;<br>
double x(x) ;<br>
x:units = "Meter" ;<br>
x:long_name = "x coordinate of projection" ;<br>
x:standard_name = "projection_x_coordinate" ;<br>
double y(y) ;<br>
y:units = "Meter" ;<br>
y:long_name = "y coordinate of projection" ;<br>
y:standard_name = "projection_y_coordinate" ;<br>
<br>
// global attributes:<br>
:creation_time = "20170729T001638" ;<br>
:Conventions = "CF-1.6" ;<br>
}<br>
</blockquote>
<br>
<p></p>
<p><br>
</p>
<br>
<div class="moz-cite-prefix">On 7/28/17 1:52 PM, James Spinks wrote:<br>
</div>
<blockquote type="cite">
<div id="divtagdefaultwrapper" dir="ltr" style="font-size:12pt; color:#000000; font-family:Calibri,Helvetica,sans-serif">
<p>Hi all,</p>
<p><br>
</p>
<p>NOAA AHPS changed their netcdf for format for the precipitation data as stated here:
<a href="http://www.nws.noaa.gov/om/notification/scn17-32ahps_pcpnaaa.htm" class="OWAAutoLink" id="LPlnk442844" previewremoved="true">
http://www.nws.noaa.gov/om/notification/scn17-32ahps_pcpnaaa.htm</a></p>
<p><br>
</p>
<p class="MsoPlainText"><span class="GramE"><span style="font-family:"Courier New"">will</span></span><span style="font-family:"Courier New""> be available in the following formats:</span></p>
<p class="MsoPlainText"><span style=""><b><span style=""> </span>- Climate and Forecast (CF)-Network Common Data Form (<span class="SpellE">netCDF</span>)</b>
</span></p>
<p class="MsoPlainText"><span style=""> </span></p>
<p class="MsoPlainText"><span style="">Precipitation products will no longer be available in the following
</span></p>
<p class="MsoPlainText"><span class="GramE"><span style="font-family:"Courier New"">formats</span></span><span style="font-family:"Courier New"">:</span></p>
<p class="MsoPlainText"><span style=""><b><span style=""> </span>- <span class="SpellE">
<span class="GramE">netCDF</span></span> not adhering to the <span class="SpellE">
netCDF</span> CF metadata convention</b></span></p>
<span style="font-family:"Courier New""><span class="GramE"></span></span>
<p><span style="font-family:"Courier New""><span class="GramE"><br>
</span></span></p>
<p><span style="font-family:"Courier New""><span class="GramE" style="font-family:Calibri,Helvetica,sans-serif">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:</span></span></p>
<p><span style="font-family:"Courier New""><span class="GramE"><br>
</span></span></p>
<p><span style="font-family:"Courier New""><span class="GramE"></span></span></p>
<div><i><span style="font-family:Calibri,Helvetica,sans-serif">DSET ^nws_precip_1day_20170728_conus.nc</span><br>
<br>
<span style="font-family:Calibri,Helvetica,sans-serif">dtype netcdf</span><br>
<span style="font-family:Calibri,Helvetica,sans-serif">title amountofprecip</span><br>
<span style="font-family:Calibri,Helvetica,sans-serif">undef -1</span><br>
<span style="font-family:Calibri,Helvetica,sans-serif">pdef 1051 813 nps 388 1591 -105.000000 4.7625</span><br>
<span style="font-family:Calibri,Helvetica,sans-serif">xdef 1500 linear -135 0.05</span><br>
<span style="font-family:Calibri,Helvetica,sans-serif">ydef 700 linear 20 0.05</span><br>
<span style="font-family:Calibri,Helvetica,sans-serif">zdef 1 linear 1 1</span><br>
<span style="font-family:Calibri,Helvetica,sans-serif">tdef 1 linear 00Z28Jul2017 1dy</span><br>
<span style="font-family:Calibri,Helvetica,sans-serif">vars 1</span><br>
<span style="font-family:Calibri,Helvetica,sans-serif">amountofprecip=>pcp 1 y,x amountofprecip</span><br>
</i><span style="font-family:Calibri,Helvetica,sans-serif"><i>endvars</i></span><br>
</div>
<p><br>
</p>
<p><br>
</p>
<p>Now this ctl does not work nor does ncdump read the new file from <a href="http://water.weather.gov/precip/downloads/" class="OWAAutoLink" id="LPlnk476671" previewremoved="true">
http://water.weather.gov/precip/downloads/</a></p>
<p><br>
</p>
<p>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?</p>
<p><br>
</p>
<p>Thanks,</p>
<p><br>
</p>
<p>James Spinks</p>
<p><br>
</p>
</div>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset> <br>
<pre>_______________________________________________
gradsusr mailing list
<a class="moz-txt-link-abbreviated" href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>
<a class="moz-txt-link-freetext" href="http://gradsusr.org/mailman/listinfo/gradsusr">http://gradsusr.org/mailman/listinfo/gradsusr</a>
</pre>
</blockquote>
<br>
</div>
</div>
</div>
</body>
</html>