<div dir="ltr"><div class="gmail_default" style="font-size:small">Hi Brad,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">There are a few problems with the script in general.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">First, as Will pointed out, you are trying to set map resources when you create a contour plot. This is why you are getting errors.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">But, I do see that later you are calling wrf_map_overlays, which is where you want to apply the map resources. However, with the wrf_xxxx plotting scripts, you have to be careful about zooming in on the data when it's being done over a map. The wrf_map_overlays procedure assumes that you are using the map projection that's provided on the WRF output file, and hence it doesn't really know anything about the lat/lon coordinates of your data.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">If you then try to zoom in on your WRF data without telling NCL what part of your data actually corresponds to the area you are zooming in on, then your results are likely going to be incorrect.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">wrf_map_overlays recognizes a special "ZoomIn" resource that allows you to do map zooming. Please see the "olr.ncl" example at:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style=""><a href="http://www2.mmm.ucar.edu/wrf/OnLineTutorial/Graphics/NCL/Examples/SPECIAL/olr.htm">http://www2.mmm.ucar.edu/wrf/OnLineTutorial/Graphics/NCL/Examples/SPECIAL/olr.htm</a><br></div><div class="gmail_default" style=""><br></div><div class="gmail_default" style="">The scripts on the "WRF debug" page also show how to use ZoomIn:</div><div class="gmail_default" style=""><br></div><div class="gmail_default" style=""><a href="http://www.ncl.ucar.edu/Applications/wrfdebug.shtml">http://www.ncl.ucar.edu/Applications/wrfdebug.shtml</a><br></div><div class="gmail_default" style=""><br></div><div class="gmail_default" style="">If you don't want to use the wrf_contour/wrf_map_overlays plotting scripts, then you can plot your WRF data with gsn_csm_contour_map. We have some examples of this, and also show how to zoom in on your data. See wrf_gsn_3.ncl at:</div><div class="gmail_default" style=""><br></div><div class="gmail_default" style=""><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" style=""><br></div><div class="gmail_default" style="">--Mary</div><div class="gmail_default" style=""><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Apr 13, 2016 at 9:47 PM, Will Hobbs <span dir="ltr"><<a href="mailto:will.hobbs@utas.edu.au" target="_blank">will.hobbs@utas.edu.au</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div style="word-wrap:break-word;color:rgb(0,0,0);font-size:14px;font-family:Calibri,sans-serif">
<div>Brad</div>
<div><br>
</div>
<div>I'm not a WRF user so any suggestion here is just a guess. However, have you tried setting the resources trYMinF, trYMaxF etc instead of mpMaxLatF etc.</div>
<div><br>
</div>
<div>I would guess that wrf_contour() is not a map plotting function, so it wouldn't recognise the mp resources. As long as the array that you're plotting has lat/lon coordinate dimensions the bounds set by the 'tr' resources should work fine.</div>
<div><br>
</div>
<div>Hope that helps</div>
<div><br>
</div>
<div>Will</div>
<div><br>
</div>
<div><br>
</div>
<span>
<div style="font-family:Calibri;font-size:11pt;text-align:left;color:black;BORDER-BOTTOM:medium none;BORDER-LEFT:medium none;PADDING-BOTTOM:0in;PADDING-LEFT:0in;PADDING-RIGHT:0in;BORDER-TOP:#b5c4df 1pt solid;BORDER-RIGHT:medium none;PADDING-TOP:3pt">
<span style="font-weight:bold">From: </span><Illston>, "Bradley G." <<a href="mailto:illston@ou.edu" target="_blank">illston@ou.edu</a>><br>
<span style="font-weight:bold">Date: </span>Thursday, 14 April 2016 1:16 PM<br>
<span style="font-weight:bold">To: </span>"<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>" <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>><br>
<span style="font-weight:bold">Subject: </span>[ncl-talk] Zooming Functions Clarification<br>
</div><div><div class="h5">
<div><br>
</div>
<div>
<div lang="EN-US" link="#0563C1" vlink="#954F72">
<div>
<p class="MsoNormal">NCL Users,<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">I am one to try and figure my problems on their own, but I am a bit stumped after working on this for a few days. I am trying to zoom my map in to a set of bounding latitudes and longitudes in my domain to no avail. I kept trying various
iterations utilizing ZoomIn or mpMinLatF/mpMinLonF/etc. I am not sure when you should use one and when you would use the other. Below is my latest script that works, but not zoomed in. When I run it, I get the following warnings:<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">warning:mpLimitMode is not a valid resource in /tair_24hr_2011_2008091502_contour at this time<u></u><u></u></p>
<p class="MsoNormal">warning:mpMinLonF is not a valid resource in /tair_24hr_2011_2008091502_contour at this time<u></u><u></u></p>
<p class="MsoNormal">warning:mpMaxLonF is not a valid resource in /tair_24hr_2011_2008091502_contour at this time<u></u><u></u></p>
<p class="MsoNormal">warning:mpMinLatF is not a valid resource in /tair_24hr_2011_2008091502_contour at this time<u></u><u></u></p>
<p class="MsoNormal">warning:mpMaxLatF is not a valid resource in /tair_24hr_2011_2008091502_contour at this time<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Am I mixing various gsn functions that is causing a confusion? Is there a preferred one for plotting basic WRF output (in NetCDF format) to a PNG file? Any thoughts? Thank you so much.<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Cheers,<u></u><u></u></p>
<p class="MsoNormal">Brad Illston<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"<u></u><u></u></p>
<p class="MsoNormal">load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"<u></u><u></u></p>
<p class="MsoNormal">load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">begin<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">;---Define colormap. First is background color, second is foreground color<u></u><u></u></p>
<p class="MsoNormal"> colors = (/(/255, 255, 255/),(/0, 0, 0/),(/255, 255, 255/),(/088, 164, 074/),(/047,124,062/),(/055, 072, 114/),(/095,129,183/),(/094,165,212/),(/116,196,205/),(/135,099,162/),(/128,076,144/)/) / 255.
<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">;---Open WRF output file<u></u><u></u></p>
<p class="MsoNormal"> a = addfile(wrf_filename,"r")<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"> print(" NetCDF file = "+wrf_filename)<u></u><u></u></p>
<p class="MsoNormal"> print(" Output image = "+output_png)<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">;---Read temperature at first time step<u></u><u></u></p>
<p class="MsoNormal"> tc = wrf_user_getvar(a,"T2",0)<u></u><u></u></p>
<p class="MsoNormal"> tc@long_name = "2m Air Temperature"<u></u><u></u></p>
<p class="MsoNormal"> tc@units = "C"<u></u><u></u></p>
<p class="MsoNormal"> tc=tc-273.15<u></u><u></u></p>
<p class="MsoNormal"><u></u><u></u></p>
<p class="MsoNormal"> wks = gsn_open_wks("png",output_png)<u></u><u></u></p>
<p class="MsoNormal"> gsn_define_colormap(wks,"WhiteBlueGreenYellowRed") ; choose colormap<u></u><u></u></p>
<p class="MsoNormal">; gsn_define_colormap(wks,colors)<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">;---Resources for filled contour plot<u></u><u></u></p>
<p class="MsoNormal"> res = True<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"> res@cnFillOn = True <u></u><u></u></p>
<p class="MsoNormal"> res@cnLinesOn = False<u></u><u></u></p>
<p class="MsoNormal"> res@gsnMaximize = True<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"> res@tiMainString = run_length + " Run (" + year + ") - " + run_date + " " + end_time + " UTC"<u></u><u></u></p>
<p class="MsoNormal"> res@Footer = False<u></u><u></u></p>
<p class="MsoNormal"> res@FieldTitle = "2m Air Temperature"<u></u><u></u></p>
<p class="MsoNormal"> res@lbTitleOn = True<u></u><u></u></p>
<p class="MsoNormal"> res@lbTitleString = ""<u></u><u></u></p>
<p class="MsoNormal"> res@NoHeaderFooter = True ; Switch headers and footers off<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"> res@cnLevelSelectionMode = "ExplicitLevels"<u></u><u></u></p>
<p class="MsoNormal"> res@cnLevels = (/10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25/)<u></u><u></u></p>
<p class="MsoNormal">; res@ContourParameters = (/10,25,1/) ; Contour levels<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"> res@mpLimitMode = "LatLon"<u></u><u></u></p>
<p class="MsoNormal"> res@mpMinLatF = 35<u></u><u></u></p>
<p class="MsoNormal"> res@mpMinLonF = -98<u></u><u></u></p>
<p class="MsoNormal"> res@mpMaxLatF = 36<u></u><u></u></p>
<p class="MsoNormal"> res@mpMaxLonF = -97<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"> contour = wrf_contour(a,wks,tc,res) ; Create contour plot<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"> pltres = True ; Basic overlay plot options<u></u><u></u></p>
<p class="MsoNormal"> pltres@PanelPlot = True ; Tells wrf_map_overlays not to remove overlays<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"> mpres = True ; Set map options<u></u><u></u></p>
<p class="MsoNormal"> mpres@mpOutlineOn = True ; Turn on map outlines<u></u><u></u></p>
<p class="MsoNormal"> mpres@mpFillOn = False ; Turn off map fill<u></u><u></u></p>
<p class="MsoNormal"> mpres@mpDataBaseVersion = "Ncarg4_1"<u></u><u></u></p>
<p class="MsoNormal"> mpres@mpDataSetName = "Earth..2" ; For counties<u></u><u></u></p>
<p class="MsoNormal"> mpres@mpOutlineBoundarySets = "GeophysicalAndUSStates"<u></u><u></u></p>
<p class="MsoNormal"> mpres@mpOutlineSpecifiers = (/"Land", "Oklahoma:counties"/)<u></u><u></u></p>
<p class="MsoNormal"> mpres@tiMainString = run_length + " Run (" + year + ") - " + run_date + " " + end_time + " UTC"<u></u><u></u></p>
<p class="MsoNormal"> mpres@gsnMaximize = True<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">;---Create the contours over the WRF map (nothing will be drawn yet).<u></u><u></u></p>
<p class="MsoNormal"> plot = wrf_map_overlays(a,wks,contour,pltres,mpres)<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">end<u></u><u></u></p>
</div>
</div>
</div>
</div></div></span>
<p style="font-size:10pt;line-height:10pt;font-family:Calibri,sans-serif"><br>
<br>
University of Tasmania Electronic Communications Policy (December, 2014). <br>
This email is confidential, and is for the intended recipient only. Access, disclosure, copying, distribution, or reliance on any of it by anyone outside the intended recipient organisation is prohibited and may be a criminal offence. Please delete if obtained
in error and email confirmation to the sender. The views expressed in this email are not necessarily the views of the University of Tasmania, unless clearly intended otherwise.
</p>
</div>
<br>_______________________________________________<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/mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>