<div dir="ltr">Hi Mary,<div><br></div><div>Thanks for getting back to me. Sorry about the delay on my end.</div><div><br></div><div><font face="monospace, monospace">printVarSummary of the rlat (or rlon) array shows</font></div><div><font face="monospace, monospace"><br></font></div><div><div><font face="monospace, monospace">Variable: rlat</font></div><div><font face="monospace, monospace">Type: float</font></div><div><font face="monospace, monospace">Total Size: 131044 bytes</font></div><div><font face="monospace, monospace"> 32761 values</font></div><div><font face="monospace, monospace">Number of Dimensions: 2</font></div><div><font face="monospace, monospace">Dimensions and sizes:<span class="" style="white-space:pre"> </span>[y0 | 181] x [x0 | 181]</font></div><div><font face="monospace, monospace">Coordinates: </font></div><div><font face="monospace, monospace"> y0: [-120..119.94]</font></div><div><font face="monospace, monospace"> x0: [-120..119.94]</font></div><div><font face="monospace, monospace">Number Of Attributes: 3</font></div><div><font face="monospace, monospace"> _FillValue :<span class="" style="white-space:pre"> </span>41.94297</font></div><div><font face="monospace, monospace"> standard_name :<span class="" style="white-space:pre"> </span>latitude</font></div><div><font face="monospace, monospace"> units :<span class="" style="white-space:pre"> </span>degrees_north</font></div></div><div><br></div><div><br></div><div>and for rd_plane shows</div><div><br></div><div><br></div><div><div><font face="monospace, monospace">Variable: rd_plane</font></div><div><font face="monospace, monospace">Type: float</font></div><div><font face="monospace, monospace">Total Size: 131044 bytes</font></div><div><font face="monospace, monospace"> 32761 values</font></div><div><font face="monospace, monospace">Number of Dimensions: 2</font></div><div><font face="monospace, monospace">Dimensions and sizes:<span class="" style="white-space:pre"> </span>[y0 | 181] x [x0 | 181]</font></div><div><font face="monospace, monospace">Coordinates: </font></div><div><font face="monospace, monospace"> y0: [-120..119.94]</font></div><div><font face="monospace, monospace"> x0: [-120..119.94]</font></div><div><font face="monospace, monospace">Number Of Attributes: 12</font></div><div><font face="monospace, monospace"> z0 :<span class="" style="white-space:pre"> </span> 2</font></div><div><font face="monospace, monospace"> time :<span class="" style="white-space:pre"> </span>1424628563</font></div><div><font face="monospace, monospace"> grid_mapping :<span class="" style="white-space:pre"> </span>grid_mapping_0</font></div><div><font face="monospace, monospace"> coordinates :<span class="" style="white-space:pre"> </span>lon0 lat0</font></div><div><font face="monospace, monospace"> units :<span class="" style="white-space:pre"> </span>dBZ</font></div><div><font face="monospace, monospace"> long_name :<span class="" style="white-space:pre"> </span>radar_reflectivity</font></div><div><font face="monospace, monospace"> standard_name :<span class="" style="white-space:pre"> </span>REF</font></div><div><font face="monospace, monospace"> _FillValue :<span class="" style="white-space:pre"> </span>9.96921e+36</font></div><div><font face="monospace, monospace"> valid_max :<span class="" style="white-space:pre"> </span>25.3757</font></div><div><font face="monospace, monospace"> valid_min :<span class="" style="white-space:pre"> </span>-9.187451</font></div></div><div><br></div><div>printMinMax shows:</div><div><br></div><div><div><font face="monospace, monospace">(0)<span class="" style="white-space:pre"> </span>latitude: min=39.779 max=41.943</font></div><div><font face="monospace, monospace">(0)<span class="" style="white-space:pre"> </span>longitude: min=-74.313 max=-71.4156</font></div></div><div><br></div><div>The values of rlat and rlon are the correct latitude and longitude values, but the trouble seems to be in matching the i,j'th point in rd_plane with the lat,lon contained in the i'th point of rlon and j'th point of rlat, so that the code can plot radar reflectivity in the correct place.</div><div><br></div><div>Thanks,</div><div>Ryan</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Sun, Jul 10, 2016 at 11:53 PM, Mary Haley <span dir="ltr"><<a href="mailto:haley@ucar.edu" target="_blank">haley@ucar.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div class="gmail_default" style="font-size:small">Hi Ryan,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">You can't use wrf_xxxx functions to plot anything but WRF-ARW data. These functions assume that you are using WRF files with names like "wrfout_d01_2008-09-28_00:00:00" which have a bunch of global attributes in them that define the native map projection of the WRF data. This is why you were getting an error about MAP_PROJ being undefined. This is one of the required global attributes.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">If you are trying to plot the data over a map, then you will need to examine your file and your data variable (rd_plane) to determine if it is also on some native map projection, or whether you need to provide lat/lon coordinates for plotting. I then suggest using the gsn_csm_contour or gsn_csm_contour_map functions to generate just contours, or contours over a map.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">Please see this page for some examples on how to plot data over a map using NCL:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default"><a href="http://www.ncl.ucar.edu/Applications/plot_data_on_map.shtml" target="_blank">http://www.ncl.ucar.edu/Applications/plot_data_on_map.shtml</a><br></div><div class="gmail_default"><br></div><div class="gmail_default">You might also want to look at this page:</div><div class="gmail_default"><br></div><div class="gmail_default"><a href="http://www.ncl.ucar.edu/Applications/radar.shtml" target="_blank">http://www.ncl.ucar.edu/Applications/radar.shtml</a><br></div><div class="gmail_default"><br></div><div class="gmail_default"><br></div><div class="gmail_default">--Mary</div><div class="gmail_default"><br></div>
</div><div class="gmail_extra"><br><div class="gmail_quote"><div><div class="h5">On Fri, Jul 8, 2016 at 4:35 PM, Ryan Connelly <span dir="ltr"><<a href="mailto:rconne01@gmail.com" target="_blank">rconne01@gmail.com</a>></span> wrote:<br></div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="h5"><div dir="ltr">Hi,<div><br></div><div>I have Level II radar data in netcdf format (using radx and the steps outlined here: <a href="https://publish.illinois.edu/snesbitt/resources/radar-data-processing-with-radx/" target="_blank">https://publish.illinois.edu/snesbitt/resources/radar-data-processing-with-radx/</a>).</div><div><br></div><div>I want to overlay it on a base map, in the same way I would with WRF data. Here's a highly condensed version of what I have:</div><div><br></div><div><div><font face="monospace, monospace">rf = addfile(diri+fili,"r")</font></div><div><font face="monospace, monospace"> rd = rf->REF(:,:,:,:) ; float REF(time, z0, y0, x0) ; RHO, ZDR</font></div><div><font face="monospace, monospace"><br></font></div><div><span style="font-family:monospace,monospace"> rd_plane = rd(0,0,:,:)</span><br></div><div><div><font face="monospace, monospace"><br></font></div><div><font face="monospace, monospace"> contour = gsn_csm_contour(wks,rd_plane,radres)</font></div><div><font face="monospace, monospace"> plot = wrf_map_overlays(rf,wks,(/contour/),radres,mpres)</font></div></div><div><br></div><div>I was just assuming the netcdf had location data packaged into it, but I get the error:</div><div><br></div><div><div><font face="monospace, monospace">(0)<span style="white-space:pre-wrap"> </span>wrf_map: Error: no MAP_PROJ attribute in input file</font></div><div><font face="monospace, monospace">fatal:overlay: missing value as input, can't continue</font></div></div><div><br></div><div>I tried adding in a wrf file I had laying around and using it for map plotting, and while that got be a basemap, by radar data was obviously not matched geographically (it showed the KOKX RDA over New York City!).</div><div><br></div><div>Thanks,</div><div>Ryan</div><span><font color="#888888"><div><br></div>-- <br><div data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><font size="1">Ryan Connelly</font></div><div dir="ltr"><div dir="ltr" style="font-size:12.8px"><font color="#000000"><font size="1">M.S. Student in Atmospheric Sciences, </font><span style="font-size:x-small">Stony Brook University</span><br></font></div><div dir="ltr" style="font-size:12.8px"><div><font color="#000000"><font size="1">B.S. in Meteorology with Minors in Mathematics and GIS, </font><span style="font-size:x-small">Valparaiso University</span></font></div></div><div><a href="mailto:rconne01@gmail.com" style="font-size:x-small" target="_blank">rconne01@gmail.com</a><br></div><div><font size="1"><a href="mailto:ryan.connelly@stonybrook.edu" target="_blank">ryan.connelly@stonybrook.edu</a></font></div></div></div></div></div></div></div></div></div></div>
</font></span></div></div>
<br></div></div>_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><font size="1">Ryan Connelly</font></div><div dir="ltr"><div dir="ltr" style="font-size:12.8000001907349px"><font color="#000000"><font size="1">M.S. Student in Atmospheric Sciences, </font><span style="font-size:x-small">Stony Brook University</span><br></font></div><div dir="ltr" style="font-size:12.8000001907349px"><div><font color="#000000"><font size="1">B.S. in Meteorology with Minors in Mathematics and GIS, </font><span style="font-size:x-small">Valparaiso University</span></font></div></div><div><a href="mailto:rconne01@gmail.com" style="font-size:x-small" target="_blank">rconne01@gmail.com</a><br></div><div><font size="1"><a href="mailto:ryan.connelly@stonybrook.edu" target="_blank">ryan.connelly@stonybrook.edu</a></font></div></div></div></div></div></div></div></div></div></div>
</div>