<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&#39;m having a hard time. I&#39;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 &quot;TriangularMesh&quot; if coordinates contain missing values<br><br></div>Entire script: <br><br>load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl&quot;<br> load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl&quot;<br> load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl&quot;<br> <br>begin<br>   a = addfile(&quot;./wrfout_d03_2014-01-11_23:00:<a href="http://00.nc">00.nc</a>&quot;,&quot;r&quot;)<br><br>  it = 0<br>  hgt       = wrf_user_getvar(a, &quot;HGT&quot;, it)<br>  hgt@lat2d = wrf_user_getvar(a, &quot;XLAT&quot;, it)<br>  hgt@lon2d = wrf_user_getvar(a, &quot;XLONG&quot;, it)<br>  u         = wrf_user_getvar(a, &quot;ua&quot;, it)<br>  v         = wrf_user_getvar(a, &quot;va&quot;, it)<br>  p         = wrf_user_getvar(a, &quot;pressure&quot;, it)<br>  qc        = wrf_user_getvar(a, &quot;QCLOUD&quot;, it)<br>  qc        = qc*1000<br>  qr        = wrf_user_getvar(a, &quot;QRAIN&quot;, it)<br>  qr        = qr*1000<br>  qs        = wrf_user_getvar(a, &quot;QSNOW&quot;, it)<br>  qs        = qs*1000<br>  qi        = wrf_user_getvar(a, &quot;QICE&quot;, it)<br>  qi        = qs*1000<br><br>  u_wind    = wrf_user_intrp3d(u, p, &quot;h&quot;, 600., 0.0, False)<br>  v_wind    = wrf_user_intrp3d(v, p, &quot;h&quot;, 600., 0.0, False)<br>  qcl       = wrf_user_intrp3d(qc, p, &quot;h&quot;, 600., 0.0, False)<br>  qrn       = wrf_user_intrp3d(qr, p, &quot;h&quot;, 600., 0.0, False)<br>  qsn       = wrf_user_intrp3d(qs, p, &quot;h&quot;, 600., 0.0, False)<br>  qice      = wrf_user_intrp3d(qi, p, &quot;h&quot;, 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 = &quot;Wind Speed&quot;<br> ; spd@units = &quot;kts&quot;<br>  <br>    wks_type = &quot;png&quot;<br>    wks_type@wkWidth = 2500<br>    wks_type@wkHeight = 2500<br>    wks = gsn_open_wks(wks_type,&quot;test&quot;)         ; send graphics to PNG file<br>  <br>  <br>; gsn_define_colormap(wks,&quot;BlAqGrYeOrRe&quot;)<br>  gsn_define_colormap(wks,&quot;matlab_jet&quot;) <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 = &quot;&quot;<br>  res@pmTickMarkDisplayMode  = &quot;Always&quot;<br>  res@mpProjection  = &quot;CylindricalEquidistant&quot;    ;The default<br>  res@mpDataBaseVersion      = &quot;MediumRes&quot;<br>  res@mpOutlineOn            =True<br>  res@lbOrientation          = &quot;Vertical&quot;<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 = &quot;&quot;<br>  res@gsnCenterString = &quot;Winds (kts) at 600 hpa on Jan 11 at 23:00UTC&quot;<br>  res@gsnStringFontHeightF = 0.025<br>  res@gsnRightString = &quot;&quot;<br>  res@gsnMaximize     = True<br>  res@mpShapeMode     = &quot;FreeAspect&quot;  <br>  res@lbTitleString   = &quot;Terrain (m)&quot;<br>  res@lbTitlePosition = &quot;Right&quot;<br>  res@lbTitleDirection = &quot;Across&quot;<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 = &quot;WindBarb&quot;<br>  res2@gsnLeftString = &quot;&quot;<br>  res2@gsnRightString = &quot;&quot;<br><br>;-------Frozen<br>res3 = True<br>res3@cnLineColor = &quot;Blue&quot;<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 = &quot;Black&quot;<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 = &quot;blue&quot;<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>