<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>