<div dir="ltr"><div dir="ltr">Thanks for your quickly reply.<span class="gmail-im" style="color:rgb(80,0,80)"><div dir="auto">1) the output of the command printMinMax (y,0) in line 62 is; </div><div dir="auto">(0)    MAM: precipitation  (mm/month) : min=-53.8846.   max=31.7119</div><div dir="auto"><br></div></span><div dir="auto">Additional: I have tried also to set res2@cnLevelSelectionMode  in line 100 to "ExplicitLevels" but still similar problem persist. </div><span class="gmail-im" style="color:rgb(80,0,80)"><div dir="auto"><br></div><div dir="auto">2) res2@mpLandFillColor  ="white" .I think this command has no problem because it fill the color outside the desired region i.e. Tanzania. </div><div dir="auto"><br></div><div dir="auto">Kindest regards,</div><div dir="auto">Dickson. </div></span></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Apr 20, 2020 at 8:34 PM <<a href="mailto:ncl-talk-request@ucar.edu">ncl-talk-request@ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Send ncl-talk mailing list submissions to<br>
        <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
<br>
To subscribe or unsubscribe via the World Wide Web, visit<br>
        <a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
or, via email, send a message with subject or body 'help' to<br>
        <a href="mailto:ncl-talk-request@ucar.edu" target="_blank">ncl-talk-request@ucar.edu</a><br>
<br>
You can reach the person managing the list at<br>
        <a href="mailto:ncl-talk-owner@ucar.edu" target="_blank">ncl-talk-owner@ucar.edu</a><br>
<br>
When replying, please edit your Subject line so it is more specific<br>
than "Re: Contents of ncl-talk digest..."<br>
<br>
<br>
Today's Topics:<br>
<br>
   1. Triple Integral (Prashanth Bhalachandran)<br>
   2. Getting time series from several files in ncl (zilore mumba)<br>
   3. Re: Getting time series from several files in ncl (Rick Brownrigg)<br>
   4. installing ncl on Debian Linux 10 (Francesco Trotta)<br>
   5. Re: Getting time series from several files in ncl (zilore mumba)<br>
   6. Re: Getting time series from several files in ncl (zilore mumba)<br>
   7. Re: Getting time series from several files in ncl (Rick Brownrigg)<br>
   8. extinction coefficient-calipso (zahra A)<br>
   9. Re: installing ncl on Debian Linux 10 (Gus Correa)<br>
  10. Re: extinction coefficient-calipso (Dennis Shea)<br>
  11. Re: Getting time series from several files in ncl (zilore mumba)<br>
  12. Re: Getting time series from several files in ncl (Barry Lynn)<br>
  13. Re: Getting time series from several files in ncl (Barry Lynn)<br>
  14. Re: Getting time series from several files in ncl (zilore mumba)<br>
  15. Re: Getting time series from several files in ncl (Barry Lynn)<br>
  16. Re: warning:ContourPlotInialize: no valid values in scalar<br>
      field; ContourPlot not possible: [errorno=1101] (dickson mbigi)<br>
  17. Re: warning:ContourPlotInialize: no valid values in scalar<br>
      field; ContourPlot not possible: [errorno=1101] (Alessandra Giannini)<br>
  18. Re: Getting time series from several files in ncl (zilore mumba)<br>
  19. Re: Getting time series from several files in ncl (Rashed Mahmood)<br>
<br>
<br>
----------------------------------------------------------------------<br>
<br>
Message: 1<br>
Date: Sun, 19 Apr 2020 12:23:41 -0700<br>
From: Prashanth Bhalachandran <<a href="mailto:prashanth.bhalachandran@gmail.com" target="_blank">prashanth.bhalachandran@gmail.com</a>><br>
To: ncl-talk <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>><br>
Subject: [ncl-talk] Triple Integral<br>
Message-ID:<br>
        <<a href="mailto:CALbaLSeHV3xbsoovysgjkqsN1p0jbqgrZB2te9R73xCnjRJmaw@mail.gmail.com" target="_blank">CALbaLSeHV3xbsoovysgjkqsN1p0jbqgrZB2te9R73xCnjRJmaw@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
Dear NCL'ers,<br>
I have a variable in cylindrical (r, theta, z) coordinates where the dz is<br>
varying as a function of z. I am interested in performing a volume average<br>
where the radial, as well as the vertical average, comprises dealing with<br>
elements of varying size. What is the most efficient way of performing the<br>
triple integral necessary for the volume averaging of this variable in<br>
cylindrical coordinates?<br>
<br>
I am sure many of you have dealt with this before.<br>
<br>
Sincerely,<br>
Prashanth<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20200419/4e0d10fa/attachment-0001.html" rel="noreferrer" target="_blank">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20200419/4e0d10fa/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 2<br>
Date: Sun, 19 Apr 2020 19:54:10 +0000 (UTC)<br>
From: zilore mumba <<a href="mailto:zmumba@yahoo.com" target="_blank">zmumba@yahoo.com</a>><br>
To: <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
Subject: [ncl-talk] Getting time series from several files in ncl<br>
Message-ID: <<a href="mailto:1108291697.2096087.1587326050592@mail.yahoo.com" target="_blank">1108291697.2096087.1587326050592@mail.yahoo.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
I am trying to get time series at given locations from several WRF output files, e.g. 48H forecast temperature at 850hPa. Each file contains a day's model run. I have used ncl for years by adapting existing scripts, without really understanding. One issue I have difficulty with is getting variables from a file, especially WRF output.I thought the script to get a time series at 15S, 30E would be as simple as the one below. I am able to print the varsummary, but printing x tells me that "a is not a valid file reference.I would appreciate ant assistance.begin<br>
files = systemfunc("ls -1 data/wrfout_d01_2020*") + ".nc"<br>
? a = addfiles(files,"r")<br>
<br>
? ListSetType (a, "cat")????????? ; concatenate (=default)<br>
? tc = wrf_user_getvar(a,"tc",0)<br>
? printVarSummary (tc)<br>
<br>
? x = a->T2(0,{850},{-15},{30})<br>
? asciiwrite("outfile.txt",x)<br>
end<br>
<br>
Zilore Mumba<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20200419/0d835f84/attachment-0001.html" rel="noreferrer" target="_blank">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20200419/0d835f84/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 3<br>
Date: Sun, 19 Apr 2020 14:16:59 -0600<br>
From: Rick Brownrigg <<a href="mailto:brownrig@ucar.edu" target="_blank">brownrig@ucar.edu</a>><br>
To: zilore mumba <<a href="mailto:zmumba@yahoo.com" target="_blank">zmumba@yahoo.com</a>><br>
Cc: Ncl-talk <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>><br>
Subject: Re: [ncl-talk] Getting time series from several files in ncl<br>
Message-ID:<br>
        <<a href="mailto:CAGKRhbGvgo6diU0q1X7hOkpYNUYx5X263Ny83vKWtmPLNeSSZA@mail.gmail.com" target="_blank">CAGKRhbGvgo6diU0q1X7hOkpYNUYx5X263Ny83vKWtmPLNeSSZA@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
HI Zilore,<br>
<br>
I think the issue is that addfiles() return a *list* of filevars, where<br>
addfile returns a single filevar.  Lists are indexed with [ ] notation;<br>
perhaps something like this will work:<br>
<br>
  x = a[:]->T2(0,{850},{-15},{30})<br>
<br>
Rick<br>
<br>
On Sun, Apr 19, 2020 at 1:54 PM zilore mumba via ncl-talk <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>><br>
wrote:<br>
<br>
> I am trying to get time series at given locations from several WRF output<br>
> files, e.g. 48H forecast temperature at 850hPa. Each file contains a day's<br>
> model run.<br>
> I have used ncl for years by adapting existing scripts, without really<br>
> understanding. One issue I have difficulty with is getting variables from a<br>
> file, especially WRF output.<br>
> I thought the script to get a time series at 15S, 30E would be as simple<br>
> as the one below. I am able to print the varsummary, but printing x tells<br>
> me that "a is not a valid file reference.<br>
> I would appreciate ant assistance.<br>
> begin<br>
> files = systemfunc("ls -1 data/wrfout_d01_2020*") + ".nc"<br>
>   a = addfiles(files,"r")<br>
><br>
>   ListSetType (a, "cat")          ; concatenate (=default)<br>
>   tc = wrf_user_getvar(a,"tc",0)<br>
>   printVarSummary (tc)<br>
><br>
>   x = a->T2(0,{850},{-15},{30})<br>
>   asciiwrite("outfile.txt",x)<br>
> end<br>
><br>
> Zilore Mumba<br>
> _______________________________________________<br>
> ncl-talk mailing list<br>
> <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
> List instructions, subscriber options, unsubscribe:<br>
> <a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20200419/c76852de/attachment-0001.html" rel="noreferrer" target="_blank">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20200419/c76852de/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 4<br>
Date: Sun, 19 Apr 2020 21:24:53 +0000<br>
From: Francesco Trotta <<a href="mailto:francesco.trotta4@unibo.it" target="_blank">francesco.trotta4@unibo.it</a>><br>
To: "<a href="mailto:ncl-install@ucar.edu" target="_blank">ncl-install@ucar.edu</a>" <<a href="mailto:ncl-install@ucar.edu" target="_blank">ncl-install@ucar.edu</a>>, ncl-talk<br>
        <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>><br>
Subject: [ncl-talk] installing ncl on Debian Linux 10<br>
Message-ID:<br>
        <<a href="mailto:AM0PR0102MB357227ECCAB3338F4860B51BC6D70@AM0PR0102MB3572.eurprd01.prod.exchangelabs.com" target="_blank">AM0PR0102MB357227ECCAB3338F4860B51BC6D70@AM0PR0102MB3572.eurprd01.prod.exchangelabs.com</a>><br>
<br>
Content-Type: text/plain; charset="iso-8859-1"<br>
<br>
Dear ncl team<br>
    I upgraded my operating system to Debian GNU/Linux 10 buster<br>
<br>
    This new Debian has GNU Fortran (Debian 8.3.0-6) 8.3.0 installed<br>
<br>
    Now, when I run ncl it gives me the error:<br>
   ncl: error while loading shared libraries: libgfortran.so.3: cannot open shared object file: No such file or directory<br>
