From zmumba at gmail.com  Mon Aug  1 15:43:57 2022
From: zmumba at gmail.com (Zilore Mumba)
Date: Mon, 1 Aug 2022 23:43:57 +0200
Subject: [ncl-talk] Country names in mask_specs
Message-ID: <CAKS2qbJCZ+ZXxaWDr_wp+NoLiV7xsuy7E_v_cZTfs0v_iOFjgg@mail.gmail.com>

Hello,
I am plotting a map of 2m temperatures over a part of Africa (see graphic
attached). In my mask_specs, the script says
"warning:MapV41DHUpdateDrawList: invalid boundary specification string:
"central-african-republic". There is the same arning for zaire.
The mask_specs I have from the ncl website (mask_4.ncl) are as below.
Can someone advise me how these countries are called for masking?
;mask_specs = (/\

 ;"algeria","angola","angola-exclave-called-cabinda","benin","botswana",\

 ;"burundi","cameroon","central-african-republic","chad","congo","djibouti",\

 ;"egypt","equatorial-guinea","ethiopia","gabon","gambia","ghana","guinea",\
     ;"guinea-bissau","ivory-coast","kenya","lesotho","liberia","libya",\
     ;"madagascar","malawi","mali","mauritania","mauritius","morocco",\

 ;"mozambique","namibia","niger","nigeria","rwanda","senegal","sierra-leone",\

 ;"somalia","south-africa","sudan","swaziland","tanzania","togo","tunisia",\
     ;"uganda","upper-volta","western-sahara","zaire","zambia","zimbabwe"/)
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220801/f51a2d97/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Tmin-da_27.png
Type: image/png
Size: 111057 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220801/f51a2d97/attachment.png>

From zmumba at gmail.com  Tue Aug  2 00:02:00 2022
From: zmumba at gmail.com (Zilore Mumba)
Date: Tue, 2 Aug 2022 08:02:00 +0200
Subject: [ncl-talk] Country names in mask_specs
In-Reply-To: <CAGKRhbFhmjypiRDHFCP9p2akMSeU8GByy608hMDUWf6L2C4rLw@mail.gmail.com>
References: <CAKS2qbJCZ+ZXxaWDr_wp+NoLiV7xsuy7E_v_cZTfs0v_iOFjgg@mail.gmail.com>
 <CAGKRhbFhmjypiRDHFCP9p2akMSeU8GByy608hMDUWf6L2C4rLw@mail.gmail.com>
Message-ID: <CAKS2qbL4x8PPef9_-Ai8DFMBj0CqB+1XztsWR6VwGOP74aNOjg@mail.gmail.com>

Thanks Rick. "central african republic" works without dashes. I am only
remaining with congo dr. I tried "democratic republic of congo", "congo
democratic republic". Both do not work.   I will use a plot without
masking.. just for interest, the script used is attached.

<https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail>
Virus-free.
www.avast.com
<https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail>
<#DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2>

On Tue, Aug 2, 2022 at 12:18 AM Rick Brownrigg <brownrig at ucar.edu> wrote:

> Hi Zilore,
>
> It might help to see the rest of your map resources.   I'm uncertain as to
> which map-names database is being used. Is it possible that Zaire is also
> "Democratic  Republic of the Congo" (forgive my ignorance).  I see "Central
> African Republic" in all the names databases. There are no dashes between
> the names, as there are in the mask_specs. But if that were the issue, then
> other names in that list, like ivory-coast and equator-guinea should also
> result in warnings. I can find no docs that indicate the dashes are needed
> or ignored, nor that the names in the mask are case insensitive?
>
> Rick
>
>
>
> On Mon, Aug 1, 2022 at 3:44 PM Zilore Mumba via ncl-talk <
> ncl-talk at mailman.ucar.edu> wrote:
>
>> Hello,
>> I am plotting a map of 2m temperatures over a part of Africa (see graphic
>> attached). In my mask_specs, the script says
>> "warning:MapV41DHUpdateDrawList: invalid boundary specification string:
>> "central-african-republic". There is the same arning for zaire.
>> The mask_specs I have from the ncl website (mask_4.ncl) are as below.
>> Can someone advise me how these countries are called for masking?
>> ;mask_specs = (/\
>>
>>  ;"algeria","angola","angola-exclave-called-cabinda","benin","botswana",\
>>
>>  ;"burundi","cameroon","central-african-republic","chad","congo","djibouti",\
>>
>>  ;"egypt","equatorial-guinea","ethiopia","gabon","gambia","ghana","guinea",\
>>      ;"guinea-bissau","ivory-coast","kenya","lesotho","liberia","libya",\
>>      ;"madagascar","malawi","mali","mauritania","mauritius","morocco",\
>>
>>  ;"mozambique","namibia","niger","nigeria","rwanda","senegal","sierra-leone",\
>>
>>  ;"somalia","south-africa","sudan","swaziland","tanzania","togo","tunisia",\
>>
>>  ;"uganda","upper-volta","western-sahara","zaire","zambia","zimbabwe"/)
>> _______________________________________________
>> ncl-talk mailing list
>> ncl-talk at mailman.ucar.edu
>> List instructions, subscriber options, unsubscribe:
>> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220802/130e3311/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: 20_PlotT2mMinMax.ncl
Type: application/octet-stream
Size: 5878 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220802/130e3311/attachment.obj>

From setareh.rahimi at gmail.com  Tue Aug  2 11:23:37 2022
From: setareh.rahimi at gmail.com (Setareh Rahimi)
Date: Tue, 2 Aug 2022 21:53:37 +0430
Subject: [ncl-talk] Dust concentration from WRF-chem
Message-ID: <CAN+DFgwBaRS7e+8G4OrFVdYy57Gx_xHJuwzfdcNbsaP+mimBDA@mail.gmail.com>

Dear all NCL users,
I tried to plot dust concentration (from WRF-chem) using the attached
script. But It just plotted for the first hour of simulation while I need
to plot all the time.
Would you please kindly help me to sort this issue out?
 wrfout_d01_2010-07-14_00:00:00
<https://drive.google.com/file/d/1Xm8qa0-oksHiTBBCLP6K42M5pcKvkr1N/view?usp=drive_web>


-- 
S.Rahimi
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220802/7c1a24c0/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Dust_5
Type: application/octet-stream
Size: 815 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220802/7c1a24c0/attachment.obj>

