<div dir="ltr"><div>Anie, your data ranges look good, and there are no missing values.  Let's test your basic NCL installation.  Try running the same example on a different data file from the NCL website:</div><div><a href="https://www.ncl.ucar.edu/Applications/Data/cdf/wrfout_d01_2005-08-28_00:00:00">https://www.ncl.ucar.edu/Applications/Data/cdf/wrfout_d01_2005-08-28_00:00:00</a><br></div><div><br></div><div>Use the original script example wrf_gsn_1.ncl, with no changes except for the file name.  This works for me.  I get a full color plot with no missing areas.  Let us know if that works.</div><div dir="ltr"><br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sat, Aug 6, 2022 at 6:36 AM Anie K Lal <<a href="mailto:Anieklal@cas.iitd.ac.in">Anieklal@cas.iitd.ac.in</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,<br>
<br>
(1) Output from<br>
<br>
"printVarSummary (hgt)<br>
printMinMax (hgt)<br>
print (num (ismissing (hgt))<br>
<br>
printMinMax (hgt@lat2d)<br>
printMinMax (hgt@lon2d)"  are:<br>
----------------------------------------------------------------------------<br>
Variable: hgt<br>
Type: float<br>
Total Size: 88804 bytes<br>
             22201 values<br>
Number of Dimensions: 2<br>
Dimensions and sizes:    [south_north | 149] x [west_east | 149]<br>
Coordinates:<br>
Number Of Attributes: 8<br>
   lon2d :    <ARRAY of 22201 elements><br>
   lat2d :    <ARRAY of 22201 elements><br>
   FieldType :    104<br>
   MemoryOrder :    XY<br>
   description :    Terrain Height<br>
   units :    m<br>
   stagger :<br>
   coordinates :    XLONG XLAT XTIME<br>
(0)<br>
(0)    Terrain Height (m) : min=0   max=6044.23<br>
(0)    0<br>
(0)<br>
(0)    min=16.1061   max=39.9454<br>
(0)<br>
(0)    min=63.4411   max=90.7589<br>
----------------------------------------------------------------------------<br>
<br>
(2) No messages where printed when I ran the program. No errors or <br>
warnings.<br>
<br>
(3) Yes, I ran the original example script by just changing the path and <br>
location of my data file. Still its showing the same output, without <br>
contours.<br>
<br>
Kindly advice<br>
Thank you<br>
<br><br>
Subject: Re: [ncl-talk] Contours not showing while plotting wrf out file<br>
Date: 05.08.2022 22:01<br>
 From: Dave Allured - NOAA Affiliate <<a href="mailto:dave.allured@noaa.gov" target="_blank">dave.allured@noaa.gov</a>><br>
To: Anie K Lal <<a href="mailto:Anieklal@cas.iitd.ac.in" target="_blank">Anieklal@cas.iitd.ac.in</a>><br>
<br>
Anie, minimal diagnostic information is needed to be able to help you.<br>
<br>
(1)  Please show output from:<br>
<br>
printVarSummary (hgt)<br>
printMinMax (hgt)<br>
print (num (ismissing (hgt))<br>
<br>
printMinMax (hgt@lat2d)<br>
printMinMax (hgt@lon2d)<br>
<br>
(2)  Please show any messages that were printed when you ran your<br>
program.<br>
<br>
(3)  Please run the original example wrf_gsn_1.ncl with no changes<br>
except the path and location of your data file.  Does that make a<br>
valid contour map?<br>
<br>
On Fri, Aug 5, 2022 at 9:23 AM Anie K Lal via ncl-talk<br>
<<a href="mailto:ncl-talk@mailman.ucar.edu" target="_blank">ncl-talk@mailman.ucar.edu</a>> wrote:<br>
<br>
> 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>
________________________________________________________________________-<br>
> d1 = addfile("/home/anie/Documents/WRF_Chem/EXP01/<a href="http://wrfout_d01.nc" rel="noreferrer" target="_blank">wrfout_d01.nc</a><br>
> [1]", "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<br>
> 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>
> <br>
> Kindly help me identify the issue.<br>
> <br>
> Looking forward to any kind of help.<br>
> <br>
> Thank you<br>
<br>
CAUTION: This email originated from outside of IIT Delhi. Do not click<br>
links or open attachments unless you recognize the sender and know the<br>
content is safe.<br><br>
Links:<br>
------<br>
[1] <a href="http://wrfout_d01.nc" rel="noreferrer" target="_blank">http://wrfout_d01.nc</a><br>
</blockquote></div></div>