<br>
   The new Debian has libgfortran.so.5<br>
<br>
   /usr/lib/x86_64-linux-gnu/libgfortran.so.5<br>
<br>
    so I upgraded ncl using to the last precompiled binary<br>
<br>
ncl_ncarg-6.6.2-Debian9.8_64bit_nodap_gnu630.tar.gz<br>
<br>
    bu I got the same error.<br>
<br>
    I can I use ncl in Debian 10?<br>
<br>
Many Thanks<br>
Francesco<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20200419/247ca69d/attachment-0001.html" rel="noreferrer" target="_blank">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20200419/247ca69d/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 5<br>
Date: Sun, 19 Apr 2020 21:27:47 +0000 (UTC)<br>
From: zilore mumba <<a href="mailto:zmumba@yahoo.com" target="_blank">zmumba@yahoo.com</a>><br>
To: Rick Brownrigg <<a href="mailto:brownrig@ucar.edu" target="_blank">brownrig@ucar.edu</a>><br>
Cc: Ncl-talk <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>><br>
Subject: Re: [ncl-talk] Getting time series from several files in ncl<br>
Message-ID: <<a href="mailto:1659568512.2117481.1587331667298@mail.yahoo.com" target="_blank">1659568512.2117481.1587331667298@mail.yahoo.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
 Thanks Rick, that is close. The error am getting now is: fatal:Number of subscripts on rhs do not match number of dimensions of aggregated cat type variable, (4) Subscripts used, (3) Subscripts expected<br>
<br>
    On Sunday, April 19, 2020, 10:17:12 PM GMT+2, Rick Brownrigg <<a href="mailto:brownrig@ucar.edu" target="_blank">brownrig@ucar.edu</a>> wrote:  <br>
<br>
 HI Zilore,<br>
I think the issue is that addfiles() return a *list* of filevars, where addfile returns a single filevar.? Lists are indexed with [ ] notation; perhaps something like this will work:<br>
? x = a[:]->T2(0,{850},{-15},{30})<br>
Rick<br>
<br>
On Sun, Apr 19, 2020 at 1:54 PM zilore mumba via ncl-talk <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>> wrote:<br>
<br>
I am trying to get time series at given locations from several WRF output files, e.g. 48H forecast temperature at 850hPa. Each file contains a day's model run. I have used ncl for years by adapting existing scripts, without really understanding. One issue I have difficulty with is getting variables from a file, especially WRF output.I thought the script to get a time series at 15S, 30E would be as simple as the one below. I am able to print the varsummary, but printing x tells me that "a is not a valid file reference.I would appreciate ant assistance.begin<br>
files = systemfunc("ls -1 data/wrfout_d01_2020*") + ".nc"<br>
? a = addfiles(files,"r")<br>
<br>
? ListSetType (a, "cat")????????? ; concatenate (=default)<br>
? tc = wrf_user_getvar(a,"tc",0)<br>
? printVarSummary (tc)<br>
<br>
? x = a->T2(0,{850},{-15},{30})<br>
? asciiwrite("outfile.txt",x)<br>
end<br>
<br>
Zilore Mumba<br>
_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20200419/ae1e1737/attachment-0001.html" rel="noreferrer" target="_blank">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20200419/ae1e1737/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 6<br>
Date: Sun, 19 Apr 2020 21:29:35 +0000 (UTC)<br>
From: zilore mumba <<a href="mailto:zmumba@yahoo.com" target="_blank">zmumba@yahoo.com</a>><br>
To: Rick Brownrigg <<a href="mailto:brownrig@ucar.edu" target="_blank">brownrig@ucar.edu</a>><br>
Cc: Ncl-talk <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>><br>
Subject: Re: [ncl-talk] Getting time series from several files in ncl<br>
Message-ID: <<a href="mailto:580311785.2109733.1587331775512@mail.yahoo.com" target="_blank">580311785.2109733.1587331775512@mail.yahoo.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
<br>
<br>
    On Sunday, April 19, 2020, 11:27:47 PM GMT+2, zilore mumba <<a href="mailto:zmumba@yahoo.com" target="_blank">zmumba@yahoo.com</a>> wrote:  <br>
<br>
  Thanks Rick, that is close. The error am getting now is: fatal:Number of subscripts on rhs do not match number of dimensions of aggregated cat type variable, (4) Subscripts used, (3) Subscripts expected<br>
<br>
    On Sunday, April 19, 2020, 10:17:12 PM GMT+2, Rick Brownrigg <<a href="mailto:brownrig@ucar.edu" target="_blank">brownrig@ucar.edu</a>> wrote:  <br>
<br>
 HI Zilore,<br>
I think the issue is that addfiles() return a *list* of filevars, where addfile returns a single filevar.? Lists are indexed with [ ] notation; perhaps something like this will work:<br>
? x = a[:]->T2(0,{850},{-15},{30})<br>
Rick<br>
<br>
On Sun, Apr 19, 2020 at 1:54 PM zilore mumba via ncl-talk <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>> wrote:<br>
<br>
I am trying to get time series at given locations from several WRF output files, e.g. 48H forecast temperature at 850hPa. Each file contains a day's model run. I have used ncl for years by adapting existing scripts, without really understanding. One issue I have difficulty with is getting variables from a file, especially WRF output.I thought the script to get a time series at 15S, 30E would be as simple as the one below. I am able to print the varsummary, but printing x tells me that "a is not a valid file reference.I would appreciate ant assistance.begin<br>
files = systemfunc("ls -1 data/wrfout_d01_2020*") + ".nc"<br>
? a = addfiles(files,"r")<br>
<br>
? ListSetType (a, "cat")????????? ; concatenate (=default)<br>
? tc = wrf_user_getvar(a,"tc",0)<br>
? printVarSummary (tc)<br>
<br>
? x = a->T2(0,{850},{-15},{30})<br>
? asciiwrite("outfile.txt",x)<br>
end<br>
<br>
Zilore Mumba<br>
_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20200419/058eb7ff/attachment-0001.html" rel="noreferrer" target="_blank">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20200419/058eb7ff/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 7<br>
Date: Sun, 19 Apr 2020 15:53:19 -0600<br>
From: Rick Brownrigg <<a href="mailto:brownrig@ucar.edu" target="_blank">brownrig@ucar.edu</a>><br>
To: zilore mumba <<a href="mailto:zmumba@yahoo.com" target="_blank">zmumba@yahoo.com</a>><br>
Cc: Ncl-talk <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>><br>
Subject: Re: [ncl-talk] Getting time series from several files in ncl<br>
Message-ID:<br>
        <CAGKRhbEtdtm7X4cd=uLuVXAZeUMLtCCZn6062s8dvg8Msi_=<a href="mailto:tw@mail.gmail.com" target="_blank">tw@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
Perhaps do a printVarSummary(a[0]->T2)  ; it sounds like its not the shape<br>
you think it might/should be; i.e. 3D vs 4D<br>
<br>
<br>
<br>
On Sun, Apr 19, 2020 at 3:29 PM zilore mumba <<a href="mailto:zmumba@yahoo.com" target="_blank">zmumba@yahoo.com</a>> wrote:<br>
<br>
><br>
><br>
> On Sunday, April 19, 2020, 11:27:47 PM GMT+2, zilore mumba <<br>
> <a href="mailto:zmumba@yahoo.com" target="_blank">zmumba@yahoo.com</a>> wrote:<br>
><br>
><br>
> Thanks Rick, that is close. The error am getting now is:<br>
> fatal:Number of subscripts on rhs do not match number of dimensions of<br>
> aggregated cat type variable, (4) Subscripts used, (3) Subscripts expected<br>
><br>
> On Sunday, April 19, 2020, 10:17:12 PM GMT+2, Rick Brownrigg <<br>
> <a href="mailto:brownrig@ucar.edu" target="_blank">brownrig@ucar.edu</a>> wrote:<br>
><br>
><br>
> HI Zilore,<br>
><br>
> I think the issue is that addfiles() return a *list* of filevars, where<br>
> addfile returns a single filevar.  Lists are indexed with [ ] notation;<br>
> perhaps something like this will work:<br>
><br>
>   x = a[:]->T2(0,{850},{-15},{30})<br>
><br>
> Rick<br>
><br>
> On Sun, Apr 19, 2020 at 1:54 PM zilore mumba via ncl-talk <<br>
> <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>> wrote:<br>
><br>
> I am trying to get time series at given locations from several WRF output<br>
> files, e.g. 48H forecast temperature at 850hPa. Each file contains a day's<br>
> model run.<br>
> I have used ncl for years by adapting existing scripts, without really<br>
> understanding. One issue I have difficulty with is getting variables from a<br>
> file, especially WRF output.<br>
> I thought the script to get a time series at 15S, 30E would be as simple<br>
> as the one below. I am able to print the varsummary, but printing x tells<br>
> me that "a is not a valid file reference.<br>
> I would appreciate ant assistance.<br>
> begin<br>
> files = systemfunc("ls -1 data/wrfout_d01_2020*") + ".nc"<br>
>   a = addfiles(files,"r")<br>
><br>
>   ListSetType (a, "cat")          ; concatenate (=default)<br>
>   tc = wrf_user_getvar(a,"tc",0)<br>
>   printVarSummary (tc)<br>
><br>
>   x = a->T2(0,{850},{-15},{30})<br>
>   asciiwrite("outfile.txt",x)<br>
> end<br>
><br>
> Zilore Mumba<br>
> _______________________________________________<br>
> ncl-talk mailing list<br>
> <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
> List instructions, subscriber options, unsubscribe:<br>
> <a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
><br>
><br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20200419/1458e762/attachment-0001.html" rel="noreferrer" target="_blank">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20200419/1458e762/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 8<br>
Date: Mon, 20 Apr 2020 02:56:00 +0430<br>
From: zahra A <<a href="mailto:zahra93578@gmail.com" target="_blank">zahra93578@gmail.com</a>><br>
To: <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
Subject: [ncl-talk] extinction coefficient-calipso<br>
Message-ID:<br>
        <CAJfT7CaF=HB6uE7_ziUKwyxOj5Ok9KUa2htH_JKhaM4i1KY=<a href="mailto:3A@mail.gmail.com" target="_blank">3A@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