From asphilli at ucar.edu  Wed Aug  3 17:05:55 2022
From: asphilli at ucar.edu (Adam Phillips)
Date: Wed, 3 Aug 2022 17:05:55 -0600
Subject: [ncl-talk] Dust concentration from WRF-chem
In-Reply-To: <CAN+DFgwBaRS7e+8G4OrFVdYy57Gx_xHJuwzfdcNbsaP+mimBDA@mail.gmail.com>
References: <CAN+DFgwBaRS7e+8G4OrFVdYy57Gx_xHJuwzfdcNbsaP+mimBDA@mail.gmail.com>
Message-ID: <CACvWuqAtPYmWuJnHBELXA7vb4e5PwydWY3atLPH7hQxHfttTiw@mail.gmail.com>

Hi Setareh,
Change this:
do nt=0,ntimes,6
to this:
do nt=0,ntimes-1

ntimes = 5, so the former line will run a do loop from 0 to 5 by 6 (=once
when nt = 0), and the latter will run the a do loop from 0 to 5 by 1. (The
1 is implied by not specifying it.)
Hope that helps! If you have any further questions please respond to the
ncl-talk email list.
Adam


On Tue, Aug 2, 2022 at 11:24 AM Setareh Rahimi via ncl-talk <
ncl-talk at mailman.ucar.edu> wrote:

> Dear all NCL users,
> I tried to plot dust concentration (from WRF-chem) using the attached
> script. But It just plotted for the first hour of simulation while I need
> to plot all the time.
> Would you please kindly help me to sort this issue out?
>  wrfout_d01_2010-07-14_00:00:00
> <https://drive.google.com/file/d/1Xm8qa0-oksHiTBBCLP6K42M5pcKvkr1N/view?usp=drive_web>
>
>
> --
> S.Rahimi
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at mailman.ucar.edu
> List instructions, subscriber options, unsubscribe:
> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
>


-- 
Adam Phillips
Associate Scientist,  Climate and Global Dynamics Laboratory, NCAR
www.cgd.ucar.edu/staff/asphilli/   303-497-1726

<http://www.cgd.ucar.edu/staff/asphilli>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220803/4f34bfb1/attachment.htm>

From asphilli at ucar.edu  Wed Aug  3 17:21:28 2022
From: asphilli at ucar.edu (Adam Phillips)
Date: Wed, 3 Aug 2022 17:21:28 -0600
Subject: [ncl-talk] hdf_to_nc_and_mergetime
In-Reply-To: <CAOiHCtTFfF1waYhYWiyPKTrkjjkxOHyO5hVZ308Fzy0d5H+CsQ@mail.gmail.com>
References: <CAOiHCtTFfF1waYhYWiyPKTrkjjkxOHyO5hVZ308Fzy0d5H+CsQ@mail.gmail.com>
Message-ID: <CACvWuqDhZC4MHQxf17D+-3hEcm+LAQdhjLVsQmpyQRqm2SUnVQ@mail.gmail.com>

Hi Kunal,
According to the ListSetType documentation: When one uses the join option,
the input variables from each file must be the same size.
https://www.ncl.ucar.edu/Document/Functions/Built-in/ListSetType.shtml

Your four files vary in the 1st and the 3rd dimension:
; [nCandidate | 8] x [nWavel | 3] x [YDim | 161] x [XDim | 1440]
; [nCandidate | 7] x [nWavel | 3] x [YDim | 161] x [XDim | 1440]
; [nCandidate | 7] x [nWavel | 3] x [YDim | 160] x [XDim | 1440]
; [nCandidate | 7] x [nWavel | 3] x [YDim | 161] x [XDim | 1440]

Hope that helps! If you have any further questions please respond to the
ncl-talk email list.
Adam

On Fri, Jul 29, 2022 at 3:09 AM Kunal Bali via ncl-talk <
ncl-talk at mailman.ucar.edu> wrote:

> Dear NCL users,
>
> I am trying to convert the .nc4 (or hdf) to .nc (or classic .nc) and then
> merge with timesteps by using the attached script (create_time.ncl).
> But, it's not working. I am getting error messages as
>
> [image: Screen Shot 2022-07-29 at 1.04.39 AM.png]
>
> If you need samples data, you can find over
> https://www.dropbox.com/s/978vyxtttutwcq0/Archive.zip?dl=0
>
> Desired variable is "FinalAerosolOpticalDepth"
>
> Thanks
> ---
> Kunal
>
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at mailman.ucar.edu
> List instructions, subscriber options, unsubscribe:
> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
>


-- 
Adam Phillips
Associate Scientist,  Climate and Global Dynamics Laboratory, NCAR
www.cgd.ucar.edu/staff/asphilli/   303-497-1726

<http://www.cgd.ucar.edu/staff/asphilli>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220803/341937dd/attachment-0001.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Screen Shot 2022-07-29 at 1.04.39 AM.png
Type: image/png
Size: 1273023 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220803/341937dd/attachment-0001.png>

From Anieklal at cas.iitd.ac.in  Fri Aug  5 09:23:07 2022
From: Anieklal at cas.iitd.ac.in (Anie K Lal)
Date: Fri, 05 Aug 2022 20:53:07 +0530
Subject: [ncl-talk] Contours not showing while plotting wrf out file
Message-ID: <23648bdd3ecec1489209d64345cb39b4@cas.iitd.ac.in>

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 --------------
A non-text attachment was scrubbed...
Name: wrf_gsn.png
Type: image/png
Size: 43949 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220805/b8e4c702/attachment.png>

From dave.allured at noaa.gov  Fri Aug  5 10:31:08 2022
From: dave.allured at noaa.gov (Dave Allured - NOAA Affiliate)
Date: Fri, 5 Aug 2022 10:31:08 -0600
Subject: [ncl-talk] Contours not showing while plotting wrf out file
In-Reply-To: <23648bdd3ecec1489209d64345cb39b4@cas.iitd.ac.in>
References: <23648bdd3ecec1489209d64345cb39b4@cas.iitd.ac.in>
Message-ID: <CALqwTFNvYZpy5tv2KrSPDJwAFncKsipbhVNCnr4XJHgwKs_VFQ@mail.gmail.com>

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>

From Anieklal at cas.iitd.ac.in  Sat Aug  6 06:36:34 2022
From: Anieklal at cas.iitd.ac.in (Anie K Lal)
Date: Sat, 06 Aug 2022 18:06:34 +0530
Subject: [ncl-talk] Contours not showing while plotting wrf out file
In-Reply-To: <CALqwTFNvYZpy5tv2KrSPDJwAFncKsipbhVNCnr4XJHgwKs_VFQ@mail.gmail.com>
References: <23648bdd3ecec1489209d64345cb39b4@cas.iitd.ac.in>
 <CALqwTFNvYZpy5tv2KrSPDJwAFncKsipbhVNCnr4XJHgwKs_VFQ@mail.gmail.com>
