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

Dave Allured - NOAA Affiliate dave.allured at noaa.gov
Wed Aug 10 14:56:59 MDT 2022


Anie, yes, that is some kind of installation problem.  After this message,
let us move the conversation to the ncl-install mailing list, and drop
ncl-talk.

I am not experienced with conda installs or debugging.  Given that, here
are some general suggestions that may or may not help.  Perhaps someone
else has another suggestion,

* Contact conda support for advice.  Probably try this first.  Mention that
you now have two different errors: can not open libgsl, and contours not
showing.  I think the first error is more important, and should be solved
first.

* Try uninstalling, then reinstalling GSL using conda.

* Uninstall, then reinstall NCL using conda.

* Deactivate your main conda environment.  Then set up a new conda
environment in a different location on the same machine.  Perform a fresh
install of NCL and all dependencies in the alternate conda environment.

* As a last resort, try an alternate installation method.


On Sun, Aug 7, 2022 at 11:00 PM Anie K Lal <Anieklal at cas.iitd.ac.in> wrote:

> Hi,
> I was able to plot it successfully in another PC where I had installed
> ncl. So it is an installation error. I got an error
>
> "ncl: error while loading shared libraries: libgsl.so.25: cannot open
> shared object file: No such file or directory"
>
> In both PCs I installed ncl using conda (I am using Anaconda),
>                     "conda create -n ncl_stable -c conda-forge ncl"
>
> I tried "conda update -n ncl_stable -c conda-forge --all", the error
> still persists.
>
> I am attaching the results of :
>      conda list -n ncl_stable
>      conda --version
>      which ncl
>      ncl -V
>      env  | grep NCARG
>
> In both PCs I installed ncl in same way, for one it worked and other
> it's showing this error!
>
> Kindly advice how to resolve this installation issue.
>
>
> Subject: Re: [ncl-talk] Contours not showing while plotting wrf out file
> Date: 06.08.2022 23:06
>  From: Dave Allured - NOAA Affiliate <dave.allured at noaa.gov>
>
> 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:
>
> https://www.ncl.ucar.edu/Applications/Data/cdf/wrfout_d01_2005-08-28_00:00:00
>
> 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.
>
> On Sat, Aug 6, 2022 at 6:36 AM Anie K Lal <Anieklal at cas.iitd.ac.in>
> wrote:
>
> > Hi,
> >
> > (1) Output from
> >
> > "printVarSummary (hgt)
> > printMinMax (hgt)
> > print (num (ismissing (hgt))
> >
> > printMinMax (hgt at lat2d)
> > printMinMax (hgt at lon2d)"  are:
> >
>
> ----------------------------------------------------------------------------
> > Variable: hgt
> > Type: float
> > Total Size: 88804 bytes
> > 22201 values
> > Number of Dimensions: 2
> > Dimensions and sizes:    [south_north | 149] x [west_east | 149]
> > Coordinates:
> > Number Of Attributes: 8
> > lon2d :    <ARRAY of 22201 elements>
> > lat2d :    <ARRAY of 22201 elements>
> > FieldType :    104
> > MemoryOrder :    XY
> > description :    Terrain Height
> > units :    m
> > stagger :
> > coordinates :    XLONG XLAT XTIME
> > (0)
> > (0)    Terrain Height (m) : min=0   max=6044.23
> > (0)    0
> > (0)
> > (0)    min=16.1061   max=39.9454
> > (0)
> > (0)    min=63.4411   max=90.7589
> >
>
> ----------------------------------------------------------------------------
> >
> > (2) No messages where printed when I ran the program. No errors or
> > warnings.
> >
> > (3) Yes, I ran the original example script by just changing the path
> > and
> > location of my data file. Still its showing the same output, without
> >
> > contours.
> >
> > Kindly advice
> > Thank you
> >
> > Subject: Re: [ncl-talk] Contours not showing while plotting wrf out
> > file
> > Date: 05.08.2022 22:01
> > From: Dave Allured - NOAA Affiliate <dave.allured at noaa.gov>
> > To: Anie K Lal <Anieklal at cas.iitd.ac.in>
> >
> > 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
> > [1]
> >> [1]", "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
> >
> > CAUTION: This email originated from outside of IIT Delhi. Do not
> > click
> > links or open attachments unless you recognize the sender and know
> > the
> > content is safe.
> >
> > Links:
> > ------
> > [1] http://wrfout_d01.nc
>
> CAUTION: This email originated from outside of IIT Delhi. Do not click
> links or open attachments unless you recognize the sender and know the
> content is safe.
>
>
> Links:
> ------
> [1] http://wrfout_d01.nc
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220810/335fff30/attachment.htm>


More information about the ncl-talk mailing list