Hi<br>
I am trying to derive the diagram of extinction coefficient and particulate<br>
depolarization ratio from calipso level 2 layer and profile products for<br>
aerosols. I would appreciate if anybody can tell me how can I find any<br>
example code in this regard.<br>
Thanks<br>
Zahra<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20200420/dd83e513/attachment-0001.html" rel="noreferrer" target="_blank">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20200420/dd83e513/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 9<br>
Date: Sun, 19 Apr 2020 21:57:51 -0400<br>
From: Gus Correa <<a href="mailto:gus@ldeo.columbia.edu" target="_blank">gus@ldeo.columbia.edu</a>><br>
To: NCL Talk <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>><br>
Subject: Re: [ncl-talk] installing ncl on Debian Linux 10<br>
Message-ID:<br>
        <CA+TEG=XQCzsNEc9W=AvdFdKWH=XBfyMLtW=<a href="mailto:k2mKQtM_D98iQNw@mail.gmail.com" target="_blank">k2mKQtM_D98iQNw@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
Hi Francesco<br>
<br>
Have you checked if your Debian distribution has a "compat" (compatibility)<br>
package for gfortran/libgfortran<br>
between libgfortran.5 and libgfortran.3 ?<br>
Apt-get may have a way to check if it is available.<br>
I don't use Debian, but CentOS-7, which has an older kernel and older<br>
libraries.<br>
However, I have libgfortran.3 along with a "compat" package that provides<br>
compatibility with libgfortran.1, for example.<br>
<br>
An alternative is to install NCL using Anaconda Python.<br>
This is currently the recommended method (argh!):<br>
<a href="https://www.ncl.ucar.edu/Download/" rel="noreferrer" target="_blank">https://www.ncl.ucar.edu/Download/</a><br>
That is self-contained, doesn't depend on the system libraries.<br>
In my experience it has other problems, particularly a sensitive dependency<br>
on other Anaconda Python packages<br>
that change by the hour, and don't give a damn to backward compatibility,<br>
which can easily break functionality.<br>
So, if you install this way, be careful with any conda packages updates, or<br>
simply don't update.<br>
<br>
I hope this helps,<br>
Gus Correa<br>
<br>
On Sun, Apr 19, 2020 at 5:25 PM Francesco Trotta via ncl-talk <<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>> wrote:<br>
<br>
> Dear ncl team<br>
>     I upgraded my operating system to Debian GNU/Linux 10 buster<br>
><br>
>     This new Debian has GNU Fortran (Debian 8.3.0-6) 8.3.0 installed<br>
><br>
>     Now, when I run ncl it gives me the error:<br>
>    ncl: error while loading shared libraries: libgfortran.so.3: cannot<br>
> open shared object file: No such file or directory<br>
><br>
>    The new Debian has libgfortran.so.5<br>
><br>
>    /usr/lib/x86_64-linux-gnu/libgfortran.so.5<br>
><br>
>     so I upgraded ncl using to the last precompiled binary<br>
><br>
> ncl_ncarg-6.6.2-Debian9.8_64bit_nodap_gnu630.tar.gz<br>
><br>
>     bu I got the same error.<br>
><br>
>     I can I use ncl in Debian 10?<br>
><br>
> Many Thanks<br>
> Francesco<br>
> _______________________________________________<br>
> ncl-talk mailing list<br>
> <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
> List instructions, subscriber options, unsubscribe:<br>
> <a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20200419/e8325ba6/attachment-0001.html" rel="noreferrer" target="_blank">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20200419/e8325ba6/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 10<br>
Date: Sun, 19 Apr 2020 20:19:02 -0600<br>
From: Dennis Shea <<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>><br>
To: zahra A <<a href="mailto:zahra93578@gmail.com" target="_blank">zahra93578@gmail.com</a>><br>
Cc: Ncl-talk <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>><br>
Subject: Re: [ncl-talk] extinction coefficient-calipso<br>
Message-ID:<br>
        <<a href="mailto:CAOF1d_5gHOwbVjpyqzhu-UiGY9HsqfDWfcaEWEpq1Y1Q1MaSQQ@mail.gmail.com" target="_blank">CAOF1d_5gHOwbVjpyqzhu-UiGY9HsqfDWfcaEWEpq1Y1Q1MaSQQ@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
There are no NCL codes to do what you want.<br>
----<br>
I suggest contacting the Calipso team:  *<a href="https://www-calipso.larc.nasa.gov/*" rel="noreferrer" target="_blank">https://www-calipso.larc.nasa.gov/*</a><br>
<<a href="https://www-calipso.larc.nasa.gov/" rel="noreferrer" target="_blank">https://www-calipso.larc.nasa.gov/</a>><br>
<br>
The following may be useful to you:<br>
<br>
<br>
*CALIOP Algorithm Theoretical Basis Document Part 4: Extinction Retrieval<br>
Algorithms*<br>
*<a href="https://www-calipso.larc.nasa.gov/resources/pdfs/PC-SCI-202_Part4_v1.0.pdf*" rel="noreferrer" target="_blank">https://www-calipso.larc.nasa.gov/resources/pdfs/PC-SCI-202_Part4_v1.0.pdf*</a><br>
<<a href="https://www-calipso.larc.nasa.gov/resources/pdfs/PC-SCI-202_Part4_v1.0.pdf" rel="noreferrer" target="_blank">https://www-calipso.larc.nasa.gov/resources/pdfs/PC-SCI-202_Part4_v1.0.pdf</a>><br>
<br>
Good Luck<br>
<br>
On Sun, Apr 19, 2020 at 4:26 PM zahra A via ncl-talk <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>><br>
wrote:<br>
<br>
> Hi<br>
> I am trying to derive the diagram of extinction coefficient and<br>
> particulate depolarization ratio from calipso level 2 layer and profile<br>
> products for aerosols. I would appreciate if anybody can tell me how can I<br>
> find any example code in this regard.<br>
> Thanks<br>
> Zahra<br>
> _______________________________________________<br>
> ncl-talk mailing list<br>
> <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
> List instructions, subscriber options, unsubscribe:<br>
> <a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20200419/02a7ce22/attachment-0001.html" rel="noreferrer" target="_blank">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20200419/02a7ce22/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 11<br>
Date: Mon, 20 Apr 2020 05:17:46 +0000 (UTC)<br>
From: zilore mumba <<a href="mailto:zmumba@yahoo.com" target="_blank">zmumba@yahoo.com</a>><br>
To: Rick Brownrigg <<a href="mailto:brownrig@ucar.edu" target="_blank">brownrig@ucar.edu</a>><br>
Cc: Ncl-talk <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>><br>
Subject: Re: [ncl-talk] Getting time series from several files in ncl<br>
Message-ID: <<a href="mailto:1685826303.2207658.1587359866358@mail.yahoo.com" target="_blank">1685826303.2207658.1587359866358@mail.yahoo.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
 Sorry Rick I put in a surface variable T2 which only 3 coordinates, so I put "x = a[:]->T2(0,{-15},{30})". Now I am getting error of coordinates west-east and south-north not having? associated coordinate variables. I have seen this error before. I will find out how to get out of it <br>
<br>
    On Sunday, April 19, 2020, 11:53:33 PM GMT+2, Rick Brownrigg <<a href="mailto:brownrig@ucar.edu" target="_blank">brownrig@ucar.edu</a>> wrote:  <br>
<br>
 Perhaps do a printVarSummary(a[0]->T2)? ; it sounds like its not the shape you think it might/should be; i.e. 3D vs 4D<br>
<br>
<br>
On Sun, Apr 19, 2020 at 3:29 PM zilore mumba <<a href="mailto:zmumba@yahoo.com" target="_blank">zmumba@yahoo.com</a>> wrote:<br>
<br>
<br>
<br>
    On Sunday, April 19, 2020, 11:27:47 PM GMT+2, zilore mumba <<a href="mailto:zmumba@yahoo.com" target="_blank">zmumba@yahoo.com</a>> wrote:  <br>
<br>
  Thanks Rick, that is close. The error am getting now is: fatal:Number of subscripts on rhs do not match number of dimensions of aggregated cat type variable, (4) Subscripts used, (3) Subscripts expected<br>
<br>
    On Sunday, April 19, 2020, 10:17:12 PM GMT+2, Rick Brownrigg <<a href="mailto:brownrig@ucar.edu" target="_blank">brownrig@ucar.edu</a>> wrote:  <br>
<br>
 HI Zilore,<br>
I think the issue is that addfiles() return a *list* of filevars, where addfile returns a single filevar.? Lists are indexed with [ ] notation; perhaps something like this will work:<br>
? x = a[:]->T2(0,{850},{-15},{30})<br>
Rick<br>
<br>
On Sun, Apr 19, 2020 at 1:54 PM zilore mumba via ncl-talk <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>> wrote:<br>
<br>
I am trying to get time series at given locations from several WRF output files, e.g. 48H forecast temperature at 850hPa. Each file contains a day's model run. I have used ncl for years by adapting existing scripts, without really understanding. One issue I have difficulty with is getting variables from a file, especially WRF output.I thought the script to get a time series at 15S, 30E would be as simple as the one below. I am able to print the varsummary, but printing x tells me that "a is not a valid file reference.I would appreciate ant assistance.begin<br>
files = systemfunc("ls -1 data/wrfout_d01_2020*") + ".nc"<br>
? a = addfiles(files,"r")<br>
<br>
? ListSetType (a, "cat")????????? ; concatenate (=default)<br>
? tc = wrf_user_getvar(a,"tc",0)<br>
? printVarSummary (tc)<br>
<br>
? x = a->T2(0,{850},{-15},{30})<br>
? asciiwrite("outfile.txt",x)<br>
end<br>
<br>
Zilore Mumba<br>
_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br>
<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20200420/9abe6e52/attachment-0001.html" rel="noreferrer" target="_blank">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20200420/9abe6e52/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 12<br>
Date: Mon, 20 Apr 2020 09:00:28 +0300<br>
From: Barry Lynn <<a href="mailto:barry.h.lynn@gmail.com" target="_blank">barry.h.lynn@gmail.com</a>><br>
To: zilore mumba <<a href="mailto:zmumba@yahoo.com" target="_blank">zmumba@yahoo.com</a>><br>
Cc: Ncl-talk <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>><br>
Subject: Re: [ncl-talk] Getting time series from several files in ncl<br>
Message-ID:<br>
        <<a href="mailto:CAN7Bv24BAL4LNp_6y6edFwt%2BSg-vr%2BeaLgvTfO0wtmnvteyn6w@mail.gmail.com" target="_blank">CAN7Bv24BAL4LNp_6y6edFwt+Sg-vr+eaLgvTfO0wtmnvteyn6w@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