Message-ID: <de9d9ed779d329a9f616daf1969e3516@cas.iitd.ac.in>

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





-------- Original Message --------
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]", "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

From dave.allured at noaa.gov  Sat Aug  6 11:36:07 2022
From: dave.allured at noaa.gov (Dave Allured - NOAA Affiliate)
Date: Sat, 6 Aug 2022 11:36:07 -0600
Subject: [ncl-talk] Contours not showing while plotting wrf out file
In-Reply-To: <de9d9ed779d329a9f616daf1969e3516@cas.iitd.ac.in>
References: <23648bdd3ecec1489209d64345cb39b4@cas.iitd.ac.in>
 <CALqwTFNvYZpy5tv2KrSPDJwAFncKsipbhVNCnr4XJHgwKs_VFQ@mail.gmail.com>
 <de9d9ed779d329a9f616daf1969e3516@cas.iitd.ac.in>
Message-ID: <CALqwTFPzg-smaDdihxUid0q75pbcU+A0gvq3G2HVaQOTadRXAw@mail.gmail.com>

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]", "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
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220806/8f7878c9/attachment.htm>

From Anieklal at cas.iitd.ac.in  Sun Aug  7 00:09:54 2022
From: Anieklal at cas.iitd.ac.in (Anie K Lal)
Date: Sun, 07 Aug 2022 11:39:54 +0530
Subject: [ncl-talk] Distorted projection
Message-ID: <f026a99b12626c04e8be9c3e9c52aba0@cas.iitd.ac.in>

Hi all,

I am trying to make a Terrain plot with domain boxes. Even though I was 
successful in getting the plot, the projection is not looking good. I am 
attaching my script and the plot i got.

Looking forward to any kind of help.

Thank you
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: domain_test01.ncl
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220807/28a2542b/attachment.ksh>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: wrf_overlay01.png
Type: image/png
Size: 94640 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220807/28a2542b/attachment.png>

From barry.h.lynn at gmail.com  Sun Aug  7 03:33:05 2022
From: barry.h.lynn at gmail.com (Barry Lynn)
Date: Sun, 7 Aug 2022 12:33:05 +0300
Subject: [ncl-talk] Distorted projection
In-Reply-To: <f026a99b12626c04e8be9c3e9c52aba0@cas.iitd.ac.in>
References: <f026a99b12626c04e8be9c3e9c52aba0@cas.iitd.ac.in>
Message-ID: <CAN7Bv24JSxjtEtfx38+-Tp=JCT0OG2rqUdjxPJhs=gyMO7UeNg@mail.gmail.com>

Hi Anie:

That's because your model data is not on a Mercator grid, but a Lambert,
probably.  I believe that the blank area indicates that you have no data in
that location.  Try shrinking your data.  The only question then, is, is
your background map being projected on the same grid.

On Sun, Aug 7, 2022 at 9:10 AM Anie K Lal via ncl-talk <
ncl-talk at mailman.ucar.edu> wrote:

> Hi all,
>
> I am trying to make a Terrain plot with domain boxes. Even though I was
> successful in getting the plot, the projection is not looking good. I am
> attaching my script and the plot i got.
>
> Looking forward to any kind of help.
>
> Thank you_______________________________________________
> ncl-talk mailing list
> ncl-talk at mailman.ucar.edu
> List instructions, subscriber options, unsubscribe:
> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
>


-- 

Barry H. Lynn, Ph.D
Senior Scientist, Lecturer,
The Institute of Earth Sciences,
The Hebrew University of Jerusalem,
Givat Ram, Jerusalem 91904, Israel
Tel: 972 547 231 170
Fax: (972)-25662581

Weather It Is, LTD
Weather and Climate Focus
https://weather-it-is.com <http://weather-it-is.com>
Jerusalem, Israel
Local: 02 930 9525
Cell: 054 7 231 170
Int-IS: x972 2 930 9525
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220807/a0c93fc1/attachment.htm>

From jpclark186000 at gmail.com  Sun Aug  7 10:02:01 2022
From: jpclark186000 at gmail.com (Joseph Clark)
Date: Sun, 7 Aug 2022 12:02:01 -0400
Subject: [ncl-talk] How to place polymarker ABOVE the x-axis
Message-ID: <CAD_pzv4u=3vd8E3a1+umyuA_RAoVbNBQawZimDyGgEu3juvDWA@mail.gmail.com>

Hi NCL users,

Below I attached a script and pdf showing the problem. In the PDF, you can
see that I am trying to overlay a point on the x-intercept. The point is
half-way behind the x-axis, however I want the point *on top* of the
x-axis. I'm using version 6.3.0 of NCL and would prefer a solution that
doesn't require me to update the software. (I'm not the administrator).

Any help is appreciated!

Joseph
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220807/dc8b8413/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: dummy_code.ncl
Type: application/octet-stream
Size: 548 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220807/dc8b8413/attachment.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: dummy.pdf
Type: application/pdf
Size: 13618 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220807/dc8b8413/attachment.pdf>

From jpclark186000 at gmail.com  Sun Aug  7 10:10:22 2022
From: jpclark186000 at gmail.com (Joseph Clark)
Date: Sun, 7 Aug 2022 12:10:22 -0400
Subject: [ncl-talk] How to place polymarker ABOVE the x-axis
In-Reply-To: <CAD_pzv4u=3vd8E3a1+umyuA_RAoVbNBQawZimDyGgEu3juvDWA@mail.gmail.com>
References: <CAD_pzv4u=3vd8E3a1+umyuA_RAoVbNBQawZimDyGgEu3juvDWA@mail.gmail.com>
Message-ID: <CAD_pzv4dHhG-sxJs3uZOGsnQ-AQkj_SR9eO6KB9p8szSY_SLng@mail.gmail.com>

Solved:

mres at gsClipOn = False

Posting here in case others have a similar problem.

On Sun, Aug 7, 2022 at 12:02 PM Joseph Clark <jpclark186000 at gmail.com>
wrote:

> Hi NCL users,
>
> Below I attached a script and pdf showing the problem. In the PDF, you can
> see that I am trying to overlay a point on the x-intercept. The point is
> half-way behind the x-axis, however I want the point *on top* of the
> x-axis. I'm using version 6.3.0 of NCL and would prefer a solution that
> doesn't require me to update the software. (I'm not the administrator).
>
> Any help is appreciated!
>
> Joseph
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220807/4c83f16a/attachment.htm>

