<div dir="ltr"><div>Anie, minimal diagnostic information is needed to be able to help you.</div><div><br></div><div>(1)  Please show output from:</div><div><br></div><div>printVarSummary (hgt)</div><div>printMinMax (hgt)</div><div>print (num (ismissing (hgt))</div><div><br></div><div>printMinMax (hgt@lat2d)<br>printMinMax (hgt@lon2d)<br></div><div><br></div><div>(2)  Please show any messages that were printed when you ran your program.</div><div><br></div><div>(3)  Please run the original example wrf_gsn_1.ncl with no changes except the path and location of your data file.  Does that make a valid contour map?</div><div dir="ltr"><br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Aug 5, 2022 at 9:23 AM Anie K Lal via ncl-talk <<a href="mailto:ncl-talk@mailman.ucar.edu" target="_blank">ncl-talk@mailman.ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex">Hi all,<br>
<br>
I am trying to plot WRF-ARW data using gsn_csm_contour_map. I have <br>
followed the script in <br>
<a href="https://www.ncl.ucar.edu/Applications/wrfgsn.shtml" rel="noreferrer" target="_blank">https://www.ncl.ucar.edu/Applications/wrfgsn.shtml</a>.<br>
<br>
But unfortunately I am not getting the contours. I am attaching the <br>
script I used and the plot I got.<br>________________________________________________________________________-<br>
d1 = addfile("/home/anie/Documents/WRF_Chem/EXP01/<a href="http://wrfout_d01.nc" rel="noreferrer" target="_blank">wrfout_d01.nc</a>", "r")<br>
<br>
   it        = 0     ; first time step<br>
   hgt       = wrf_user_getvar(d1,"HGT",it)    ; Terrain elevation<br>
   hgt@lat2d = wrf_user_getvar(d1,"XLAT",it)   ; latitude/longitude<br>
   hgt@lon2d = wrf_user_getvar(d1,"XLONG",it)  ; required for plotting<br>
<br>
   wks = gsn_open_wks("png","wrf_gsn")<br>
<br>
;---Set some basic plot options<br>
   res               = True<br>
<br>
   res@gsnMaximize   = True   ; maximize plot in frame<br>
<br>
  ; res@tiMainString  = filename<br>
<br>
   res@cnFillOn      = True<br>
   res@cnFillPalette = "OceanLakeLandSnow"<br>
   res@cnLinesOn     = False<br>
<br>
   res@mpProjection  = "CylindricalEquidistant"    ; The default<br>
   res@mpDataBaseVersion = "MediumRes"<br>
<br>
   res@gsnAddCyclic      = False<br>
<br>
;---Zoom in on plot<br>
   res@mpMinLatF     = min(hgt@lat2d)<br>
   res@mpMaxLatF     = max(hgt@lat2d)<br>
   res@mpMinLonF     = min(hgt@lon2d)<br>
   res@mpMaxLonF     = max(hgt@lon2d)<br>
<br>
   contour = gsn_csm_contour_map(wks,hgt,res)<br>
_________________________________________________________________________________________<br>
<br>
Kindly help me identify the issue.<br>
<br>
Looking forward to any kind of help.<br>
<br>
Thank you<br>
</blockquote></div></div>