<div dir="ltr"><div><div><div><div>Hi Alan Brammer,<br><br></div><div>I appreciate your help in this matter!<br></div><br>I still got a few errors:<br><br></div>1. When i use overlay, I get: warning:wkImageFileName is not a resource in the given object<br>warning:wkImageFormat is not a resource in the given object<br>warning:PlotManagerSetValues: TickMark annotation cannot be added after NhlCreate<br>warning:mpProjection is not a valid resource in test_contour at this time<br>warning:mpMinLonF is not a valid resource in test_contour at this time<br>warning:mpMaxLonF is not a valid resource in test_contour at this time<br>warning:mpMinLatF is not a valid resource in test_contour at this time<br>warning:mpMaxLatF is not a valid resource in test_contour at this time<br>warning:mpShapeMode is not a valid resource in test_contour at this time<br>warning:mpFillOn is not a valid resource in test_contour at this time<br>warning:mpOutlineOn is not a valid resource in test_contour at this time<br>warning:mpDataBaseVersion is not a valid resource in test_contour at this time<br>warning:cnLineOn is not a valid resource in test_contour at this time<br>(0) wrf_map_overlay: Warning: This procedure is obsolete. Consider<br>(0) using wrf_map_overlays instead.<br>warning:ContourPlotDraw: data boundary is out of range<br><br></div>2. I do not understand part of no.5 "you want to transfer the mp resources to a new variable and pass that to wrf_map() instead of True"<br><br></div>The revised script:<br><br>;*************************************************<br>; overlay_1.ncl<br>;<br>; Concepts illustrated:<br>; - Overlaying line contours on filled contours<br>; - Explicitly setting contour levels<br>; - Selecting a different color map<br>;<br>;*************************************************<br>;<br>; These files are loaded by default in NCL V6.2.0 and newer<br> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"<br> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"<br> load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"<br><br>begin<br> a = addfile("./wrfout_d03_2014-01-11_21:00:<a href="http://00.nc">00.nc</a>","r")<br> <br> it = 0<br> hgt = wrf_user_getvar(a, "HGT", it)<br> hgt@lat2d = wrf_user_getvar(a, "XLAT", it)<br> hgt@lon2d = wrf_user_getvar(a, "XLONG", it)<br> u = wrf_user_getvar(a, "ua", it)<br> v = wrf_user_getvar(a, "va", it)<br> p = wrf_user_getvar(a, "pressure", it)<br><br> u_wind = wrf_user_intrp3d(u, p, "h", 650., 0.0, False)<br> v_wind = wrf_user_intrp3d(v, p, "h", 650., 0.0, False)<br> <br><br> spd = (u_wind*u_wind + v_wind*v_wind)^(0.5) ;m/s<br> ; spd = spd*1.94384449<br> ; spd@units = "Wind Speed"<br> ; spd@units = "m/s"<br> <br> <br><br> wks = gsn_open_wks("png","test") ; send graphics to PNG file<br> gsn_define_colormap(wks,"BlAqGrYeOrRe")<br> <br> <br> res = True<br> res@mpFillOn = False<br> res@mpMaxLatF = 38.20 ; specify the plot domain<br> res@mpMinLatF = 37.75 ; <br> res@mpMinLonF = -120.00 ;<br> res@mpMaxLonF = -119.10 ;<br> res@mpOutlineOn = True ; turn the map outline on<br> res@gsnDraw = False ; do not draw the plot<br> res@gsnFrame = False ; do not advance the frame<br>; res@csnLinesOn = True<br>; res@tiMainString = a<br><br> res@mpProjection = "CylindricalEquidistant" ;The default<br> res@pmTickMarkDisplayMode = "Always"<br> res@mpDataBaseVersion = "MediumRes"<br> res@mpOutlineOn =True<br> res@lbOrientation = "Vertical"<br> res@tiMainOffsetYF = -0.03<br> res@mpShapeMode = "FreeAspect"<br> res@vpWidthF = 0.9<br> res@vpHeightF = 0.9<br><br>; res@cnLevelSelectionMode = "ExplicitLevels" ; use explicit levels<br>; res@cnLevels = ispan(215,265,5) ; set the contour levels<br> res@cnLineLabelsOn = False ; do not use line labels<br> res@cnFillOn = True ; color fill<br> res@cnLineOn = False ; do not draw contour lines<br>; res@cnFillPalette = "BlueDarkRed18"<br><br> res@gsnAddCyclic = False<br><br> map = wrf_map(wks, a, True)<br><br>u_wind = u_wind*1.94384449<br>v_wind = v_wind*1.94384449<br><br> contour = gsn_csm_contour(wks,hgt,res) <br> <br> ;------wind vectors<br> res2 = True<br> res2@gsnDraw = True<br> res2@gsnFrame = True<br> res2@FieldTitle = "Winds"<br> res2@NumVectors = 47<br> res2@vcWindBarbLineThicknessF= 2.5<br><br>vector = wrf_vector(a,wks,u_wind,v_wind,res2)<br> wrf_map_overlay(wrf,map,(/contour,vector/),True)<br>;overlay(contour,vector)<br><br><br> <br>end<br> <br> <br><br><br></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Mar 14, 2016 at 9:04 AM, Alan Brammer <span dir="ltr"><<a href="mailto:abrammer@albany.edu" target="_blank">abrammer@albany.edu</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"><div><br></div><div>Not sure where you’re getting those error messages from but there are many problems in the script. Here are the errors I got on running your code on my own file and the step by step to fix them. Make the below edits. If there are still problems send the whole new script to the list so we can see your edits. </div><div><br></div><div><br></div><div><div></div><blockquote type="cite"><div>warning:attsetvalues: the value associated with (tiMainString) does not have an HLU representation</div><div>warning:csnLinesOn is not a valid resource in gsnapp_contour at this time</div><div>(0)<span style="white-space:pre-wrap"> </span>gsn_csm_map_ce: Fatal: The resources mpMinLatF/mpLeftCornerLatF must be less than the resources mpMaxLatF/mpRightCornerF.</div><div>(0)<span style="white-space:pre-wrap"> </span>Execution halted.</div></blockquote></div><div><br></div><div>;====</div><div>1.</div><div>so the top error:</div><div><blockquote type="cite">"warning:attsetvalues: the value associated with (tiMainString) does not have an HLU representation"</blockquote></div><div><br></div><div>Some thing is wrong with tiMainString. " res1@tiMainString = a “ a is a pointer to a file. NCL has no idea how to turn that into a string. </div><div>If you upgrade to 6.3.0 you could do the following.</div><div>res1@tiMainString = getfilepath(a) ;; needs 6.3.0!!!</div><div>or </div><div>res1@tiMainString = “Some string value for a title"</div><div><br></div><div><br></div><div>;====</div><div>2.</div><div><blockquote type="cite">warning:csnLinesOn is not a valid resource in gsnapp_contour at this time</blockquote></div><div>Self explanatory, csnLineOn isn’t a resource. This is a typo. Find that resource and fix it. </div><div><br></div><div><br></div><div>;====</div><div>3.</div><div><div></div><blockquote type="cite"><div>(0)<span style="white-space:pre-wrap"> </span>gsn_csm_map_ce: Fatal: The resources mpMinLatF/mpLeftCornerLatF must be less than the resources mpMaxLatF/mpRightCornerF.</div><div>(0)<span style="white-space:pre-wrap"> </span>Execution halted.</div></blockquote></div><div>Again, does what it says on the tin. minLat must be less than maxLat. Yours are not. Switch your min/max lat resources and you’ll be set. </div><div><br></div><div><br></div><div>;====</div><div>4.</div><div>Once those are fixed, you’ll get the following error.</div><div><div></div><blockquote type="cite"><div>(0)<span style="white-space:pre-wrap"> </span>wrf_map_overlay: Warning: This procedure is obsolete. Consider</div><div>(0)<span style="white-space:pre-wrap"> </span> using wrf_map_overlays instead.</div><div>fatal:NhlAddOverlay: plot class mapPlotClass cannot be overlay plot member</div><div>warning:NhlRemoveOverlay: plot not found in overlay sequence</div></blockquote></div><div><br></div><div>You can’t overlay 2 maps on top of each other. wrf_map() and gsn_csm_contour_map() so just remove the _map part and make a plain contour plot. </div><div><br></div><div><br></div><div>;====</div><div>5.</div><div>Lastly, you want to transfer the mp resources to a new variable and pass that to wrf_map() instead of True. </div><div>and change the order in the overlay function … map, (/contour, vector/)...</div><div><br></div><div><br></div><div><br></div><div><br></div><br><div><blockquote type="cite"><div><div class="h5"><div>On 13 Mar 2016, at 18:12, Kerwyn Texeira <<a href="mailto:ktish86@gmail.com" target="_blank">ktish86@gmail.com</a>> wrote:</div><br></div></div><div><div><div class="h5"><div dir="ltr"><div>Hi All,<br><br>For a month now I have been having trouble overlay windbards on zoomed terrain with specified lat and lon in ncl.<span> I do not want the entire domain plotted with winds. If so, I would have used the multiOverlay.ncl script which works perfectly for many overlay for the entire domain. No images are being plotted with the script below. I'm not getting any fatal errors. I'm using NCL version 6.2.1. All I'm getting is the following; warning:wkImageFileName is not a resource in the given object<br>warning:wkImageFormat is not a resource in the given object.<br><br> </span>I would like some help with this script please.<span> </span>I looked online for a script to help with this (overlay on zoomed in terrain with specified lat and lon) but I’m having trouble finding it.<span>
</span><span>
</span>Any advice will be greatly appreciated. Thanks, Kerwyn<br><br></div>Script:<br><br>;*************************************************<br>; overlay_1.ncl<br>;<br>; Concepts illustrated:<br>; - Overlaying line contours on filled contours<br>; - Explicitly setting contour levels<br>; - Selecting a different color map<br>;<br>;*************************************************<br>;<br>; These files are loaded by default in NCL V6.2.0 and newer<br> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"<br> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"<br> load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"<br><br>begin<br> a = addfile("./wrfout_d03_2014-01-11_21:00:<a href="http://00.nc/" target="_blank">00.nc</a>","r")<br> <br> it = 0<br> hgt = wrf_user_getvar(a, "HGT", it)<br> hgt@lat2d = wrf_user_getvar(a, "XLAT", it)<br> hgt@lon2d = wrf_user_getvar(a, "XLONG", it)<br> u = wrf_user_getvar(a, "ua", it)<br> v = wrf_user_getvar(a, "va", it)<br> p = wrf_user_getvar(a, "pressure", it)<br><br> u_wind = wrf_user_intrp3d(u, p, "h", 650., 0.0, False)<br> v_wind = wrf_user_intrp3d(v, p, "h", 650., 0.0, False)<br> <br><br>; spd = (u_wind*u_wind + v_wind*v_wind)^(0.5) ;m/s<br>; spd@units = "Wind Speed"<br>; spd@units = "m/s"<br><br><br> wks = gsn_open_wks("png","test") ; send graphics to PNG file<br> gsn_define_colormap(wks,"BlAqGrYeOrRe")<br> <br> <br> res1 = True<br> res1@mpFillOn = False<br> res1@mpMaxLatF = 37.75 ; specify the plot domain<br> res1@mpMinLatF = 38.20 ; <br> res1@mpMinLonF = -120.00 ;<br> res1@mpMaxLonF = -119.10 ;<br> res1@mpOutlineOn = True ; turn the map outline on<br> res1@gsnDraw = False ; do not draw the plot<br> res1@gsnFrame = False ; do not advance the frame<br> res1@csnLinesOn = True<br> res1@tiMainString = a<br><br> res1@mpProjection = "CylindricalEquidistant" ;The default<br> res1@pmTickMarkDisplayMode = "Always"<br> res1@mpDataBaseVersion = "MediumRes"<br> res1@mpOutlineOn =True<br> res1@lbOrientation = "Vertical"<br> res1@tiMainOffsetYF = -0.03<br> res1@mpShapeMode = "FreeAspect"<br> res1@vpWidthF = 0.9<br> res1@vpHeightF = 0.9<br><br>; res@cnLevelSelectionMode = "ExplicitLevels" ; use explicit levels<br>; res@cnLevels = ispan(215,265,5) ; set the contour levels<br> res1@cnLineLabelsOn = False ; do not use line labels<br> res1@cnFillOn = True ; color fill<br> res1@cnLinesOn = False ; do not draw contour lines<br>; res@cnFillPalette = "BlueDarkRed18"<br><br> res1@gsnAddCyclic = False<br><br> map = wrf_map(wks, a, True)<br><br> contour = gsn_csm_contour_map(wks,hgt,res1) <br> <br> ;------wind vectors<br> res2 = True<br> res2@FieldTitle = "Winds"<br> res2@NumVectors = 47<br> res2@vcWindBarbLineThicknessF= 2.5<br> vector = wrf_vector(a, wks, u_wind, v_wind, res2)<br><br> wrf_map_overlay(wks,map,(/vector, contour/),True)<br> <br>end<br></div></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" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br></div></blockquote></div><br></div></blockquote></div><br></div>