<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">The WRF plotting routines want to plot data in the native projection on the file, and won't let you set a different map projection. This is actually why the file handle ("a" in your case) is the first argument to wrf_map_overlays, because it queries the file for the map projection information and sets it.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">It looks like your WRF data may be on an LCC projection, and so that's what it's trying to use.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">If you really want to set your own map projection, then I suggest you use gsn_csm_contour_map. To do this, you will need to read in XLAT and XLONG off the file, and attach them as attributes called "lat2d" and "lon2d".<br></div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">See example #10 at:</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#ex10">http://www.ncl.ucar.edu/Applications/plot_data_on_map.shtml#ex10</a><br></div><div class="gmail_default"><br></div><div class="gmail_default">The first frame shows how to plot data using the native projection on the WRF output file. The second frame, which is the one you want, shows how to read XLAT/XLONG off the file. You want to follow this example, and then set your map resources as you indicated.</div><div class="gmail_default"><br></div><div class="gmail_default">A more simple example is the first one at:</div><div class="gmail_default"><br></div><div class="gmail_default"><a href="http://www.ncl.ucar.edu/Applications/wrfgsn.shtml">http://www.ncl.ucar.edu/Applications/wrfgsn.shtml</a><br></div><div class="gmail_default"><br></div><div class="gmail_default">See the first frame of this example (not the second one, which uses wrf_xxxx to plot).</div><div class="gmail_default"><br></div><div class="gmail_default">It plots the data on a Cylindrical Equidistant map, but you can change this to mercator.</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">On Tue, Jan 17, 2017 at 12:29 PM, Ryan Connelly <span dir="ltr"><<a href="mailto:rconne01@gmail.com" target="_blank">rconne01@gmail.com</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">I'm following this page as a guide: <a href="http://www.ncl.ucar.edu/Applications/wrfmerc.shtml" target="_blank">http://www.ncl.ucar.<wbr>edu/Applications/wrfmerc.shtml</a><div><br></div><div>I've added the following lines to my script:</div><div><div><font face="monospace, monospace"> </font></div><div><font face="monospace, monospace">; Plotting options for limited area Mercator projection plot consistent across domains</font></div><div><font face="monospace, monospace"><br></font></div><div><font face="monospace, monospace"> mpres@tfDoNDCOverlay = True ; set True for native (direct) mapping</font></div><div><font face="monospace, monospace"><br></font></div><div><font face="monospace, monospace"> mpres@mpProjection = "Mercator"</font></div><div><font face="monospace, monospace"> mpres@mpLimitMode = "Corners"</font></div><div><font face="monospace, monospace"> mpres@mpRightCornerLatF = 41.22405</font></div><div><font face="monospace, monospace"> mpres@mpLeftCornerLonF = -72.25338</font></div><div><font face="monospace, monospace"> mpres@mpLeftCornerLatF = 36.38864</font></div><div><font face="monospace, monospace"> mpres@mpRightCornerLonF = -78.62663</font></div><div><br></div><div>And I call the plot like this:</div><div><font face="monospace, monospace"> plot = wrf_map_overlays(a,wks,(/<wbr>contour/),pltres,mpres)</font><br></div><div><br></div><div>I get the correct lat,lon corners, but on an LCC plot, not a Mercator. Please advise.</div><span class="HOEnZb"><font color="#888888"><div><br></div><div>Ryan</div><div><br></div>-- <br><div class="m_-9026195970611737775gmail_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>______________________________<wbr>_________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu">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/<wbr>mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>