From Anieklal at cas.iitd.ac.in  Sun Aug  7 22:43:18 2022
From: Anieklal at cas.iitd.ac.in (Anie K Lal)
Date: Mon, 08 Aug 2022 10:13:18 +0530
Subject: [ncl-talk] Distorted projection
In-Reply-To: <CAN7Bv24JSxjtEtfx38+-Tp=JCT0OG2rqUdjxPJhs=gyMO7UeNg@mail.gmail.com>
References: <f026a99b12626c04e8be9c3e9c52aba0@cas.iitd.ac.in>
 <CAN7Bv24JSxjtEtfx38+-Tp=JCT0OG2rqUdjxPJhs=gyMO7UeNg@mail.gmail.com>
Message-ID: <4a653a65a27335559357c6e19bc12e90@cas.iitd.ac.in>

My problem got solved after using "res = 
set_mp_wrf_map_resources(a,res)".

Thank you.


-------- Original Message --------
Subject: Re: [ncl-talk] Distorted projection
Date: 07.08.2022 15:03
 From: Barry Lynn <barry.h.lynn at gmail.com>
To: Anie K Lal <Anieklal at cas.iitd.ac.in>

Hi Anie:

That's because your model data is not on a Mercator grid, but a
Lambert, probably.  I believe that the blank area indicates that you
have no data in that location.  Try shrinking your data.  The only
question then, is, is your background map being projected on the same
grid.

On Sun, Aug 7, 2022 at 9:10 AM Anie K Lal via ncl-talk
<ncl-talk at mailman.ucar.edu> wrote:

> Hi all,
> 
> I am trying to make a Terrain plot with domain boxes. Even though I
> was
> successful in getting the plot, the projection is not looking good.
> I am
> attaching my script and the plot i got.
> 
> Looking forward to any kind of help.
> 
> Thank you_______________________________________________
> ncl-talk mailing list
> ncl-talk at mailman.ucar.edu
> List instructions, subscriber options, unsubscribe:
> https://mailman.ucar.edu/mailman/listinfo/ncl-talk

-- 

Barry H. Lynn, Ph.D

Senior Scientist, Lecturer,

The Institute of Earth Sciences,
The Hebrew University of Jerusalem,
Givat Ram, Jerusalem 91904, Israel
Tel: 972 547 231 170
Fax: (972)-25662581

Weather It Is, LTD
Weather and Climate Focus
https://weather-it-is.com [1]
Jerusalem, Israel
Local: 02 930 9525
Cell: 054 7 231 170
Int-IS: x972 2 930 9525

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://weather-it-is.com

From Anieklal at cas.iitd.ac.in  Sun Aug  7 23:00:28 2022
From: Anieklal at cas.iitd.ac.in (Anie K Lal)
Date: Mon, 08 Aug 2022 10:30:28 +0530
Subject: [ncl-talk] Contours not showing while plotting wrf out file
In-Reply-To: <CALqwTFPzg-smaDdihxUid0q75pbcU+A0gvq3G2HVaQOTadRXAw@mail.gmail.com>
References: <23648bdd3ecec1489209d64345cb39b4@cas.iitd.ac.in>
 <CALqwTFNvYZpy5tv2KrSPDJwAFncKsipbhVNCnr4XJHgwKs_VFQ@mail.gmail.com>
 <de9d9ed779d329a9f616daf1969e3516@cas.iitd.ac.in>
 <CALqwTFPzg-smaDdihxUid0q75pbcU+A0gvq3G2HVaQOTadRXAw@mail.gmail.com>
Message-ID: <41c1b3172c8f7d2424e2d8e758b89daf@cas.iitd.ac.in>

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.




-------- Original Message --------
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>
To: Anie K Lal <Anieklal at cas.iitd.ac.in>

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 embedded and charset-unspecified text was scrubbed...
Name: Ncl_test
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220808/d938c73d/attachment.ksh>

From barry.h.lynn at gmail.com  Wed Aug 10 02:09:38 2022
From: barry.h.lynn at gmail.com (Barry Lynn)
Date: Wed, 10 Aug 2022 11:09:38 +0300
Subject: [ncl-talk] Is this a bug (or a mirage)?
Message-ID: <CAN7Bv26qrKw3Svm2H+uUBaXa0w-Rb+7MFr_ORnFLzR=d+adiQw@mail.gmail.com>

Hi:

I can't figure out how this can be.  Why isn't "max"working and why are the
values divided by 1000?

cold_hours_cost = new(2,float)

hot_hours_cost = new(2,float)

p_cold = T*cold_hours

p_hot = T*hot_hours

cold_hours_cost(1) = p_cold

hot_hours_cost(1) = p_hot

 print("cold_hours = " + cold_hours_cost)

 print("hot_hours = " + hot_hours_cost)

max_c = max(cold_hours)

max_h = max(hot_hours)

print("max(cold_hours) = " + max_c)

print("max(hot_hours) = " + max_h)

(0) cold_hours = 0

(1) cold_hours = -4111.11

(0) hot_hours = 0

(1) hot_hours = 4111.11

(0) max(cold_hours) = -4.11111

(0) max(hot_hours) = 4.11111

-- 

Barry H. Lynn, Ph.D
Senior Scientist, Lecturer,
The Institute of Earth Sciences,
The Hebrew University of Jerusalem,
Givat Ram, Jerusalem 91904, Israel
Tel: 972 547 231 170
Fax: (972)-25662581

Weather It Is, LTD
Weather and Climate Focus
https://weather-it-is.com <http://weather-it-is.com>
Jerusalem, Israel
Local: 02 930 9525
Cell: 054 7 231 170
Int-IS: x972 2 930 9525
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220810/29c403be/attachment.htm>

From brownrig at ucar.edu  Wed Aug 10 07:35:45 2022
From: brownrig at ucar.edu (Rick Brownrigg)
Date: Wed, 10 Aug 2022 07:35:45 -0600
Subject: [ncl-talk] Is this a bug (or a mirage)?
In-Reply-To: <CAN7Bv26qrKw3Svm2H+uUBaXa0w-Rb+7MFr_ORnFLzR=d+adiQw@mail.gmail.com>
References: <CAN7Bv26qrKw3Svm2H+uUBaXa0w-Rb+7MFr_ORnFLzR=d+adiQw@mail.gmail.com>
Message-ID: <CAGKRhbEe=ydAzy1CkWR_dzLj5C2k-DSyDicXM0ShBhCqkztk5w@mail.gmail.com>

What are the values of T, cold_hours, and hot_hours?

On Wed, Aug 10, 2022 at 2:10 AM Barry Lynn via ncl-talk <
ncl-talk at mailman.ucar.edu> wrote:

> Hi:
>
> I can't figure out how this can be.  Why isn't "max"working and why are
> the values divided by 1000?
>
> cold_hours_cost = new(2,float)
>
> hot_hours_cost = new(2,float)
>
> p_cold = T*cold_hours
>
> p_hot = T*hot_hours
>
> cold_hours_cost(1) = p_cold
>
> hot_hours_cost(1) = p_hot
>
>  print("cold_hours = " + cold_hours_cost)
>
>  print("hot_hours = " + hot_hours_cost)
>
> max_c = max(cold_hours)
>
> max_h = max(hot_hours)
>
> print("max(cold_hours) = " + max_c)
>
> print("max(hot_hours) = " + max_h)
>
> (0) cold_hours = 0
>
> (1) cold_hours = -4111.11
>
> (0) hot_hours = 0
>
> (1) hot_hours = 4111.11
>
> (0) max(cold_hours) = -4.11111
>
> (0) max(hot_hours) = 4.11111
>
> --
>
> Barry H. Lynn, Ph.D
> Senior Scientist, Lecturer,
> The Institute of Earth Sciences,
> The Hebrew University of Jerusalem,
> Givat Ram, Jerusalem 91904, Israel
> Tel: 972 547 231 170
> Fax: (972)-25662581
>
> Weather It Is, LTD
> Weather and Climate Focus
> https://weather-it-is.com <http://weather-it-is.com>
> Jerusalem, Israel
> Local: 02 930 9525
> Cell: 054 7 231 170
> Int-IS: x972 2 930 9525
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at mailman.ucar.edu
> List instructions, subscriber options, unsubscribe:
> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220810/ce2233b7/attachment.htm>

From dave.allured at noaa.gov  Wed Aug 10 14:56:59 2022
From: dave.allured at noaa.gov (Dave Allured - NOAA Affiliate)
Date: Wed, 10 Aug 2022 14:56:59 -0600
Subject: [ncl-talk] Contours not showing while plotting wrf out file
In-Reply-To: <41c1b3172c8f7d2424e2d8e758b89daf@cas.iitd.ac.in>
References: <23648bdd3ecec1489209d64345cb39b4@cas.iitd.ac.in>
 <CALqwTFNvYZpy5tv2KrSPDJwAFncKsipbhVNCnr4XJHgwKs_VFQ@mail.gmail.com>
 <de9d9ed779d329a9f616daf1969e3516@cas.iitd.ac.in>
 <CALqwTFPzg-smaDdihxUid0q75pbcU+A0gvq3G2HVaQOTadRXAw@mail.gmail.com>
 <41c1b3172c8f7d2424e2d8e758b89daf@cas.iitd.ac.in>
Message-ID: <CALqwTFPg0Jz45jEgPnJ7MMwSKNJkK9zeT00kG6qLHXG60xS0XQ@mail.gmail.com>

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>

From setareh.rahimi at gmail.com  Wed Aug 10 15:57:14 2022
From: setareh.rahimi at gmail.com (Setareh Rahimi)
Date: Thu, 11 Aug 2022 02:27:14 +0430
Subject: [ncl-talk] Aerosol optical depth from WRF-chem
Message-ID: <CAN+DFgzHYAwBZEapge5z4O=saVjCoymvoYVxcxqgDYU5OrFa9w@mail.gmail.com>

Dear NCL users,

I need to plot aerosol optical depth(AOD) from WRF-chem outputs. I follow
the structure( script attached), however,
 wrfout_d01_2010-07-14_00:00:00
<https://drive.google.com/file/d/1-9zPQ-P6cmvrYs86CCeiNB-VU66Nuluw/view?usp=drive_web>
the NCL returned this error:
Error: scalar_field: If the input data is 1-dimensional, you must set
sfXArray and sfYArray to 1-dimensional arrays of the same length.
My question is how can I set sfXArray and sfYArray to 1-dimensional arrays
of the same length, please?
Many thanks in advance,
Best wishes,
-- 
S.Rahimi
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220811/f89393d6/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: dust-ncl
Type: application/octet-stream
Size: 1301 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220811/f89393d6/attachment.obj>

From asphilli at ucar.edu  Thu Aug 11 14:38:21 2022
From: asphilli at ucar.edu (Adam Phillips)
Date: Thu, 11 Aug 2022 14:38:21 -0600
Subject: [ncl-talk] Aerosol optical depth from WRF-chem
In-Reply-To: <CAN+DFgzHYAwBZEapge5z4O=saVjCoymvoYVxcxqgDYU5OrFa9w@mail.gmail.com>
References: <CAN+DFgzHYAwBZEapge5z4O=saVjCoymvoYVxcxqgDYU5OrFa9w@mail.gmail.com>
Message-ID: <CACvWuqDaPEsbYk45Ermn_V7k5RdduZ0nUb_fHs=-Y5DVTNtBLw@mail.gmail.com>

Hi Setareh,
The array that you are plotting, AOD550_2D,  is a scalar value. I do not
believe that gsn_csm_contour_map can plot a single value.
The end of your script:
AOD550_3D = a2* ((400/550)^angstrom_exponent)
 printVarSummary(AOD550_3D)
AOD550_2D = sum(AOD550_3D)
 printVarSummary(AOD550_2D)
res at tiMainString = "WRF-CHEM (DUST_5) " + times(nt)
res at gsnLeftString = x at description
plot = gsn_csm_contour_map(wks,AOD550_2D,res)

Your AO550_3D array is dimensioned 35 x 40 x 40, but there is no coordinate
variable information. I would suggest adding
copy_VarMeta(a2,AOD550_3D)
to copy metadata from your a2 array to AOD550_3D. To plot the spatial
array, you will also need to read in the XLONG/XLAT arrays (which are 3D)
and set them to lat1d/lon1d arrays attached as attributes:
AOD550_3D at lat2d = f->XLAT(nt,:,:)
AOD550_3D at lon2d = f->XLONG(nt,:,:)

Hope that gets you going down the right path. If you have any further
questions, please respond to the ncl-talk email list.
Adam




On Wed, Aug 10, 2022 at 3:58 PM Setareh Rahimi via ncl-talk <
ncl-talk at mailman.ucar.edu> wrote:

> Dear NCL users,
>
> I need to plot aerosol optical depth(AOD) from WRF-chem outputs. I follow
> the structure( script attached), however,
>  wrfout_d01_2010-07-14_00:00:00
> <https://drive.google.com/file/d/1-9zPQ-P6cmvrYs86CCeiNB-VU66Nuluw/view?usp=drive_web>
> the NCL returned this error:
> Error: scalar_field: If the input data is 1-dimensional, you must set
> sfXArray and sfYArray to 1-dimensional arrays of the same length.
> My question is how can I set sfXArray and sfYArray to 1-dimensional
> arrays of the same length, please?
> Many thanks in advance,
> Best wishes,
> --
> S.Rahimi
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at mailman.ucar.edu
> List instructions, subscriber options, unsubscribe:
> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
>