Zilore:<br>
<br>
T2 is a three dimensional variable.  Also, T2 is the two meter temperature.<br>
<br>
You want tc or tk (see:<br>
<a href="https://www.ncl.ucar.edu/Document/Functions/WRF_arw/wrf_user_getvar.shtml" rel="noreferrer" target="_blank">https://www.ncl.ucar.edu/Document/Functions/WRF_arw/wrf_user_getvar.shtml</a>).<br>
<br>
Could you please explain how you are interpolating to 850 mb?  There is a<br>
WRF program to do this, and I am not sure how the indexing you are using<br>
works. How are you indexing the 48 time slot?  WRF has a program to list<br>
times, and you can match the time you want to the listed times to find the<br>
index.<br>
<br>
Barry<br>
<br>
On Mon, Apr 20, 2020 at 8:17 AM zilore mumba via ncl-talk <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>><br>
wrote:<br>
<br>
> Sorry Rick I put in a surface variable T2 which only 3 coordinates, so I<br>
> put "x = a[:]->T2(0,{-15},{30})". Now I am getting error of coordinates<br>
> west-east and south-north not having  associated coordinate variables. I<br>
> have seen this error before. I will find out how to get out of it<br>
><br>
> On Sunday, April 19, 2020, 11:53:33 PM GMT+2, Rick Brownrigg <<br>
> <a href="mailto:brownrig@ucar.edu" target="_blank">brownrig@ucar.edu</a>> wrote:<br>
><br>
><br>
> Perhaps do a printVarSummary(a[0]->T2)  ; it sounds like its not the shape<br>
> you think it might/should be; i.e. 3D vs 4D<br>
><br>
><br>
><br>
> On Sun, Apr 19, 2020 at 3:29 PM zilore mumba <<a href="mailto:zmumba@yahoo.com" target="_blank">zmumba@yahoo.com</a>> wrote:<br>
><br>
><br>
><br>
> On Sunday, April 19, 2020, 11:27:47 PM GMT+2, zilore mumba <<br>
> <a href="mailto:zmumba@yahoo.com" target="_blank">zmumba@yahoo.com</a>> wrote:<br>
><br>
><br>
> Thanks Rick, that is close. The error am getting now is:<br>
> fatal:Number of subscripts on rhs do not match number of dimensions of<br>
> aggregated cat type variable, (4) Subscripts used, (3) Subscripts expected<br>
><br>
> On Sunday, April 19, 2020, 10:17:12 PM GMT+2, Rick Brownrigg <<br>
> <a href="mailto:brownrig@ucar.edu" target="_blank">brownrig@ucar.edu</a>> wrote:<br>
><br>
><br>
> HI Zilore,<br>
><br>
> I think the issue is that addfiles() return a *list* of filevars, where<br>
> addfile returns a single filevar.  Lists are indexed with [ ] notation;<br>
> perhaps something like this will work:<br>
><br>
>   x = a[:]->T2(0,{850},{-15},{30})<br>
><br>
> Rick<br>
><br>
> On Sun, Apr 19, 2020 at 1:54 PM zilore mumba via ncl-talk <<br>
> <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>> wrote:<br>
><br>
> I am trying to get time series at given locations from several WRF output<br>
> files, e.g. 48H forecast temperature at 850hPa. Each file contains a day's<br>
> model run.<br>
> I have used ncl for years by adapting existing scripts, without really<br>
> understanding. One issue I have difficulty with is getting variables from a<br>
> file, especially WRF output.<br>
> I thought the script to get a time series at 15S, 30E would be as simple<br>
> as the one below. I am able to print the varsummary, but printing x tells<br>
> me that "a is not a valid file reference.<br>
> I would appreciate ant assistance.<br>
> begin<br>
> files = systemfunc("ls -1 data/wrfout_d01_2020*") + ".nc"<br>
>   a = addfiles(files,"r")<br>
><br>
>   ListSetType (a, "cat")          ; concatenate (=default)<br>
>   tc = wrf_user_getvar(a,"tc",0)<br>
>   printVarSummary (tc)<br>
><br>
>   x = a->T2(0,{850},{-15},{30})<br>
>   asciiwrite("outfile.txt",x)<br>
> end<br>
><br>
> Zilore Mumba<br>
> _______________________________________________<br>
> ncl-talk mailing list<br>
> <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
> List instructions, subscriber options, unsubscribe:<br>
> <a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
><br>
> _______________________________________________<br>
> ncl-talk mailing list<br>
> <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
> List instructions, subscriber options, unsubscribe:<br>
> <a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br>
<br>
<br>
-- <br>
Barry H. Lynn, Ph.D<br>
Senior Associate Scientist, Lecturer,<br>
The Institute of the Earth Science,<br>
The Hebrew University of Jerusalem,<br>
Givat Ram, Jerusalem 91904, Israel<br>
Tel: 972 547 231 170<br>
Fax: (972)-25662581<br>
<br>
C.E.O, Weather It Is, LTD<br>
Weather and Climate Focus<br>
<a href="http://weather-it-is.com" rel="noreferrer" target="_blank">http://weather-it-is.com</a><br>
Jerusalem, Israel<br>
Local: 02 930 9525<br>
Cell: 054 7 231 170<br>
Int-IS: x972 2 930 9525<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20200420/3712537e/attachment-0001.html" rel="noreferrer" target="_blank">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20200420/3712537e/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 13<br>
Date: Mon, 20 Apr 2020 09:05:50 +0300<br>
From: Barry Lynn <<a href="mailto:barry.h.lynn@gmail.com" target="_blank">barry.h.lynn@gmail.com</a>><br>
To: zilore mumba <<a href="mailto:zmumba@yahoo.com" target="_blank">zmumba@yahoo.com</a>><br>
Cc: Ncl-talk <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>><br>
Subject: Re: [ncl-talk] Getting time series from several files in ncl<br>
Message-ID:<br>
        <CAN7Bv273mEvwVu=<a href="mailto:4QTXyoMEbLN16fn6aHe6GKYQfuXuMgBQbjw@mail.gmail.com" target="_blank">4QTXyoMEbLN16fn6aHe6GKYQfuXuMgBQbjw@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
Just to clarify: T2 is 2 dimensional.   It is written every time step as<br>
such.  Tk and Tc are three dimensional.<br>
<br>
Barry<br>
<br>
On Mon, Apr 20, 2020 at 9:00 AM Barry Lynn <<a href="mailto:barry.h.lynn@gmail.com" target="_blank">barry.h.lynn@gmail.com</a>> wrote:<br>
<br>
> Zilore:<br>
><br>
> T2 is a three dimensional variable.  Also, T2 is the two meter<br>
> temperature.<br>
><br>
> You want tc or tk (see:<br>
> <a href="https://www.ncl.ucar.edu/Document/Functions/WRF_arw/wrf_user_getvar.shtml" rel="noreferrer" target="_blank">https://www.ncl.ucar.edu/Document/Functions/WRF_arw/wrf_user_getvar.shtml</a><br>
> ).<br>
><br>
> Could you please explain how you are interpolating to 850 mb?  There is a<br>
> WRF program to do this, and I am not sure how the indexing you are using<br>
> works. How are you indexing the 48 time slot?  WRF has a program to list<br>
> times, and you can match the time you want to the listed times to find the<br>
> index.<br>
><br>
> Barry<br>
><br>
> On Mon, Apr 20, 2020 at 8:17 AM zilore mumba via ncl-talk <<br>
> <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>> wrote:<br>
><br>
>> Sorry Rick I put in a surface variable T2 which only 3 coordinates, so I<br>
>> put "x = a[:]->T2(0,{-15},{30})". Now I am getting error of coordinates<br>
>> west-east and south-north not having  associated coordinate variables. I<br>
>> have seen this error before. I will find out how to get out of it<br>
>><br>
>> On Sunday, April 19, 2020, 11:53:33 PM GMT+2, Rick Brownrigg <<br>
>> <a href="mailto:brownrig@ucar.edu" target="_blank">brownrig@ucar.edu</a>> wrote:<br>
>><br>
>><br>
>> Perhaps do a printVarSummary(a[0]->T2)  ; it sounds like its not the<br>
>> shape you think it might/should be; i.e. 3D vs 4D<br>
>><br>
>><br>
>><br>
>> On Sun, Apr 19, 2020 at 3:29 PM zilore mumba <<a href="mailto:zmumba@yahoo.com" target="_blank">zmumba@yahoo.com</a>> wrote:<br>
>><br>
>><br>
>><br>
>> On Sunday, April 19, 2020, 11:27:47 PM GMT+2, zilore mumba <<br>
>> <a href="mailto:zmumba@yahoo.com" target="_blank">zmumba@yahoo.com</a>> wrote:<br>
>><br>
>><br>
>> Thanks Rick, that is close. The error am getting now is:<br>
>> fatal:Number of subscripts on rhs do not match number of dimensions of<br>
>> aggregated cat type variable, (4) Subscripts used, (3) Subscripts expected<br>
>><br>
>> On Sunday, April 19, 2020, 10:17:12 PM GMT+2, Rick Brownrigg <<br>
>> <a href="mailto:brownrig@ucar.edu" target="_blank">brownrig@ucar.edu</a>> wrote:<br>
>><br>
>><br>
>> HI Zilore,<br>
>><br>
>> I think the issue is that addfiles() return a *list* of filevars, where<br>
>> addfile returns a single filevar.  Lists are indexed with [ ] notation;<br>
>> perhaps something like this will work:<br>
>><br>
>>   x = a[:]->T2(0,{850},{-15},{30})<br>
>><br>
>> Rick<br>
>><br>
>> On Sun, Apr 19, 2020 at 1:54 PM zilore mumba via ncl-talk <<br>
>> <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>> wrote:<br>
>><br>
>> I am trying to get time series at given locations from several WRF output<br>
>> files, e.g. 48H forecast temperature at 850hPa. Each file contains a day's<br>
>> model run.<br>
>> I have used ncl for years by adapting existing scripts, without really<br>
>> understanding. One issue I have difficulty with is getting variables from a<br>
>> file, especially WRF output.<br>
>> I thought the script to get a time series at 15S, 30E would be as simple<br>
>> as the one below. I am able to print the varsummary, but printing x tells<br>
>> me that "a is not a valid file reference.<br>
>> I would appreciate ant assistance.<br>
>> begin<br>
>> files = systemfunc("ls -1 data/wrfout_d01_2020*") + ".nc"<br>
>>   a = addfiles(files,"r")<br>
>><br>
>>   ListSetType (a, "cat")          ; concatenate (=default)<br>
>>   tc = wrf_user_getvar(a,"tc",0)<br>
>>   printVarSummary (tc)<br>
>><br>
>>   x = a->T2(0,{850},{-15},{30})<br>
>>   asciiwrite("outfile.txt",x)<br>
>> end<br>
>><br>
>> Zilore Mumba<br>
>> _______________________________________________<br>
>> ncl-talk mailing list<br>
>> <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
>> List instructions, subscriber options, unsubscribe:<br>
>> <a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
>><br>
>> _______________________________________________<br>
>> ncl-talk mailing list<br>
>> <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
>> List instructions, subscriber options, unsubscribe:<br>
>> <a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
><br>
><br>
><br>
> --<br>
> Barry H. Lynn, Ph.D<br>
> Senior Associate Scientist, Lecturer,<br>
> The Institute of the Earth Science,<br>
> The Hebrew University of Jerusalem,<br>
> Givat Ram, Jerusalem 91904, Israel<br>
> Tel: 972 547 231 170<br>
> Fax: (972)-25662581<br>
><br>
> C.E.O, Weather It Is, LTD<br>
> Weather and Climate Focus<br>
> <a href="http://weather-it-is.com" rel="noreferrer" target="_blank">http://weather-it-is.com</a><br>
> Jerusalem, Israel<br>
> Local: 02 930 9525<br>
> Cell: 054 7 231 170<br>
> Int-IS: x972 2 930 9525<br>
><br>
><br>
<br>
-- <br>
Barry H. Lynn, Ph.D<br>
Senior Associate Scientist, Lecturer,<br>
The Institute of the Earth Science,<br>
The Hebrew University of Jerusalem,<br>
Givat Ram, Jerusalem 91904, Israel<br>
Tel: 972 547 231 170<br>
Fax: (972)-25662581<br>
<br>
C.E.O, Weather It Is, LTD<br>
Weather and Climate Focus<br>
<a href="http://weather-it-is.com" rel="noreferrer" target="_blank">http://weather-it-is.com</a><br>
Jerusalem, Israel<br>
Local: 02 930 9525<br>
Cell: 054 7 231 170<br>
Int-IS: x972 2 930 9525<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20200420/f8a7122c/attachment-0001.html" rel="noreferrer" target="_blank">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20200420/f8a7122c/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 14<br>
Date: Mon, 20 Apr 2020 08:55:31 +0000 (UTC)<br>
From: zilore mumba <<a href="mailto:zmumba@yahoo.com" target="_blank">zmumba@yahoo.com</a>><br>
To: Barry Lynn <<a href="mailto:barry.h.lynn@gmail.com" target="_blank">barry.h.lynn@gmail.com</a>><br>
Cc: Ncl-talk <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>><br>
Subject: Re: [ncl-talk] Getting time series from several files in ncl<br>
Message-ID: <<a href="mailto:318273061.2228766.1587372931388@mail.yahoo.com" target="_blank">318273061.2228766.1587372931388@mail.yahoo.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
 Barry,My intention is to have RH and U, V time series. That's what I had in mind when I put 48. T2 only has time, lon ,lat. I have several daily output files, and I want an n-day series of 48 hour forecasts of T2<br>
<br>
    On Monday, April 20, 2020, 8:00:41 AM GMT+2, Barry Lynn <<a href="mailto:barry.h.lynn@gmail.com" target="_blank">barry.h.lynn@gmail.com</a>> wrote:  <br>
<br>
 Zilore:<br>
T2 is a three dimensional variable.? Also, T2 is the two meter temperature.?<br>
You want tc or tk (see:?<a href="https://www.ncl.ucar.edu/Document/Functions/WRF_arw/wrf_user_getvar.shtml" rel="noreferrer" target="_blank">https://www.ncl.ucar.edu/Document/Functions/WRF_arw/wrf_user_getvar.shtml</a>).<br>
Could you please explain how you are interpolating to 850 mb?? There is a WRF program to do this, and I am not sure how the indexing you are using works. How are you indexing the 48 time slot?? WRF has a program to list times, and you can match the time you want to the listed times to find the index.<br>
Barry<br>
On Mon, Apr 20, 2020 at 8:17 AM zilore mumba via ncl-talk <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>> wrote:<br>
<br>
 Sorry Rick I put in a surface variable T2 which only 3 coordinates, so I put "x = a[:]->T2(0,{-15},{30})". Now I am getting error of coordinates west-east and south-north not having? associated coordinate variables. I have seen this error before. I will find out how to get out of it <br>
<br>
    On Sunday, April 19, 2020, 11:53:33 PM GMT+2, Rick Brownrigg <<a href="mailto:brownrig@ucar.edu" target="_blank">brownrig@ucar.edu</a>> wrote:  <br>
<br>
 Perhaps do a printVarSummary(a[0]->T2)? ; it sounds like its not the shape you think it might/should be; i.e. 3D vs 4D<br>
<br>
<br>
On Sun, Apr 19, 2020 at 3:29 PM zilore mumba <<a href="mailto:zmumba@yahoo.com" target="_blank">zmumba@yahoo.com</a>> wrote:<br>
<br>
<br>
<br>
    On Sunday, April 19, 2020, 11:27:47 PM GMT+2, zilore mumba <<a href="mailto:zmumba@yahoo.com" target="_blank">zmumba@yahoo.com</a>> wrote:  <br>
<br>
  Thanks Rick, that is close. The error am getting now is: fatal:Number of subscripts on rhs do not match number of dimensions of aggregated cat type variable, (4) Subscripts used, (3) Subscripts expected<br>
<br>
    On Sunday, April 19, 2020, 10:17:12 PM GMT+2, Rick Brownrigg <<a href="mailto:brownrig@ucar.edu" target="_blank">brownrig@ucar.edu</a>> wrote:  <br>
<br>
 HI Zilore,<br>
I think the issue is that addfiles() return a *list* of filevars, where addfile returns a single filevar.? Lists are indexed with [ ] notation; perhaps something like this will work:<br>
? x = a[:]->T2(0,{850},{-15},{30})<br>
Rick<br>
<br>
On Sun, Apr 19, 2020 at 1:54 PM zilore mumba via ncl-talk <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>> wrote:<br>
<br>
I am trying to get time series at given locations from several WRF output files, e.g. 48H forecast temperature at 850hPa. Each file contains a day's model run. I have used ncl for years by adapting existing scripts, without really understanding. One issue I have difficulty with is getting variables from a file, especially WRF output.I thought the script to get a time series at 15S, 30E would be as simple as the one below. I am able to print the varsummary, but printing x tells me that "a is not a valid file reference.I would appreciate ant assistance.begin<br>
files = systemfunc("ls -1 data/wrfout_d01_2020*") + ".nc"<br>
? a = addfiles(files,"r")<br>
<br>
? ListSetType (a, "cat")????????? ; concatenate (=default)<br>
? tc = wrf_user_getvar(a,"tc",0)<br>
? printVarSummary (tc)<br>
<br>
? x = a->T2(0,{850},{-15},{30})<br>
? asciiwrite("outfile.txt",x)<br>
end<br>
<br>
Zilore Mumba<br>
_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br>
  _______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br>
<br>
-- <br>
Barry H. Lynn, Ph.DSenior Associate Scientist, Lecturer,The Institute of the Earth Science,?<br>
The Hebrew University of Jerusalem,?<br>
Givat Ram, Jerusalem 91904, Israel?<br>
Tel: 972 547 231 170<br>
Fax: (972)-25662581<br>
C.E.O, Weather It Is, LTD<br>
Weather and Climate Focus<br>
<a href="http://weather-it-is.com" rel="noreferrer" target="_blank">http://weather-it-is.com</a><br>
Jerusalem, Israel<br>
Local: 02 930 9525<br>
Cell: 054 7 231 170<br>
Int-IS: x972 2 930 9525<br>
<br>
<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20200420/3367f406/attachment-0001.html" rel="noreferrer" target="_blank">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20200420/3367f406/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 15<br>
Date: Mon, 20 Apr 2020 13:02:16 +0300<br>
From: Barry Lynn <<a href="mailto:barry.h.lynn@gmail.com" target="_blank">barry.h.lynn@gmail.com</a>><br>
To: zilore mumba <<a href="mailto:zmumba@yahoo.com" target="_blank">zmumba@yahoo.com</a>><br>
Cc: Ncl-talk <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>><br>
Subject: Re: [ncl-talk] Getting time series from several files in ncl<br>
Message-ID:<br>
        <CAN7Bv24Kz1STaSW9mH5=61XX5yjUtnYaZgUPX=<a href="mailto:gMgXeGqniXNg@mail.gmail.com" target="_blank">gMgXeGqniXNg@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
If you know when these values occur within the time series, then you can<br>
index them quite easily.<br>
<br>
Why don't you use wrf_user_getvar and then the wrf interpolation program.<br>
<br>
If you want them at a specific point, I think you know how to do this.<br>
There is also a get location program<br>
<br>
<a href="https://www.ncl.ucar.edu/Document/Functions/WRF_arw/wrf_user_ll_to_ij.shtml" rel="noreferrer" target="_blank">https://www.ncl.ucar.edu/Document/Functions/WRF_arw/wrf_user_ll_to_ij.shtml</a><br>
<br>
Barry<br>
<br>
On Mon, Apr 20, 2020 at 11:55 AM zilore mumba <<a href="mailto:zmumba@yahoo.com" target="_blank">zmumba@yahoo.com</a>> wrote:<br>
<br>
> Barry,<br>
> My intention is to have RH and U, V time series. That's what I had in mind<br>
> when I put 48. T2 only has time, lon ,lat. I have several daily output<br>
> files, and I want an n-day series of 48 hour forecasts of T2<br>
><br>
> On Monday, April 20, 2020, 8:00:41 AM GMT+2, Barry Lynn <<br>
> <a href="mailto:barry.h.lynn@gmail.com" target="_blank">barry.h.lynn@gmail.com</a>> wrote:<br>
><br>
><br>
> Zilore:<br>
><br>
> T2 is a three dimensional variable.  Also, T2 is the two meter<br>
> temperature.<br>
><br>
> You want tc or tk (see:<br>
> <a href="https://www.ncl.ucar.edu/Document/Functions/WRF_arw/wrf_user_getvar.shtml" rel="noreferrer" target="_blank">https://www.ncl.ucar.edu/Document/Functions/WRF_arw/wrf_user_getvar.shtml</a><br>
> ).<br>
><br>
> Could you please explain how you are interpolating to 850 mb?  There is a<br>
> WRF program to do this, and I am not sure how the indexing you are using<br>
> works. How are you indexing the 48 time slot?  WRF has a program to list<br>
> times, and you can match the time you want to the listed times to find the<br>
> index.<br>
><br>
> Barry<br>
><br>
> On Mon, Apr 20, 2020 at 8:17 AM zilore mumba via ncl-talk <<br>
> <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>> wrote:<br>
><br>
> Sorry Rick I put in a surface variable T2 which only 3 coordinates, so I<br>
> put "x = a[:]->T2(0,{-15},{30})". Now I am getting error of coordinates<br>
> west-east and south-north not having  associated coordinate variables. I<br>
> have seen this error before. I will find out how to get out of it<br>
><br>
> On Sunday, April 19, 2020, 11:53:33 PM GMT+2, Rick Brownrigg <<br>
> <a href="mailto:brownrig@ucar.edu" target="_blank">brownrig@ucar.edu</a>> wrote:<br>
><br>
><br>
> Perhaps do a printVarSummary(a[0]->T2)  ; it sounds like its not the shape<br>
> you think it might/should be; i.e. 3D vs 4D<br>
><br>
><br>
><br>
> On Sun, Apr 19, 2020 at 3:29 PM zilore mumba <<a href="mailto:zmumba@yahoo.com" target="_blank">zmumba@yahoo.com</a>> wrote:<br>
><br>
><br>
><br>
> On Sunday, April 19, 2020, 11:27:47 PM GMT+2, zilore mumba <<br>
> <a href="mailto:zmumba@yahoo.com" target="_blank">zmumba@yahoo.com</a>> wrote:<br>
><br>
><br>
> Thanks Rick, that is close. The error am getting now is:<br>
> fatal:Number of subscripts on rhs do not match number of dimensions of<br>
> aggregated cat type variable, (4) Subscripts used, (3) Subscripts expected<br>
><br>
> On Sunday, April 19, 2020, 10:17:12 PM GMT+2, Rick Brownrigg <<br>
> <a href="mailto:brownrig@ucar.edu" target="_blank">brownrig@ucar.edu</a>> wrote:<br>
><br>
><br>
> HI Zilore,<br>
><br>
> I think the issue is that addfiles() return a *list* of filevars, where<br>
> addfile returns a single filevar.  Lists are indexed with [ ] notation;<br>
> perhaps something like this will work:<br>
><br>
>   x = a[:]->T2(0,{850},{-15},{30})<br>
><br>
> Rick<br>
><br>
> On Sun, Apr 19, 2020 at 1:54 PM zilore mumba via ncl-talk <<br>
> <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>> wrote:<br>
><br>
> I am trying to get time series at given locations from several WRF output<br>
> files, e.g. 48H forecast temperature at 850hPa. Each file contains a day's<br>
> model run.<br>
> I have used ncl for years by adapting existing scripts, without really<br>
> understanding. One issue I have difficulty with is getting variables from a<br>
> file, especially WRF output.<br>
> I thought the script to get a time series at 15S, 30E would be as simple<br>
> as the one below. I am able to print the varsummary, but printing x tells<br>
> me that "a is not a valid file reference.<br>
> I would appreciate ant assistance.<br>
> begin<br>
> files = systemfunc("ls -1 data/wrfout_d01_2020*") + ".nc"<br>
>   a = addfiles(files,"r")<br>
><br>
>   ListSetType (a, "cat")          ; concatenate (=default)<br>
>   tc = wrf_user_getvar(a,"tc",0)<br>
>   printVarSummary (tc)<br>
><br>
>   x = a->T2(0,{850},{-15},{30})<br>
>   asciiwrite("outfile.txt",x)<br>
> end<br>
><br>
> Zilore Mumba<br>
> _______________________________________________<br>
> ncl-talk mailing list<br>
> <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
> List instructions, subscriber options, unsubscribe:<br>
> <a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
><br>
> _______________________________________________<br>
> ncl-talk mailing list<br>
> <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
> List instructions, subscriber options, unsubscribe:<br>
> <a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
><br>
><br>
><br>
> --<br>
> Barry H. Lynn, Ph.D<br>
> Senior Associate Scientist, Lecturer,<br>
> The Institute of the Earth Science,<br>
> The Hebrew University of Jerusalem,<br>
> Givat Ram, Jerusalem 91904, Israel<br>
> Tel: 972 547 231 170<br>
> Fax: (972)-25662581<br>
><br>
> C.E.O, Weather It Is, LTD<br>
> Weather and Climate Focus<br>
> <a href="http://weather-it-is.com" rel="noreferrer" target="_blank">http://weather-it-is.com</a><br>
> Jerusalem, Israel<br>
> Local: 02 930 9525<br>
> Cell: 054 7 231 170<br>
> Int-IS: x972 2 930 9525<br>
><br>
><br>
<br>
-- <br>
Barry H. Lynn, Ph.D<br>
Senior Associate Scientist, Lecturer,<br>
The Institute of the Earth Science,<br>
The Hebrew University of Jerusalem,<br>
Givat Ram, Jerusalem 91904, Israel<br>
Tel: 972 547 231 170<br>
Fax: (972)-25662581<br>
<br>
C.E.O, Weather It Is, LTD<br>
Weather and Climate Focus<br>
<a href="http://weather-it-is.com" rel="noreferrer" target="_blank">http://weather-it-is.com</a><br>
Jerusalem, Israel<br>
Local: 02 930 9525<br>
Cell: 054 7 231 170<br>
Int-IS: x972 2 930 9525<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20200420/9841d672/attachment-0001.html" rel="noreferrer" target="_blank">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20200420/9841d672/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 16<br>
Date: Mon, 20 Apr 2020 19:04:41 +0800<br>
From: dickson mbigi <<a href="mailto:dickson.mbigi@gmail.com" target="_blank">dickson.mbigi@gmail.com</a>><br>
To: <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
Subject: Re: [ncl-talk] warning:ContourPlotInialize: no valid values<br>
        in scalar field; ContourPlot not possible: [errorno=1101]<br>
Message-ID:<br>
        <CALUZNzPT9traAfZgh8fvQD6i=AS=<a href="mailto:HwxzW_9WUqkgCbLw_qjuag@mail.gmail.com" target="_blank">HwxzW_9WUqkgCbLw_qjuag@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
Kindly help/directives,<br>
I am plotting map of  pre anomalies 1998-2018 with reference to the base<br>
period 1981-2010 at decadal time scale. I am getting the headed above<br>
warning along with no contours in the resulting plot. I think the problem<br>
occur at line 56 which I have used "wgt_runave_n_Wrap" function for<br>
filtering.<br>
1) I have checked my values to plot, they are all not  missing values.<br>
2) the coordinate variables of y (filtered values) are  also okay.<br>
<br>
I have attached scrip and resulting plot for your reference. I am using CRU<br>
pre data that are too big to attach along with this mail.<br>
<br>
Regards,<br>
Dickson<br>
Institute of Atmospheric Physics, CAS.<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20200420/28ae36e0/attachment-0001.html" rel="noreferrer" target="_blank">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20200420/28ae36e0/attachment-0001.html</a>><br>
-------------- next part --------------<br>
A non-text attachment was scrubbed...<br>
Name: preanomaly_decadal.ncl<br>
Type: application/octet-stream<br>
Size: 5259 bytes<br>
Desc: not available<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20200420/28ae36e0/attachment-0001.obj" rel="noreferrer" target="_blank">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20200420/28ae36e0/attachment-0001.obj</a>><br>
-------------- next part --------------<br>
A non-text attachment was scrubbed...<br>
Name: preanom_mam.png<br>
Type: image/png<br>
Size: 90674 bytes<br>
Desc: not available<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20200420/28ae36e0/attachment-0001.png" rel="noreferrer" target="_blank">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20200420/28ae36e0/attachment-0001.png</a>><br>
<br>
------------------------------<br>
<br>
Message: 17<br>
Date: Mon, 20 Apr 2020 07:24:41 -0400<br>
From: Alessandra Giannini <<a href="mailto:alesall@iri.columbia.edu" target="_blank">alesall@iri.columbia.edu</a>><br>
To: dickson mbigi <<a href="mailto:dickson.mbigi@gmail.com" target="_blank">dickson.mbigi@gmail.com</a>><br>
Cc: <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
Subject: Re: [ncl-talk] warning:ContourPlotInialize: no valid values<br>
        in scalar field; ContourPlot not possible: [errorno=1101]<br>
