<div dir="ltr"><div class="gmail_default" style="font-size:small">Kerwyn,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">My guess is that there's an issue with the coordinates of the "Froz_total" variable, or a lack of them.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">I see that you defined the lat2d/lon2d attributes for u_wind, v_wind, and hgt, but I didn't see anything similar for Froz_total. If you do "printVarSummary" on this variable, does it indicate that it has any kind of lat/lon coordinates at all? This is going to be necessary in order for the "overlay" call to work properly.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">--Mary</div><div class="gmail_default" style="font-size:small"><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Sat, Apr 2, 2016 at 9:08 PM, Kerwyn Texeira <span dir="ltr"><<a href="mailto:ktish86@gmail.com" target="_blank">ktish86@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"><div><div><div><div>Hi ncl-talk,<br><br></div> I trying to overlay both winds and frozen condensate with winds and I'm having a hard time. I'm using gsn_gsm_contour and gsn_csm_vector for the winds and gsn_csm_contour_map for the terrain. Any help will be greatly appreciated.<br><br></div>Thanks,<br><br></div>Error message: warning:ContourPlotDraw: out of range coordinates encountered; standard ORV rendering method may be unreliable;<br> consider setting the resource trGridType to "TriangularMesh" if coordinates contain missing values<br><br></div>Entire script: <br><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_23: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> qc = wrf_user_getvar(a, "QCLOUD", it)<br> qc = qc*1000<br> qr = wrf_user_getvar(a, "QRAIN", it)<br> qr = qr*1000<br> qs = wrf_user_getvar(a, "QSNOW", it)<br> qs = qs*1000<br> qi = wrf_user_getvar(a, "QICE", it)<br> qi = qs*1000<br><br> u_wind = wrf_user_intrp3d(u, p, "h", 600., 0.0, False)<br> v_wind = wrf_user_intrp3d(v, p, "h", 600., 0.0, False)<br> qcl = wrf_user_intrp3d(qc, p, "h", 600., 0.0, False)<br> qrn = wrf_user_intrp3d(qr, p, "h", 600., 0.0, False)<br> qsn = wrf_user_intrp3d(qs, p, "h", 600., 0.0, False)<br> qice = wrf_user_intrp3d(qi, p, "h", 600., 0.0, False)<br><br> Melt_total = qcl + qrn;<br> Froz_total = qsn + qice; <br><br> spd = (u_wind*u_wind + v_wind*v_wind)^(0.5) ;m/s<br> u_wind = u_wind*1.94384449<br> v_wind = v_wind*1.94384449<br><br> u_wind@lon2d = hgt@lon2d<br> u_wind@lat2d = hgt@lat2d<br><br> v_wind@lon2d = hgt@lon2d<br> v_wind@lat2d = hgt@lat2d<br><br> ; spd = spd*1.94384449<br> ; spd@units = "Wind Speed"<br> ; spd@units = "kts"<br> <br> wks_type = "png"<br> wks_type@wkWidth = 2500<br> wks_type@wkHeight = 2500<br> wks = gsn_open_wks(wks_type,"test") ; send graphics to PNG file<br> <br> <br>; gsn_define_colormap(wks,"BlAqGrYeOrRe")<br> gsn_define_colormap(wks,"matlab_jet") <br> <br> res = True<br> res@gsnDraw = False ; do not draw the plot<br> res@gsnFrame = False ; do not advance the frame<br> res@cnLineLabelsOn = False ; do not use line labels<br> res@cnFillOn = True ; color fill<br> res@cnLinesOn = False ; do not draw contour lines<br> <br><br><br> res@tiMainString = ""<br> res@pmTickMarkDisplayMode = "Always"<br> res@mpProjection = "CylindricalEquidistant" ;The default<br> res@mpDataBaseVersion = "MediumRes"<br> res@mpOutlineOn =True<br> res@lbOrientation = "Vertical"<br> res@tiMainOffsetYF = -0.03<br> res@mpFillOn = False<br> res@mpOutlineOn = True ; turn the map outline on<br> res@mpMinLatF = 37.60<br> res@mpMaxLatF = 38.50<br> res@mpMinLonF = -120.0<br> res@mpMaxLonF = -119.10<br> res@gsnLeftString = ""<br> res@gsnCenterString = "Winds (kts) at 600 hpa on Jan 11 at 23:00UTC"<br> res@gsnStringFontHeightF = 0.025<br> res@gsnRightString = ""<br> res@gsnMaximize = True<br> res@mpShapeMode = "FreeAspect" <br> res@lbTitleString = "Terrain (m)"<br> res@lbTitlePosition = "Right"<br> res@lbTitleDirection = "Across"<br> res@lbTitleAngleF = 90.<br> res@lbTitleFontHeightF = 0.020<br> <br><br> ;------wind vectors<br> res2 = True<br> res2@gsnDraw = False<br> res2@gsnFrame = False<br> res2@vcWindBarbLineThicknessF= 3.0<br> res2@vcRefLengthF= 0.018<br> res2@vcRefMagnitudeF= 10<br> res2@vcMinDistanceF = 0.05<br> res2@vcGlyphStyle = "WindBarb"<br> res2@gsnLeftString = ""<br> res2@gsnRightString = ""<br><br>;-------Frozen<br>res3 = True<br>res3@cnLineColor = "Blue"<br>;res3@ContourParameters = (/0.1/)<br>res3@gsnContourLineThicknessesScale = 5.0<br>contour_solid = gsn_csm_contour(wks,Froz_total,res3)<br><br>;------Melted<br>res4 = True<br>res4@cnLineColor = "Black"<br>res4@gsnContourLineThicknessesScale = 5.0<br>contour_liquid = gsn_csm_contour(wks,Melt_total,res4)<br><br><br><br>contour = gsn_csm_contour_map(wks,hgt,res)<br>vector = gsn_csm_vector(wks,u_wind,v_wind,res2)<br>overlay(contour, vector)<br>overlay(contour, contour_solid)<br><br>plres = True<br>plres@gsLineColor = "blue"<br>plres@gsLineThicknessF = 3.0<br><br>lat = (/38.05, 38.05/)<br>lon = (/-120.10, -119.10/)<br><br>plot = gsn_add_polyline(wks,contour,lon,lat,plres)<br><br>draw(contour)<br>frame(wks)<br><br>;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<br>;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<br><br><br>end<br> <br><br></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>