-- 
Adam Phillips
Associate Scientist,  Climate and Global Dynamics Laboratory, NCAR
www.cgd.ucar.edu/staff/asphilli/   303-497-1726

<http://www.cgd.ucar.edu/staff/asphilli>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220811/50b78b85/attachment.htm>

From debasish.hazra5 at gmail.com  Wed Aug 17 04:20:43 2022
From: debasish.hazra5 at gmail.com (Debasish Hazra)
Date: Wed, 17 Aug 2022 06:20:43 -0400
Subject: [ncl-talk] missing value and RMSE
Message-ID: <CAJ0V9jqYwGjipP6-Rhs_76aVGsvxO3FdnJE1TSJng3WLb1oREQ@mail.gmail.com>

Hi,
I am calculating RMSE timeseries for series of model forecasts and
observations. However, observation is satellite measurement and contains
missing values. If I use ncl dim_rmsd for RMSE calculation, it does not
include data gaps in observations, right ?
Is there any other way to calculate RMSE, particularly applicable here. Any
help is appreciated

Thanks
Debasish
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220817/855b9f20/attachment.htm>

From shea at ucar.edu  Wed Aug 17 21:00:09 2022
From: shea at ucar.edu (Dennis Shea)
Date: Wed, 17 Aug 2022 21:00:09 -0600
Subject: [ncl-talk] missing value and RMSE
In-Reply-To: <CAJ0V9jqYwGjipP6-Rhs_76aVGsvxO3FdnJE1TSJng3WLb1oREQ@mail.gmail.com>
References: <CAJ0V9jqYwGjipP6-Rhs_76aVGsvxO3FdnJE1TSJng3WLb1oREQ@mail.gmail.com>
Message-ID: <CAOF1d_7uRPRWhAucqs7XG2gpy-WG=_h6ZrvLHUBiLpCUu=N5KA@mail.gmail.com>

"data gaps"  not sure what you mean.

You can see the dim_rmsd code in contributed.ncl

On Wed, Aug 17, 2022 at 4:21 AM Debasish Hazra via ncl-talk <
ncl-talk at mailman.ucar.edu> wrote:

> Hi,
> I am calculating RMSE timeseries for series of model forecasts and
> observations. However, observation is satellite measurement and contains
> missing values. If I use ncl dim_rmsd for RMSE calculation, it does not
> include data gaps in observations, right ?
> Is there any other way to calculate RMSE, particularly applicable here.
> Any help is appreciated
>
> Thanks
> Debasish
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at mailman.ucar.edu
> List instructions, subscriber options, unsubscribe:
> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220817/44cedf04/attachment.htm>

From manisha.ganeshan86 at gmail.com  Thu Aug 18 13:57:58 2022
From: manisha.ganeshan86 at gmail.com (Manisha Ganeshan)
Date: Thu, 18 Aug 2022 15:57:58 -0400
Subject: [ncl-talk] saturation vapor pressure for mixhum_ptrh function
Message-ID: <CA+o9MVt9HSg92+HQC8w2rg7FTW9woQxgvNhW7Pm=E3fQG3r9fg@mail.gmail.com>

Hello,

I'd like to use the *mixhum_ptrh* function to obtain the specific humidity
(given temperature, pressure, and relative humidity). I'm wondering what
equation this function assumes to estimate the saturation vapor pressure,
and whether there are any differences in the assumption regarding
saturation w.r.t liquid water versus saturation w.r.t. ice?

Thanks,
Manisha

-- 
Dr. Manisha Ganeshan
Morgan State University
Maryland, U.S.A.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220818/4464d588/attachment.htm>

From dave.allured at noaa.gov  Thu Aug 18 14:26:14 2022
From: dave.allured at noaa.gov (Dave Allured - NOAA Affiliate)
Date: Thu, 18 Aug 2022 14:26:14 -0600
Subject: [ncl-talk] saturation vapor pressure for mixhum_ptrh function
In-Reply-To: <CA+o9MVt9HSg92+HQC8w2rg7FTW9woQxgvNhW7Pm=E3fQG3r9fg@mail.gmail.com>
References: <CA+o9MVt9HSg92+HQC8w2rg7FTW9woQxgvNhW7Pm=E3fQG3r9fg@mail.gmail.com>
Message-ID: <CALqwTFNPXz8dmH6sxbjEKKJ7b77Xn9huN99zY8yDyY1b8adnzQ@mail.gmail.com>

General advice:

(1)  Check the NCL on-line function docs first.

(2)  In some cases, the NCL examples pages have supplemental documentation
for NCL functions.

(3)  Search the NCL source code.  Frequently there is extra documentation
within source code routines for scientific functions, because these were
often imported from other documented scientific packages.  There is also
the code itself, of course.

The internal code for this function is
*ncl_ncarg-6.6.2/ni/src/lib/nfpfort/mixhum_ptrh.f*.  This is called by a
wrapper function in *ncl_ncarg-6.6.2/ni/src/lib/nfp/mixhumW.c*.


On Thu, Aug 18, 2022 at 1:58 PM Manisha Ganeshan via ncl-talk <
ncl-talk at mailman.ucar.edu> wrote:

> Hello,
>
> I'd like to use the *mixhum_ptrh* function to obtain the specific
> humidity (given temperature, pressure, and relative humidity). I'm
> wondering what equation this function assumes to estimate the saturation
> vapor pressure, and whether there are any differences in the assumption
> regarding saturation w.r.t liquid water versus saturation w.r.t. ice?
>
> Thanks,
> Manisha
> --
> Dr. Manisha Ganeshan
> Morgan State University
> Maryland, U.S.A.
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220818/56514812/attachment.htm>

From setareh.rahimi at gmail.com  Fri Aug 19 13:36:48 2022
From: setareh.rahimi at gmail.com (Setareh Rahimi)
Date: Sat, 20 Aug 2022 00:06:48 +0430
Subject: [ncl-talk] Aerosol optical depth from WRF-chem
In-Reply-To: <CACvWuqDaPEsbYk45Ermn_V7k5RdduZ0nUb_fHs=-Y5DVTNtBLw@mail.gmail.com>
References: <CAN+DFgzHYAwBZEapge5z4O=saVjCoymvoYVxcxqgDYU5OrFa9w@mail.gmail.com>
 <CACvWuqDaPEsbYk45Ermn_V7k5RdduZ0nUb_fHs=-Y5DVTNtBLw@mail.gmail.com>
