[ncl-talk] Contours not showing while plotting wrf out file

[Dave Allured - NOAA Affiliate]

Anie, minimal diagnostic information is needed to be able to help you.

(1)  Please show output from:

printVarSummary (hgt)
printMinMax (hgt)
print (num (ismissing (hgt))

printMinMax (hgt at lat2d)
printMinMax (hgt at lon2d)

(2)  Please show any messages that were printed when you ran your program.

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


On Fri, Aug 5, 2022 at 9:23 AM Anie K Lal via ncl-talk <
ncl-talk at mailman.ucar.edu> wrote:

> Hi all,
>
> I am trying to plot WRF-ARW data using gsn_csm_contour_map. I have
> followed the script in
> https://www.ncl.ucar.edu/Applications/wrfgsn.shtml.
>
> But unfortunately I am not getting the contours. I am attaching the
> script I used and the plot I got.
> ________________________________________________________________________-
> d1 = addfile("/home/anie/Documents/WRF_Chem/EXP01/wrfout_d01.nc", "r")
>
>    it        = 0     ; first time step
>    hgt       = wrf_user_getvar(d1,"HGT",it)    ; Terrain elevation
>    hgt at lat2d = wrf_user_getvar(d1,"XLAT",it)   ; latitude/longitude
>    hgt at lon2d = wrf_user_getvar(d1,"XLONG",it)  ; required for plotting
>
>    wks = gsn_open_wks("png","wrf_gsn")
>
> ;---Set some basic plot options
>    res               = True
>
>    res at gsnMaximize   = True   ; maximize plot in frame
>
>   ; res at tiMainString  = filename
>
>    res at cnFillOn      = True
>    res at cnFillPalette = "OceanLakeLandSnow"
>    res at cnLinesOn     = False
>
>    res at mpProjection  = "CylindricalEquidistant"    ; The default
>    res at mpDataBaseVersion = "MediumRes"
>
>    res at gsnAddCyclic      = False
>
> ;---Zoom in on plot
>    res at mpMinLatF     = min(hgt at lat2d)
>    res at mpMaxLatF     = max(hgt at lat2d)
>    res at mpMinLonF     = min(hgt at lon2d)
>    res at mpMaxLonF     = max(hgt at lon2d)
>
>    contour = gsn_csm_contour_map(wks,hgt,res)
>
> _________________________________________________________________________________________
>
> Kindly help me identify the issue.
>
> Looking forward to any kind of help.
>
> Thank you
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220805/8dce9807/attachment.htm>