Message-ID: <<a href="mailto:B2CC8A15-4C41-4BD1-A3ED-545669B5DD3D@iri.columbia.edu" target="_blank">B2CC8A15-4C41-4BD1-A3ED-545669B5DD3D@iri.columbia.edu</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
<br>
hi there, Dickson, <br>
<br>
1) could be that the way you?ve set your contour levels does not match the values of the field you want to plot.<br>
In that regard, what is the output of the command line:<br>
<br>
printMinMax(y, 0)<br>
<br>
[line 62 in your script]<br>
<br>
2) but I also note that you?ve set:<br>
<br>
res2@mpLandFillColor       = ?white"<br>
<br>
[twice, at lines 123 and 131]<br>
<br>
What happens if you comment these two lines out?<br>
<br>
warm regards, alessandra<br>
<br>
<br>
<br>
<br>
<br>
<br>
?<br>
Alessandra Giannini<br>
IRI for Climate and Society - The Earth Institute at Columbia University<br>
P.O. Box 1000, Palisades NY 10964-8000<br>
phone/fax: +1 845 680-4473/4864 - email: <a href="mailto:alesall@iri.columbia.edu" target="_blank">alesall@iri.columbia.edu</a><br>
<br>
<br>
<br>
<br>
> On Apr 20, 2020, at 7:04 AM, dickson mbigi via ncl-talk <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>> wrote:<br>
> <br>
> Kindly help/directives,<br>
> I am plotting map of  pre anomalies 1998-2018 with reference to the base period 1981-2010 at decadal time scale. I am getting the headed above warning along with no contours in the resulting plot. I think the problem occur at line 56 which I have used "wgt_runave_n_Wrap" function for filtering. <br>
> 1) I have checked my values to plot, they are all not  missing values.<br>
> 2) the coordinate variables of y (filtered values) are  also okay.<br>
> <br>
> I have attached scrip and resulting plot for your reference. I am using CRU pre data that are too big to attach along with this mail.<br>
> <br>
> Regards,<br>
> Dickson<br>
> Institute of Atmospheric Physics, CAS.<br>
> <preanomaly_decadal.ncl><preanom_mam.png>_______________________________________________<br>
> ncl-talk mailing list<br>
> <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
> List instructions, subscriber options, unsubscribe:<br>
> <a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20200420/1b9c9d60/attachment-0001.html" rel="noreferrer" target="_blank">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20200420/1b9c9d60/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 18<br>
Date: Mon, 20 Apr 2020 11:56:56 +0000 (UTC)<br>
From: zilore mumba <<a href="mailto:zmumba@yahoo.com" target="_blank">zmumba@yahoo.com</a>><br>
To: Barry Lynn <<a href="mailto:barry.h.lynn@gmail.com" target="_blank">barry.h.lynn@gmail.com</a>><br>
Cc: Ncl-talk <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>><br>
Subject: Re: [ncl-talk] Getting time series from several files in ncl<br>
Message-ID: <<a href="mailto:469006943.2254548.1587383816611@mail.yahoo.com" target="_blank">469006943.2254548.1587383816611@mail.yahoo.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
 Thanks I have looked at examples of wrf_user_ll_to_ij, It does not help. Let me clarify two things1. I used tc just to test that I can at least access tc, but that syntax does not give me tc at a specific location.2. When I talk of a specific location, I do not need interpolation. I know WRF will give me the value at the nearest grid location. That is enough for me. My core problem is: how do I put it to get one value, instead of the whole grid?<br>