Message-ID: <CAN+DFgw+qTPubD0uGpqJsiwo6qctd4Ud4QPbssrqz4-2W6YCWw@mail.gmail.com>

Dear Adam,
So many thanks for your help and advice. And also sorry for the slow reply.
Best wishes,

On Fri, Aug 12, 2022 at 1:08 AM Adam Phillips <asphilli at ucar.edu> wrote:

> Hi Setareh,
> The array that you are plotting, AOD550_2D,  is a scalar value. I do not
> believe that gsn_csm_contour_map can plot a single value.
> The end of your script:
> AOD550_3D = a2* ((400/550)^angstrom_exponent)
>  printVarSummary(AOD550_3D)
> AOD550_2D = sum(AOD550_3D)
>  printVarSummary(AOD550_2D)
> res at tiMainString = "WRF-CHEM (DUST_5) " + times(nt)
> res at gsnLeftString = x at description
> plot = gsn_csm_contour_map(wks,AOD550_2D,res)
>
> Your AO550_3D array is dimensioned 35 x 40 x 40, but there is no
> coordinate variable information. I would suggest adding
> copy_VarMeta(a2,AOD550_3D)
> to copy metadata from your a2 array to AOD550_3D. To plot the spatial
> array, you will also need to read in the XLONG/XLAT arrays (which are 3D)
> and set them to lat1d/lon1d arrays attached as attributes:
> AOD550_3D at lat2d = f->XLAT(nt,:,:)
> AOD550_3D at lon2d = f->XLONG(nt,:,:)
>
> Hope that gets you going down the right path. If you have any further
> questions, please respond to the ncl-talk email list.
> Adam
>
>
>
>
> On Wed, Aug 10, 2022 at 3:58 PM Setareh Rahimi via ncl-talk <
> ncl-talk at mailman.ucar.edu> wrote:
>
>> Dear NCL users,
>>
>> I need to plot aerosol optical depth(AOD) from WRF-chem outputs. I follow
>> the structure( script attached), however,
>>  wrfout_d01_2010-07-14_00:00:00
>> <https://drive.google.com/file/d/1-9zPQ-P6cmvrYs86CCeiNB-VU66Nuluw/view?usp=drive_web>
>> the NCL returned this error:
>> Error: scalar_field: If the input data is 1-dimensional, you must set
>> sfXArray and sfYArray to 1-dimensional arrays of the same length.
>> My question is how can I set sfXArray and sfYArray to 1-dimensional
>> arrays of the same length, please?
>> Many thanks in advance,
>> Best wishes,
>> --
>> S.Rahimi
>>
>> _______________________________________________
>> ncl-talk mailing list
>> ncl-talk at mailman.ucar.edu
>> List instructions, subscriber options, unsubscribe:
>> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>
>
> --
> Adam Phillips
> Associate Scientist,  Climate and Global Dynamics Laboratory, NCAR
> www.cgd.ucar.edu/staff/asphilli/   303-497-1726
>
> <http://www.cgd.ucar.edu/staff/asphilli>
>


-- 
S.Rahimi
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220820/6de848a7/attachment.htm>

From setareh.rahimi at gmail.com  Fri Aug 19 13:41:28 2022
From: setareh.rahimi at gmail.com (Setareh Rahimi)
Date: Sat, 20 Aug 2022 00:11:28 +0430
Subject: [ncl-talk] Color bar
Message-ID: <CAN+DFgx=yxjEg_02fzp+=pakNv6rE0+6ttaB5y5pF3PUhMkZZw@mail.gmail.com>

Dear all NCL users,
I need to plot my map so that the color bar starts from a specific color. I
attached my plot. How can I set the color bar to start from light blue?
So many thanks in advance,
Best wishes,

-- 
S.Rahimi
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220820/75357978/attachment-0001.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: aod_20110705.000002.png
Type: image/png
Size: 159031 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220820/75357978/attachment-0001.png>

From olagueralyndonmark429 at gmail.com  Mon Aug 22 18:02:45 2022
From: olagueralyndonmark429 at gmail.com (Lyndz)
Date: Tue, 23 Aug 2022 08:02:45 +0800
Subject: [ncl-talk] How to extract diabatic heating from Potential Vorticity
 function in NCL
Message-ID: <CA+V+8dWtPd28L8HAM+QfVeW0uHs+MY9ZRkuy6AN9DQxquUCWNQ@mail.gmail.com>

Dear NCL-experts,

I saw this script that computes the Potential vorticity in isobaric levels.
https://www.ncl.ucar.edu/Applications/pot_vort.shtml

Based on the examples, the static stability and potential temperatures can
be extracted. I would like to ask how to extract the diabatic heating
(based on potential temperature) from this?

I'll appreciate any help on this.

Sincerely,
Lyndz
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220823/a17b4589/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Screenshot from 2022-08-23 08-02-04.png
Type: image/png
Size: 73446 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220823/a17b4589/attachment.png>

From shea at ucar.edu  Mon Aug 22 20:10:14 2022
From: shea at ucar.edu (Dennis Shea)
Date: Mon, 22 Aug 2022 20:10:14 -0600
Subject: [ncl-talk] How to extract diabatic heating from Potential
 Vorticity function in NCL
In-Reply-To: <CA+V+8dWtPd28L8HAM+QfVeW0uHs+MY9ZRkuy6AN9DQxquUCWNQ@mail.gmail.com>
References: <CA+V+8dWtPd28L8HAM+QfVeW0uHs+MY9ZRkuy6AN9DQxquUCWNQ@mail.gmail.com>
Message-ID: <CAOF1d_517=L324bzwU96gtmZZe0uOyF-D0HE3c_AD+8wDYWy9g@mail.gmail.com>

NCL's function is for isobaric (pressure) surfaces.

Your attachment uses a formulation for sigma  coordinates.

I don't think the components calculated in NCL's function can readily be
used to
"extract the diabatic heating (based on potential temperature)"

=============
Sorry



On Mon, Aug 22, 2022 at 6:08 PM Lyndz via ncl-talk <
ncl-talk at mailman.ucar.edu> wrote:

> Dear NCL-experts,
>
> I saw this script that computes the Potential vorticity in isobaric levels.
> https://www.ncl.ucar.edu/Applications/pot_vort.shtml
>
> Based on the examples, the static stability and potential temperatures can
> be extracted. I would like to ask how to extract the diabatic heating
> (based on potential temperature) from this?
>
> I'll appreciate any help on this.
>
> Sincerely,
> Lyndz
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at mailman.ucar.edu
> List instructions, subscriber options, unsubscribe:
> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220822/9d16630e/attachment.htm>

