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