<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri",sans-serif;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:#0563C1;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:#954F72;
text-decoration:underline;}
span.EmailStyle17
{mso-style-type:personal-compose;
font-family:"Calibri",sans-serif;
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri",sans-serif;}
@page WordSection1
{size:8.5in 11.0in;
margin:1.0in 1.0in 1.0in 1.0in;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="EN-US" link="#0563C1" vlink="#954F72">
<div class="WordSection1">
<p class="MsoNormal">NCL Users,<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></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:<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">warning:mpLimitMode is not a valid resource in /tair_24hr_2011_2008091502_contour at this time<o:p></o:p></p>
<p class="MsoNormal">warning:mpMinLonF is not a valid resource in /tair_24hr_2011_2008091502_contour at this time<o:p></o:p></p>
<p class="MsoNormal">warning:mpMaxLonF is not a valid resource in /tair_24hr_2011_2008091502_contour at this time<o:p></o:p></p>
<p class="MsoNormal">warning:mpMinLatF is not a valid resource in /tair_24hr_2011_2008091502_contour at this time<o:p></o:p></p>
<p class="MsoNormal">warning:mpMaxLatF is not a valid resource in /tair_24hr_2011_2008091502_contour at this time<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></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.<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Cheers,<o:p></o:p></p>
<p class="MsoNormal">Brad Illston<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"<o:p></o:p></p>
<p class="MsoNormal">load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"<o:p></o:p></p>
<p class="MsoNormal">load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">begin<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">;---Define colormap. First is background color, second is foreground color<o:p></o:p></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.
<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">;---Open WRF output file<o:p></o:p></p>
<p class="MsoNormal"> a = addfile(wrf_filename,"r")<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"> print(" NetCDF file = "+wrf_filename)<o:p></o:p></p>
<p class="MsoNormal"> print(" Output image = "+output_png)<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">;---Read temperature at first time step<o:p></o:p></p>
<p class="MsoNormal"> tc = wrf_user_getvar(a,"T2",0)<o:p></o:p></p>
<p class="MsoNormal"> tc@long_name = "2m Air Temperature"<o:p></o:p></p>
<p class="MsoNormal"> tc@units = "C"<o:p></o:p></p>
<p class="MsoNormal"> tc=tc-273.15<o:p></o:p></p>
<p class="MsoNormal"><o:p></o:p></p>
<p class="MsoNormal"> wks = gsn_open_wks("png",output_png)<o:p></o:p></p>
<p class="MsoNormal"> gsn_define_colormap(wks,"WhiteBlueGreenYellowRed") ; choose colormap<o:p></o:p></p>
<p class="MsoNormal">; gsn_define_colormap(wks,colors)<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">;---Resources for filled contour plot<o:p></o:p></p>
<p class="MsoNormal"> res = True<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"> res@cnFillOn = True <o:p></o:p></p>
<p class="MsoNormal"> res@cnLinesOn = False<o:p></o:p></p>
<p class="MsoNormal"> res@gsnMaximize = True<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"> res@tiMainString = run_length + " Run (" + year + ") - " + run_date + " " + end_time + " UTC"<o:p></o:p></p>
<p class="MsoNormal"> res@Footer = False<o:p></o:p></p>
<p class="MsoNormal"> res@FieldTitle = "2m Air Temperature"<o:p></o:p></p>
<p class="MsoNormal"> res@lbTitleOn = True<o:p></o:p></p>
<p class="MsoNormal"> res@lbTitleString = ""<o:p></o:p></p>
<p class="MsoNormal"> res@NoHeaderFooter = True ; Switch headers and footers off<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"> res@cnLevelSelectionMode = "ExplicitLevels"<o:p></o:p></p>
<p class="MsoNormal"> res@cnLevels = (/10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25/)<o:p></o:p></p>
<p class="MsoNormal">; res@ContourParameters = (/10,25,1/) ; Contour levels<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"> res@mpLimitMode = "LatLon"<o:p></o:p></p>
<p class="MsoNormal"> res@mpMinLatF = 35<o:p></o:p></p>
<p class="MsoNormal"> res@mpMinLonF = -98<o:p></o:p></p>
<p class="MsoNormal"> res@mpMaxLatF = 36<o:p></o:p></p>
<p class="MsoNormal"> res@mpMaxLonF = -97<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"> contour = wrf_contour(a,wks,tc,res) ; Create contour plot<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"> pltres = True ; Basic overlay plot options<o:p></o:p></p>
<p class="MsoNormal"> pltres@PanelPlot = True ; Tells wrf_map_overlays not to remove overlays<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"> mpres = True ; Set map options<o:p></o:p></p>
<p class="MsoNormal"> mpres@mpOutlineOn = True ; Turn on map outlines<o:p></o:p></p>
<p class="MsoNormal"> mpres@mpFillOn = False ; Turn off map fill<o:p></o:p></p>
<p class="MsoNormal"> mpres@mpDataBaseVersion = "Ncarg4_1"<o:p></o:p></p>
<p class="MsoNormal"> mpres@mpDataSetName = "Earth..2" ; For counties<o:p></o:p></p>
<p class="MsoNormal"> mpres@mpOutlineBoundarySets = "GeophysicalAndUSStates"<o:p></o:p></p>
<p class="MsoNormal"> mpres@mpOutlineSpecifiers = (/"Land", "Oklahoma:counties"/)<o:p></o:p></p>
<p class="MsoNormal"> mpres@tiMainString = run_length + " Run (" + year + ") - " + run_date + " " + end_time + " UTC"<o:p></o:p></p>
<p class="MsoNormal"> mpres@gsnMaximize = True<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">;---Create the contours over the WRF map (nothing will be drawn yet).<o:p></o:p></p>
<p class="MsoNormal"> plot = wrf_map_overlays(a,wks,contour,pltres,mpres)<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">end<o:p></o:p></p>
</div>
</body>
</html>