<br>
<br>
| <br>
| <br>
|  | <br>
wrf_user_ll_to_ij<br>
<br>
NCL WRFUserARW.ncl functions (lat/lon functions)<br>
 |<br>
<br>
 |<br>
<br>
 |<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
    On Monday, April 20, 2020, 12:02:28 PM GMT+2, Barry Lynn <<a href="mailto:barry.h.lynn@gmail.com" target="_blank">barry.h.lynn@gmail.com</a>> wrote:  <br>
<br>
 If you know when these values occur within the time series, then you can index them quite easily.<br>
Why don't you use wrf_user_getvar and then the wrf interpolation program.<br>
If you want them at a specific point, I think you know how to do this.? There is also a get location program<br>
<a href="https://www.ncl.ucar.edu/Document/Functions/WRF_arw/wrf_user_ll_to_ij.shtml" rel="noreferrer" target="_blank">https://www.ncl.ucar.edu/Document/Functions/WRF_arw/wrf_user_ll_to_ij.shtml</a><br>
<br>
Barry<br>
On Mon, Apr 20, 2020 at 11:55 AM zilore mumba <<a href="mailto:zmumba@yahoo.com" target="_blank">zmumba@yahoo.com</a>> wrote:<br>
<br>
 Barry,My intention is to have RH and U, V time series. That's what I had in mind when I put 48. T2 only has time, lon ,lat. I have several daily output files, and I want an n-day series of 48 hour forecasts of T2<br>