From olagueralyndonmark429 at gmail.com  Mon Aug 22 20:42:42 2022
From: olagueralyndonmark429 at gmail.com (Lyndz)
Date: Tue, 23 Aug 2022 10:42:42 +0800
Subject: [ncl-talk] How to extract diabatic heating from Potential
 Vorticity function in NCL
In-Reply-To: <CAOF1d_517=L324bzwU96gtmZZe0uOyF-D0HE3c_AD+8wDYWy9g@mail.gmail.com>
References: <CAOF1d_517=L324bzwU96gtmZZe0uOyF-D0HE3c_AD+8wDYWy9g@mail.gmail.com>
Message-ID: <C463FD3E-4445-48D6-8E24-E41D78BC5A65@gmail.com>

Hi Sir Dennis,

Thanks for this. 
I just remembered you recommended a post for for this before:
http://hannahlab.org/a-comparison-of-methods-for-estimating-diabatic-heating/

Please ignore the email that I sent.

--Lyndz

Sent from my iPhone

> On Aug 23, 2022, at 10:10 AM, Dennis Shea <shea at ucar.edu> wrote:
> 
> ?
> NCL's function is for isobaric (pressure) surfaces.
> 
> Your attachment uses a formulation for sigma  coordinates. 
> 
> I don't think the components calculated in NCL's function can readily be used to
> "extract the diabatic heating (based on potential temperature)"
> 
> =============
> Sorry
> 
> 
> 
>> On Mon, Aug 22, 2022 at 6:08 PM Lyndz via ncl-talk <ncl-talk at mailman.ucar.edu> wrote:
>> Dear NCL-experts,
>> 
>> I saw this script that computes the Potential vorticity in isobaric levels.
>> https://www.ncl.ucar.edu/Applications/pot_vort.shtml
>> 
>> Based on the examples, the static stability and potential temperatures can be extracted. I would like to ask how to extract the diabatic heating (based on potential temperature) from this?
>> 
>> I'll appreciate any help on this.
>> 
>> Sincerely,
>> Lyndz
>> _______________________________________________
>> ncl-talk mailing list
>> ncl-talk at mailman.ucar.edu
>> List instructions, subscriber options, unsubscribe:
>> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220823/c7e296be/attachment.htm>

From olagueralyndonmark429 at gmail.com  Tue Aug 23 22:27:05 2022
From: olagueralyndonmark429 at gmail.com (Lyndz)
Date: Wed, 24 Aug 2022 12:27:05 +0800
Subject: [ncl-talk] About Unit conversion handling
Message-ID: <CA+V+8dUVe-ETEKXBKB+VjCX3pSE5E49HHpBCz-HB3zHQLukt7w@mail.gmail.com>

Dear NCL-experts,

I am calculating potential temperature using the following ncl commands:
 theta = pot_temp(p,temp,1,False) ; [time | 365] x [level | 17] x [lat |
73] x [lon | 144]

However, my inputs *p* and *temp *have levels that are in *hPa/mbar*.
According to the documentation, this should be in *Pa (i.e., levels*100)*.

I would like to ask how do I force the variables to have levels in Pa in
the code below?




*temp_f = addfile("air.day.ltm.1981-2010.nc
<http://air.day.ltm.1981-2010.nc>", "r"   p  = temp_f->level   temp =
temp_f->air(:,{p},:,:)*









*;;Print contents of file   printVarSummary(p)
 printVarSummary(temp);;**************************************************;;Calculates
Potential Temperature;;**************************************************
theta = pot_temp(p,temp,1,False)  print(theta&level)*

I'll appreciate any help on this.

Sincerely,
Lyndz
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220824/3a202e35/attachment.htm>

From shea at ucar.edu  Wed Aug 24 09:06:56 2022
From: shea at ucar.edu (Dennis Shea)
Date: Wed, 24 Aug 2022 09:06:56 -0600
Subject: [ncl-talk] About Unit conversion handling
In-Reply-To: <CA+V+8dUVe-ETEKXBKB+VjCX3pSE5E49HHpBCz-HB3zHQLukt7w@mail.gmail.com>
References: <CA+V+8dUVe-ETEKXBKB+VjCX3pSE5E49HHpBCz-HB3zHQLukt7w@mail.gmail.com>
Message-ID: <CAOF1d_5ZOS9_M-QeOp7hxAxSKZc47jG7TnjCjAykR2jCTpqBBw@mail.gmail.com>

* p  = temp_f->level                   ; mb*
* temp = temp_f->air*

* theta = pot_temp(p*100,temp,1,False)*

*printVarSummary(theta)*

On Tue, Aug 23, 2022 at 10:28 PM Lyndz via ncl-talk <
ncl-talk at mailman.ucar.edu> wrote:

> Dear NCL-experts,
>
> I am calculating potential temperature using the following ncl commands:
>  theta = pot_temp(p,temp,1,False) ; [time | 365] x [level | 17] x [lat |
> 73] x [lon | 144]
>
> However, my inputs *p* and *temp *have levels that are in *hPa/mbar*.
> According to the documentation, this should be in *Pa (i.e., levels*100)*.
>
> I would like to ask how do I force the variables to have levels in Pa in
> the code below?
>
>
>
>
> *temp_f = addfile("air.day.ltm.1981-2010.nc
> <http://air.day.ltm.1981-2010.nc>", "r"   p  = temp_f->level   temp =
> temp_f->air(:,{p},:,:)*
>
>
>
>
>
>
>
>
>
> *;;Print contents of file   printVarSummary(p)
>  printVarSummary(temp);;**************************************************;;Calculates
> Potential Temperature;;**************************************************
> theta = pot_temp(p,temp,1,False)  print(theta&level)*
>
> I'll appreciate any help on this.
>
> Sincerely,
> Lyndz
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at mailman.ucar.edu
> List instructions, subscriber options, unsubscribe:
> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220824/a5438479/attachment.htm>

From francesco.trotta4 at unibo.it  Wed Aug 24 16:03:36 2022
From: francesco.trotta4 at unibo.it (Francesco Trotta)
Date: Wed, 24 Aug 2022 22:03:36 +0000
Subject: [ncl-talk] from ncl to python (masking data with
 shapefile_mask_data func)
Message-ID: <PR3PR01MB7938917AD67156573447EBACC6739@PR3PR01MB7938.eurprd01.prod.exchangelabs.com>

Dear ncl-experts
   I have a ncl script which is using the function 'shapefile_mask_data' to mask the bathymetry data according to the outlines of a shapefile

   I need to convert this script to python. Is it possible to use this masking utility in python? Or some alternative?

many thanks

Francesco
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220824/ddd5a28a/attachment.htm>