<br>
    On Monday, April 20, 2020, 8:00:41 AM GMT+2, Barry Lynn <<a href="mailto:barry.h.lynn@gmail.com" target="_blank">barry.h.lynn@gmail.com</a>> wrote:  <br>
<br>
 Zilore:<br>
T2 is a three dimensional variable.? Also, T2 is the two meter temperature.?<br>
You want tc or tk (see:?<a href="https://www.ncl.ucar.edu/Document/Functions/WRF_arw/wrf_user_getvar.shtml" rel="noreferrer" target="_blank">https://www.ncl.ucar.edu/Document/Functions/WRF_arw/wrf_user_getvar.shtml</a>).<br>
Could you please explain how you are interpolating to 850 mb?? There is a WRF program to do this, and I am not sure how the indexing you are using works. How are you indexing the 48 time slot?? WRF has a program to list times, and you can match the time you want to the listed times to find the index.<br>
Barry<br>
On Mon, Apr 20, 2020 at 8:17 AM zilore mumba via ncl-talk <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>> wrote:<br>
<br>
 Sorry Rick I put in a surface variable T2 which only 3 coordinates, so I put "x = a[:]->T2(0,{-15},{30})". Now I am getting error of coordinates west-east and south-north not having? associated coordinate variables. I have seen this error before. I will find out how to get out of it <br>
<br>
    On Sunday, April 19, 2020, 11:53:33 PM GMT+2, Rick Brownrigg <<a href="mailto:brownrig@ucar.edu" target="_blank">brownrig@ucar.edu</a>> wrote:  <br>
<br>
 Perhaps do a printVarSummary(a[0]->T2)? ; it sounds like its not the shape you think it might/should be; i.e. 3D vs 4D<br>
<br>
<br>
On Sun, Apr 19, 2020 at 3:29 PM zilore mumba <<a href="mailto:zmumba@yahoo.com" target="_blank">zmumba@yahoo.com</a>> wrote:<br>
<br>
<br>
<br>
    On Sunday, April 19, 2020, 11:27:47 PM GMT+2, zilore mumba <<a href="mailto:zmumba@yahoo.com" target="_blank">zmumba@yahoo.com</a>> wrote:  <br>
<br>
  Thanks Rick, that is close. The error am getting now is: fatal:Number of subscripts on rhs do not match number of dimensions of aggregated cat type variable, (4) Subscripts used, (3) Subscripts expected<br>
<br>
    On Sunday, April 19, 2020, 10:17:12 PM GMT+2, Rick Brownrigg <<a href="mailto:brownrig@ucar.edu" target="_blank">brownrig@ucar.edu</a>> wrote:  <br>
<br>
 HI Zilore,<br>
I think the issue is that addfiles() return a *list* of filevars, where addfile returns a single filevar.? Lists are indexed with [ ] notation; perhaps something like this will work:<br>
? x = a[:]->T2(0,{850},{-15},{30})<br>
Rick<br>
<br>
On Sun, Apr 19, 2020 at 1:54 PM zilore mumba via ncl-talk <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>> wrote:<br>
<br>
I am trying to get time series at given locations from several WRF output files, e.g. 48H forecast temperature at 850hPa. Each file contains a day's model run. I have used ncl for years by adapting existing scripts, without really understanding. One issue I have difficulty with is getting variables from a file, especially WRF output.I thought the script to get a time series at 15S, 30E would be as simple as the one below. I am able to print the varsummary, but printing x tells me that "a is not a valid file reference.I would appreciate ant assistance.begin<br>
files = systemfunc("ls -1 data/wrfout_d01_2020*") + ".nc"<br>
? a = addfiles(files,"r")<br>
<br>
? ListSetType (a, "cat")????????? ; concatenate (=default)<br>
? tc = wrf_user_getvar(a,"tc",0)<br>
? printVarSummary (tc)<br>
<br>
? x = a->T2(0,{850},{-15},{30})<br>
? asciiwrite("outfile.txt",x)<br>
end<br>
<br>
Zilore Mumba<br>
_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br>
  _______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br>
<br>
-- <br>
Barry H. Lynn, Ph.DSenior Associate Scientist, Lecturer,The Institute of the Earth Science,?<br>
The Hebrew University of Jerusalem,?<br>
Givat Ram, Jerusalem 91904, Israel?<br>
Tel: 972 547 231 170<br>
Fax: (972)-25662581<br>
C.E.O, Weather It Is, LTD<br>
Weather and Climate Focus<br>
<a href="http://weather-it-is.com" rel="noreferrer" target="_blank">http://weather-it-is.com</a><br>
Jerusalem, Israel<br>
Local: 02 930 9525<br>
Cell: 054 7 231 170<br>
Int-IS: x972 2 930 9525<br>
<br>
<br>
<br>
<br>
-- <br>
Barry H. Lynn, Ph.DSenior Associate Scientist, Lecturer,The Institute of the Earth Science,?<br>
The Hebrew University of Jerusalem,?<br>
Givat Ram, Jerusalem 91904, Israel?<br>
Tel: 972 547 231 170<br>
Fax: (972)-25662581<br>
C.E.O, Weather It Is, LTD<br>
Weather and Climate Focus<br>
<a href="http://weather-it-is.com" rel="noreferrer" target="_blank">http://weather-it-is.com</a><br>
Jerusalem, Israel<br>
Local: 02 930 9525<br>
Cell: 054 7 231 170<br>
Int-IS: x972 2 930 9525<br>
<br>
<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20200420/c9492d31/attachment-0001.html" rel="noreferrer" target="_blank">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20200420/c9492d31/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 19<br>
Date: Mon, 20 Apr 2020 05:33:48 -0700<br>
From: Rashed Mahmood <<a href="mailto:rashidcomsis@gmail.com" target="_blank">rashidcomsis@gmail.com</a>><br>
To: zilore mumba <<a href="mailto:zmumba@yahoo.com" target="_blank">zmumba@yahoo.com</a>><br>
Cc: Ncl-talk <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>><br>
Subject: Re: [ncl-talk] Getting time series from several files in ncl<br>
Message-ID:<br>
        <<a href="mailto:CAOmfQPp7FmB_XGk0etguc9i70dhawS4knd88y2gpuaePd9uJWg@mail.gmail.com" target="_blank">CAOmfQPp7FmB_XGk0etguc9i70dhawS4knd88y2gpuaePd9uJWg@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
Since WRF grid is a 2D grid (i.e. curvilinear grid), you would need to use<br>
special sub-setting function/s. See example 2 on the following page:<br>
<a href="https://www.ncl.ucar.edu/Applications/latlon_subset.shtml" rel="noreferrer" target="_blank">https://www.ncl.ucar.edu/Applications/latlon_subset.shtml</a><br>
<br>
On Mon, Apr 20, 2020 at 4:57 AM zilore mumba via ncl-talk <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>><br>
wrote:<br>
<br>
> Thanks I have looked at examples of wrf_user_ll_to_ij,<br>
> <<a href="https://www.ncl.ucar.edu/Document/Functions/WRF_arw/wrf_user_ll_to_ij.shtml" rel="noreferrer" target="_blank">https://www.ncl.ucar.edu/Document/Functions/WRF_arw/wrf_user_ll_to_ij.shtml</a>><br>
> It does not help. Let me clarify two things<br>
> 1. I used tc just to test that I can at least access tc, but that syntax<br>
> does not give me tc at a specific location.<br>
> 2. When I talk of a specific location, I do not need interpolation. I know<br>
> WRF will give me the value at the nearest grid location. That is enough for<br>
> me. My core problem is: how do I put it to get one value, instead of the<br>
> whole grid?<br>
><br>
> wrf_user_ll_to_ij<br>
><br>
> NCL WRFUserARW.ncl functions (lat/lon functions)<br>
><br>
> <<a href="https://www.ncl.ucar.edu/Document/Functions/WRF_arw/wrf_user_ll_to_ij.shtml" rel="noreferrer" target="_blank">https://www.ncl.ucar.edu/Document/Functions/WRF_arw/wrf_user_ll_to_ij.shtml</a>><br>
><br>
><br>
><br>
><br>
><br>
><br>
> On Monday, April 20, 2020, 12:02:28 PM GMT+2, Barry Lynn <<br>
> <a href="mailto:barry.h.lynn@gmail.com" target="_blank">barry.h.lynn@gmail.com</a>> wrote:<br>
><br>
><br>
> If you know when these values occur within the time series, then you can<br>
> index them quite easily.<br>
><br>
> Why don't you use wrf_user_getvar and then the wrf interpolation program.<br>
><br>
> If you want them at a specific point, I think you know how to do this.<br>
> There is also a get location program<br>
><br>
> <a href="https://www.ncl.ucar.edu/Document/Functions/WRF_arw/wrf_user_ll_to_ij.shtml" rel="noreferrer" target="_blank">https://www.ncl.ucar.edu/Document/Functions/WRF_arw/wrf_user_ll_to_ij.shtml</a><br>
><br>
> Barry<br>
><br>
> On Mon, Apr 20, 2020 at 11:55 AM zilore mumba <<a href="mailto:zmumba@yahoo.com" target="_blank">zmumba@yahoo.com</a>> wrote:<br>
><br>
> Barry,<br>
> My intention is to have RH and U, V time series. That's what I had in mind<br>
> when I put 48. T2 only has time, lon ,lat. I have several daily output<br>
> files, and I want an n-day series of 48 hour forecasts of T2<br>
><br>
> On Monday, April 20, 2020, 8:00:41 AM GMT+2, Barry Lynn <<br>
> <a href="mailto:barry.h.lynn@gmail.com" target="_blank">barry.h.lynn@gmail.com</a>> wrote:<br>
><br>
><br>
> Zilore:<br>
><br>
> T2 is a three dimensional variable.  Also, T2 is the two meter<br>
> temperature.<br>
><br>
> You want tc or tk (see:<br>
> <a href="https://www.ncl.ucar.edu/Document/Functions/WRF_arw/wrf_user_getvar.shtml" rel="noreferrer" target="_blank">https://www.ncl.ucar.edu/Document/Functions/WRF_arw/wrf_user_getvar.shtml</a><br>
> ).<br>
><br>
> Could you please explain how you are interpolating to 850 mb?  There is a<br>
> WRF program to do this, and I am not sure how the indexing you are using<br>
> works. How are you indexing the 48 time slot?  WRF has a program to list<br>
> times, and you can match the time you want to the listed times to find the<br>
> index.<br>
><br>
> Barry<br>
><br>
> On Mon, Apr 20, 2020 at 8:17 AM zilore mumba via ncl-talk <<br>
> <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>> wrote:<br>
><br>
> Sorry Rick I put in a surface variable T2 which only 3 coordinates, so I<br>
> put "x = a[:]->T2(0,{-15},{30})". Now I am getting error of coordinates<br>
> west-east and south-north not having  associated coordinate variables. I<br>
> have seen this error before. I will find out how to get out of it<br>
><br>
> On Sunday, April 19, 2020, 11:53:33 PM GMT+2, Rick Brownrigg <<br>
> <a href="mailto:brownrig@ucar.edu" target="_blank">brownrig@ucar.edu</a>> wrote:<br>
><br>
><br>
> Perhaps do a printVarSummary(a[0]->T2)  ; it sounds like its not the shape<br>
> you think it might/should be; i.e. 3D vs 4D<br>
><br>
><br>
><br>
> On Sun, Apr 19, 2020 at 3:29 PM zilore mumba <<a href="mailto:zmumba@yahoo.com" target="_blank">zmumba@yahoo.com</a>> wrote:<br>
><br>
><br>
><br>
> On Sunday, April 19, 2020, 11:27:47 PM GMT+2, zilore mumba <<br>
> <a href="mailto:zmumba@yahoo.com" target="_blank">zmumba@yahoo.com</a>> wrote:<br>
><br>
><br>
> Thanks Rick, that is close. The error am getting now is:<br>
> fatal:Number of subscripts on rhs do not match number of dimensions of<br>
> aggregated cat type variable, (4) Subscripts used, (3) Subscripts expected<br>
><br>
> On Sunday, April 19, 2020, 10:17:12 PM GMT+2, Rick Brownrigg <<br>
> <a href="mailto:brownrig@ucar.edu" target="_blank">brownrig@ucar.edu</a>> wrote:<br>
><br>
><br>
> HI Zilore,<br>
><br>
> I think the issue is that addfiles() return a *list* of filevars, where<br>
> addfile returns a single filevar.  Lists are indexed with [ ] notation;<br>
> perhaps something like this will work:<br>
><br>
>   x = a[:]->T2(0,{850},{-15},{30})<br>
><br>
> Rick<br>
><br>
> On Sun, Apr 19, 2020 at 1:54 PM zilore mumba via ncl-talk <<br>
> <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>> wrote:<br>
><br>
> I am trying to get time series at given locations from several WRF output<br>
> files, e.g. 48H forecast temperature at 850hPa. Each file contains a day's<br>
> model run.<br>
> I have used ncl for years by adapting existing scripts, without really<br>
> understanding. One issue I have difficulty with is getting variables from a<br>
> file, especially WRF output.<br>
> I thought the script to get a time series at 15S, 30E would be as simple<br>
> as the one below. I am able to print the varsummary, but printing x tells<br>
> me that "a is not a valid file reference.<br>
> I would appreciate ant assistance.<br>
> begin<br>
> files = systemfunc("ls -1 data/wrfout_d01_2020*") + ".nc"<br>
>   a = addfiles(files,"r")<br>
><br>
>   ListSetType (a, "cat")          ; concatenate (=default)<br>
>   tc = wrf_user_getvar(a,"tc",0)<br>
>   printVarSummary (tc)<br>
><br>
>   x = a->T2(0,{850},{-15},{30})<br>
>   asciiwrite("outfile.txt",x)<br>
> end<br>
><br>
> Zilore Mumba<br>
> _______________________________________________<br>
> ncl-talk mailing list<br>
> <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
> List instructions, subscriber options, unsubscribe:<br>
> <a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
><br>
> _______________________________________________<br>
> ncl-talk mailing list<br>
> <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
> List instructions, subscriber options, unsubscribe:<br>
> <a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
><br>
><br>
><br>
> --<br>
> Barry H. Lynn, Ph.D<br>
> Senior Associate Scientist, Lecturer,<br>
> The Institute of the Earth Science,<br>
> The Hebrew University of Jerusalem,<br>
> Givat Ram, Jerusalem 91904, Israel<br>
> Tel: 972 547 231 170<br>
> Fax: (972)-25662581<br>
><br>
> C.E.O, Weather It Is, LTD<br>
> Weather and Climate Focus<br>
> <a href="http://weather-it-is.com" rel="noreferrer" target="_blank">http://weather-it-is.com</a><br>
> Jerusalem, Israel<br>
> Local: 02 930 9525<br>
> Cell: 054 7 231 170<br>
> Int-IS: x972 2 930 9525<br>
><br>
><br>
><br>
> --<br>
> Barry H. Lynn, Ph.D<br>
> Senior Associate Scientist, Lecturer,<br>
> The Institute of the Earth Science,<br>
> The Hebrew University of Jerusalem,<br>
> Givat Ram, Jerusalem 91904, Israel<br>
> Tel: 972 547 231 170<br>
> Fax: (972)-25662581<br>
><br>
> C.E.O, Weather It Is, LTD<br>
> Weather and Climate Focus<br>
> <a href="http://weather-it-is.com" rel="noreferrer" target="_blank">http://weather-it-is.com</a><br>
> Jerusalem, Israel<br>
> Local: 02 930 9525<br>
> Cell: 054 7 231 170<br>
> Int-IS: x972 2 930 9525<br>
><br>
> _______________________________________________<br>
> ncl-talk mailing list<br>
> <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
> List instructions, subscriber options, unsubscribe:<br>
> <a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20200420/b976f71d/attachment.html" rel="noreferrer" target="_blank">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20200420/b976f71d/attachment.html</a>><br>
<br>
------------------------------<br>
<br>
Subject: Digest Footer<br>
<br>
_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br>
------------------------------<br>
<br>
End of ncl-talk Digest, Vol 197, Issue 20<br>
*****************************************<br>
</blockquote></div></div>