[ncl-talk] [NCL Digest] python transition (ali mughal)
ali mughal
mughalali655 at gmail.com
Sat Dec 5 08:39:11 MST 2020
Dear Marco
Thanks alot. But can you tell me if I can install PyNGL via pip install. I
am trying to install on Linux and user permissions from our administrator
donot allow me to install conda.
Also can I find some video tutorials and are there any trainings online
that I can attend
Thank you once again
On Saturday, December 5, 2020, Marco Miani via ncl-talk <
ncl-talk at mailman.ucar.edu> wrote:
> Hi Ali,
>
>
>
> Python has something called PyNGL which, as the name suggests, is a python
> based NCL platform.
>
> The structure is very similar to pure NCL, and so is the slicing of data:
>
>
>
> At the bottom of this message, I enclose a little snippet.
>
> More infos can be found here, on the official webpage:
> https://www.pyngl.ucar.edu/
>
>
>
> Did this help?
> Bye
>
>
>
> Marco
>
>
>
>
>
>
>
> #
>
> # File:
>
> # newcolor1.py
>
> #
>
> # Synopsis:
>
> # Illustrates new color capabilities in PyNGL 1.5.0.
>
> #
>
> # Categories:
>
> # Contouring
>
> #
>
> # Author:
>
> # Mary Haley
>
> #
>
> # Date of initial publication:
>
> # November 2012
>
> #
>
> # Description:
>
> # This example shows how to use the new cnFillPalette resource.
>
> #
>
> # Effects illustrated:
>
> # o Using the new "cnFillPalette" resource.
>
> # o Using a named color without having to add it to color map
>
> #
>
> # Output:
>
> # This example produces one visualization
>
> #
>
> # Notes:
>
> #
>
> from __future__ import print_function
>
> import os, numpy
>
> import Ngl, Nio
>
>
>
> #
>
> # Create some dummy data for the contour plot.
>
> #
>
> dirc = Ngl.pynglpath("data")
>
> f = Nio.open_file(os.path.join(dirc,"cdf","uv300.nc"))
>
> u = f.variables["U"][1,:,:] # ß---Ali, this is what you are looking for!
>
> lat = f.variables["lat"][:]
>
> lon = f.variables["lon"][:]
>
> wks_type = "png"
>
> wks = Ngl.open_wks(wks_type,"newcolor1")
>
>
>
> cnres = Ngl.Resources()
>
>
>
> # Contour resources
>
> cnres.cnFillOn = True
>
> cnres.cnFillPalette = "BlueYellowRed" # New in PyNGL 1.5.0
>
> cnres.cnLinesOn = False
>
> cnres.cnLineLabelsOn = False
>
>
>
> # Labelbar resource
>
> cnres.lbOrientation = "horizontal"
>
>
>
> # Scalar field resources
>
> cnres.sfXArray = lon
>
> cnres.sfYArray = lat
>
>
>
> # Map resources
>
> cnres.mpFillOn = True
>
> cnres.mpFillDrawOrder = "PostDraw"
>
> cnres.mpLandFillColor = "Gray"
>
> cnres.mpOceanFillColor = "Transparent"
>
> cnres.mpInlandWaterFillColor = "Transparent"
>
>
>
> contour = Ngl.contour_map(wks,u,cnres)
>
>
>
> Ngl.end()
>
>
>
>
>
> *Marco Miani*
> Computational Support Specialist
> Climate and Atmosphere Research Center (CARE-C)
> The Cyprus Institute
> 20 Konstantinou Kavafi Street, 2121, Nicosia, Cyprus
> *Tel:* +357 22 397 561 <+357%2022%20397%20561> | *Email:*
> m.miani at cyi.ac.cy
> *Web:* cyi.ac.cy | emme-care.cyi.ac.cy
>
> Follow us on <https://www.facebook.com/CyprusInstitute/>
> <https://twitter.com/CyprusInstitute>
> <https://www.linkedin.com/company/the-cyprus-institute>
>
> P *Please consider the environment before printing this email*
>
>
>
>
>
> *From: *ncl-talk-request at mailman.ucar.edu
> *Sent: *04 December 2020 21:00
> *To: *ncl-talk at mailman.ucar.edu
> *Subject: *ncl-talk Digest, Vol 205, Issue 5
>
>
>
> Send ncl-talk mailing list submissions to
> ncl-talk at mailman.ucar.edu
>
> To subscribe or unsubscribe via the World Wide Web, visit
> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
> or, via email, send a message with subject or body 'help' to
> ncl-talk-request at mailman.ucar.edu
>
> You can reach the person managing the list at
> ncl-talk-owner at mailman.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. Re: Setting X-axis for proper data plotting (Adam Phillips)
> 2. Effective Degree of Freedom (dickson mbigi)
> 3. Re: Plot csv file (Anahita Amiri Farahani)
> 4. python transition (ali mughal)
> 5. wrfout file reading (Anil Kumar - NOAA Affiliate)
> 6. Re: wrfout file reading (Ehsan Taghizadeh)
> 7. Re: wrfout file reading (Anil Kumar - NOAA Affiliate)
> 8. Re: wrfout file reading (Anil Kumar - NOAA Affiliate)
>
>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Thu, 3 Dec 2020 18:12:48 -0700
> From: Adam Phillips <asphilli at ucar.edu>
> To: "Herb, Jason" <jherb at albany.edu>
> Cc: "ncl-talk at ucar.edu" <ncl-talk at ucar.edu>
> Subject: Re: [ncl-talk] Setting X-axis for proper data plotting
> Message-ID:
> <CACvWuqB=NnFQjPePoRUrrB+Okd74DLs+Bj3-ttjcdTzjv_ONcg@
> mail.gmail.com>
> Content-Type: text/plain; charset="utf-8"
>
> Hi Jason,
> Thank you for sending an easy-to-run script. I have attached a modified
> script along with the resulting graphic.
> Two things to note:
> 1 - I altered your asciiread statement to read in the data as a double,
> that way you can see the correct time. I then manually created my own
> yyyymmddhh array from the original time. Note that I assumed that the data
> is regular and that there are not any missing timesteps from
> 20061221-20191231.95833. Once I passed the yyyymmddhh array into
> yyyymmddhh_to_yyyyfrac I now have a regular time coordinate to pass into
> gsn_csm_xy. This alleviates the jumping you were getting in your plot.
> Please look at the print statement output and coding and make sure the
> assumptions I made are correct.
>
> 2 - To get better control of the tick mark labels, I recommend setting
> tmXBMode = "Explicit", and setting tmXBValues/tmXBLabels. See the
> tickmarks application
> page <http://www.ncl.ucar.edu/Applications/tickmarks.shtml> for other
> examples of using these resources.
>
> Hope that helps!
> Adam
>
> On Wed, Dec 2, 2020 at 2:39 PM Herb, Jason via ncl-talk <
> ncl-talk at mailman.ucar.edu> wrote:
>
> > Hello All,
> >
> > I have been working with a continuous data set from late 2006 to end of
> > 2019. I have been able to force the X-axis labeling as close as I can to
> > get year to whole numbers, but I keep getting ".0" at the end of the
> year.
> > Also, I have a plotting error in the actual data with the space
> containing
> > "flat lines" between the end or one year to the start of the next. I
> > attempted to setup another array for the year but get the following error
> > message :
> >
> > fatal:Illegal subscript. Subscripts must be integer when not using
> > coordinate indexing
> > fatal:["Execute.c":8640]:Execute: Error occurred at or near line 36 in
> > file Basiclineplot.ncl
> >
> > I have attached the plot, script and input file.
> > Does anyone have any suggestions to assist with this plotting issue?
> >
> > Thanks,
> > Jason
> > _______________________________________________
> > ncl-talk mailing list
> > ncl-talk at mailman.ucar.edu
> > List instructions, subscriber options, unsubscribe:
> > https://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
>
> --
> Adam Phillips
> Associate Scientist, Climate and Global Dynamics Laboratory, NCAR
> www.cgd.ucar.edu/staff/asphilli/ 303-497-1726
>
> <http://www.cgd.ucar.edu/staff/asphilli>
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: <https://mailman.ucar.edu/pipermail/ncl-talk/
> attachments/20201203/6db121df/attachment-0001.html>
> -------------- next part --------------
> A non-text attachment was scrubbed...
> Name: test.png
> Type: image/png
> Size: 75791 bytes
> Desc: not available
> URL: <https://mailman.ucar.edu/pipermail/ncl-talk/
> attachments/20201203/6db121df/attachment-0001.png>
> -------------- next part --------------
> A non-text attachment was scrubbed...
> Name: Basiclineplot.ncl
> Type: application/octet-stream
> Size: 1663 bytes
> Desc: not available
> URL: <https://mailman.ucar.edu/pipermail/ncl-talk/
> attachments/20201203/6db121df/attachment-0001.obj>
>
> ------------------------------
>
> Message: 2
> Date: Fri, 4 Dec 2020 11:23:49 +0800
> From: dickson mbigi <dickson.mbigi at gmail.com>
> To: ncl-talk at mailman.ucar.edu
> Subject: [ncl-talk] Effective Degree of Freedom
> Message-ID:
> <CALUZNzMQqxTPetw2OtW_nNt9px6Jh0DHeNPdDupykrvqSPxSDA
> @mail.gmail.com>
> Content-Type: text/plain; charset="utf-8"
>
> Okay got it.
> Thanks
>
> On Fri, Dec 4, 2020 at 2:21 AM <ncl-talk-request at mailman.ucar.edu> wrote:
>
> > Send ncl-talk mailing list submissions to
> > ncl-talk at mailman.ucar.edu
> >
> > To subscribe or unsubscribe via the World Wide Web, visit
> > https://mailman.ucar.edu/mailman/listinfo/ncl-talk
> > or, via email, send a message with subject or body 'help' to
> > ncl-talk-request at mailman.ucar.edu
> >
> > You can reach the person managing the list at
> > ncl-talk-owner at mailman.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. Effective Degree of Freedom (dickson mbigi)
> > 2. Re: NCARG_COLORMAPS in conda (Jian Zheng)
> > 3. Re: Effective Degree of Freedom (Dennis Shea)
> > 4. Plot csv file (Anahita Amiri Farahani)
> >
> >
> > ----------------------------------------------------------------------
> >
> > Message: 1
> > Date: Thu, 3 Dec 2020 08:25:59 +0800
> > From: dickson mbigi <dickson.mbigi at gmail.com>
> > To: NCL-talk <ncl-talk at ucar.edu>
> > Subject: [ncl-talk] Effective Degree of Freedom
> > Message-ID:
> > <CALUZNzOzszoPZz5xRtZeLVbZYaNkXpVy=
> > q+RRBY5sVm4uxUiqw at mail.gmail.com>
> > Content-Type: text/plain; charset="utf-8"
> >
> > Hi everyone,
> > I have sst(time) and pre(lon,lat,time) low pass filtered data sets. I can
> > calculate spatial correlation, but for the significance test using
> > Student's t- test I need to use effective degrees of freedom because my
> > data have been filtered. I went through NCL website and found
> ''equivalent
> > sample size'' function which I think is quite different from the
> effective
> > degrees of freedom. Is there a function or a way to use Effective Degree
> of
> > freedom to perform significance tests in spatial correlation analysis?
> >
> > Ndof =N(1-r1r2)/(1+r1r2), where N is sample size, r1&r2 are the lag-1
> > autocorrelations of the two time series.
> >
> > Thanks in advance.
> > Dickson.
> > -------------- next part --------------
> > An HTML attachment was scrubbed...
> > URL: <
> > https://mailman.ucar.edu/pipermail/ncl-talk/
> attachments/20201203/cfbd359e/attachment-0001.html
> > >
> >
> > ------------------------------
> >
> > Message: 2
> > Date: Thu, 3 Dec 2020 11:13:39 +0800
> > From: Jian Zheng <jzheng8606 at gmail.com>
> > To: Rick Brownrigg <brownrig at ucar.edu>
> > Cc: Ncl-talk <ncl-talk at ucar.edu>
> > Subject: Re: [ncl-talk] NCARG_COLORMAPS in conda
> > Message-ID:
> > <
> > CAFMEVXN2tHnLXEYVyhdMb4iR3OL0Swd+3FEZNp8PETiNL09SxA at mail.gmail.com>
> > Content-Type: text/plain; charset="utf-8"
> >
> > Hi Rick,
> >
> > I find that defining NCARG_COLORMAPS after the conda-initializing part in
> > the .bash_profile can solve the problem.
> >
> > Thanks a lot for all your help.
> > Jian
> >
> > On Sat, Nov 21, 2020 at 12:22 AM Rick Brownrigg <brownrig at ucar.edu>
> wrote:
> >
> > > Hi Jian,
> > >
> > > When you say it is not working, are you getting messages such as this
> > when
> > > you try to reference your custom colormap:
> > >
> > > fatal:CvtStringToCmap:Unable to convert string "foo" to ColorMap
> > > warning:Error retrieving resource wkColorMap from args - Ignoring Arg
> > >
> > > I don't know where the OLD_NCARG_COLORMAPS setting would be coming from
> > in
> > > your environment. I would note that the second path in it --
> > > /lib/ncarg/colormaps -- would be incorrect. I don't know that it would
> > > break things if the first path is correct however.
> > >
> > > I have a stock ncl installation from conda. I copied an arbitrary
> > colormap
> > > from $NCARG_ROOT/lib/colormaps into /tmp, renamed it "foo.rgb", and set
> > > NCARG_COLORMAPS=/tmp. Then when I run this little script, everything
> > works:
> > >
> > > ncl 1> wks = gsn_open_wks("x11","test")
> > > ncl 2> gsn_define_colormap(wks,"foo")
> > > ncl 3> gsn_draw_colormap(wks)
> > >
> > > if I unset NCARG_COLORMAPS, and re-run, then I get the fatal/warning
> > > errors mentioned above.
> > >
> > > So I would have to conclude something is not quite right in your
> > > environment. Perhaps before trying to set envvars via .bash_profile,
> > > perhaps try setting them on the command line. Are you certain your
> > colormap
> > > is in the directory you cite? Are you certain its format is valid for
> > NCL?
> > > Perhaps try what I did -- copying a known colormap into your custom
> > > directory, rename it to something that does not collide with NCL's
> > colormap
> > > names, and see if you can get the little scriptlet above to work.
> > >
> > > I hope that helps...
> > > Rick
> > >
> > >
> > > On Fri, Nov 20, 2020 at 8:54 AM Jian Zheng <jzheng8606 at gmail.com>
> wrote:
> > >
> > >> Hi Rick,
> > >>
> > >> When I ran "env | grep NCARG", I got
> > >>
> > >>
> > >>
> > OLD_NCARG_COLORMAPS=/Users/jzheng/Documents/nclcolormap:/
> lib/ncarg/colormaps
> > >>
> > >> NCARG_ROOT=/Users/jzheng/opt/miniconda3
> > >>
> > >>
> > >> Does it show you more information?
> > >>
> > >>
> > >> On Fri, Nov 20, 2020 at 12:49 PM Rick Brownrigg <brownrig at ucar.edu>
> > >> wrote:
> > >>
> > >>> Hi Jian,
> > >>>
> > >>> Offhand that all looks reasonable. Only you can verify the custom
> path.
> > >>> Did you "source" your .bash_profile after you modified it? I'm not
> > sure
> > >>> what else to suggest.
> > >>>
> > >>> Rick
> > >>>
> > >>> On Thu, Nov 19, 2020 at 9:34 PM Jian Zheng <jzheng8606 at gmail.com>
> > wrote:
> > >>>
> > >>>> Hi Rick,
> > >>>>
> > >>>> NCARG_COLORMAPS (export
> > >>>>
> > NCARG_COLORMAPS=/Users/jzheng/Documents/nclcolormap:$NCARG_
> ROOT/lib/ncarg/colormaps)
> > >>>> contains the directory where I put my user-defined colormap files.
> > When I
> > >>>> run ncl, it doesn't see these colormaps.
> > >>>>
> > >>>> Thanks,
> > >>>> Jian
> > >>>>
> > >>>> On Thu, Nov 19, 2020 at 9:33 PM Rick Brownrigg <brownrig at ucar.edu>
> > >>>> wrote:
> > >>>>
> > >>>>> Hi,
> > >>>>>
> > >>>>> What do you mean "it does not work"? What happens? Normally one
> > >>>>> should not need to set that environment variable unless there's
> > something
> > >>>>> unusual about the way the code is installed. It is common however
> > to set
> > >>>>> NCARG_ROOT, and that should point to the directory containing the
> NCL
> > >>>>> installation; i.e., the directory with the bin/ lib/ include/ and
> > man/
> > >>>>> directories belonging to NCL. I don't recall offhand where a stock
> > conda
> > >>>>> install places those directories.
> > >>>>>
> > >>>>> Rick
> > >>>>>
> > >>>>>
> > >>>>> On Thu, Nov 19, 2020 at 1:47 AM Jian Zheng via ncl-talk <
> > >>>>> ncl-talk at mailman.ucar.edu> wrote:
> > >>>>>
> > >>>>>> Hi,
> > >>>>>>
> > >>>>>> I installed ncl6.6.2 in conda root environment. When I set
> > NCARG_COLORMAPS
> > >>>>>> in my ~/.bash_profile file, it did not work. How should I set
> > NCARG_COLORMAPS
> > >>>>>> with conda correctly?
> > >>>>>>
> > >>>>>> Thanks,
> > >>>>>> Jian
> > >>>>>> _______________________________________________
> > >>>>>> ncl-talk mailing list
> > >>>>>> ncl-talk at mailman.ucar.edu
> > >>>>>> List instructions, subscriber options, unsubscribe:
> > >>>>>> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
> > >>>>>
> > >>>>>
> > -------------- next part --------------
> > An HTML attachment was scrubbed...
> > URL: <
> > https://mailman.ucar.edu/pipermail/ncl-talk/
> attachments/20201203/e73faef6/attachment-0001.html
> > >
> >
> > ------------------------------
> >
> > Message: 3
> > Date: Thu, 3 Dec 2020 09:44:35 -0700
> > From: Dennis Shea <shea at ucar.edu>
> > To: dickson mbigi <dickson.mbigi at gmail.com>
> > Cc: NCL-talk <ncl-talk at ucar.edu>
> > Subject: Re: [ncl-talk] Effective Degree of Freedom
> > Message-ID:
> > <CAOF1d_6CO1djz=
> > ZCQgjRGhqtuw1yEoxe4n62APfcgJHaXJ7sgw at mail.gmail.com>
> > Content-Type: text/plain; charset="utf-8"
> >
> > I believe you have asked this question before.
> >
> > ****ncl-talk is not a statistics forum. ****
> > * ***Please talk to a statistician****
> >
> > ====
> > (1) Possibly:
> > You could use *esacr*
> > <http://www.ncl.ucar.edu/Document/Functions/Built-in/esacr.shtml> or
> > *esacr_n*
> > <http://www.ncl.ucar.edu/Document/Functions/Built-in/esacr_n.shtml>
> > (esacr_n is available in NCL 6.5.0 and later) on the filtered series to
> get
> > the lag_1 autocorrelation coefficients.
> >
> > r_sst = esacr(sst,1)
> > r1_sst = r_sst(1) ; scalar; r1_sst(1)
> > r1_sst at long_name = "Lag-1 autocorrelation: sst
> > print(r1_sst)
> > ; pre(lat,lon,time) ... time is
> > rightmost dimension
> > r_pre = esacr(pre,1) ; r_pre(2,lat,lon,time)
> > r1_pre = r_pre(1,:,:) ; r1_pre(nlat,mlom)
> > r1_pre at long_name = "Lag-1 autocorrelation: pre"
> > printVarSummary(r1_pre)
> > printMinMax(r1_pre)
> >
> >
> > r1_sst *:=* conform(r1_pre, r1_sst, -1) ; := overwrite
> > printVarSummary(r1_sst) ; (nlat,mlon)
> >
> > r1_sst_pre = r1_sst*r1_pre ; array (grid point-by-grid
> > point) multiply
> >
> > Ndof = N*(1-r1_sst_pre)/(1+r1_sst_pre), where N is the original sample
> > size of time series
> >
> > printVarSummary(Ndof) ; (nlat,mlon)
> >
> > printMinMax(Ndof,0)
> >
> >
> > ====
> > (2)
> > Perhaps, a better approach:
> >
> > DOI: 10.1016/S0275-1062(04)90014-8
> > (2004) Determination of the degree of freedom of digital filtered time
> > series with an application to the correlation analysis between the length
> > of day and the Southern oscillation index
> > ====
> >
> > Again, it is best to talk with a statistician.
> >
> >
> > On Wed, Dec 2, 2020 at 5:26 PM dickson mbigi via ncl-talk <
> > ncl-talk at mailman.ucar.edu> wrote:
> >
> > > Hi everyone,
> > > I have sst(time) and pre(lon,lat,time) low pass filtered data sets. I
> can
> > > calculate spatial correlation, but for the significance test using
> > > Student's t- test I need to use effective degrees of freedom because
> my
> > > data have been filtered. I went through NCL website and found
> > ''equivalent
> > > sample size'' function which I think is quite different from the
> > effective
> > > degrees of freedom. Is there a function or a way to use Effective
> Degree
> > of
> > > freedom to perform significance tests in spatial correlation analysis?
> > >
> > > Ndof =N(1-r1r2)/(1+r1r2), where N is sample size, r1&r2 are the lag-1
> > > autocorrelations of the two time series.
> > >
> > > Thanks in advance.
> > > Dickson.
> > > _______________________________________________
> > > ncl-talk mailing list
> > > ncl-talk at mailman.ucar.edu
> > > List instructions, subscriber options, unsubscribe:
> > > https://mailman.ucar.edu/mailman/listinfo/ncl-talk
> > -------------- next part --------------
> > An HTML attachment was scrubbed...
> > URL: <
> > https://mailman.ucar.edu/pipermail/ncl-talk/
> attachments/20201203/7b26972c/attachment-0001.html
> > >
> >
> > ------------------------------
> >
> > Message: 4
> > Date: Thu, 3 Dec 2020 13:21:04 -0500
> > From: Anahita Amiri Farahani <aamir003 at ucr.edu>
> > To: "ncl-talk [ncl-talk]" <ncl-talk at ucar.edu>
> > Subject: [ncl-talk] Plot csv file
> > Message-ID:
> > <
> > CAABdVry-zwn-y4omPNi2VNYVJPQOUP_W5ZqEUb7ZhsBTphScvA at mail.gmail.com>
> > Content-Type: text/plain; charset="utf-8"
> >
> > Hello all,
> >
> > I tried to read and plot the csv data, but I faced some errors. I've
> > attached both data and the code.
> >
> > Thank you,
> > ana
> > -------------- next part --------------
> > An HTML attachment was scrubbed...
> > URL: <
> > https://mailman.ucar.edu/pipermail/ncl-talk/
> attachments/20201203/a266e1eb/attachment.html
> > >
> > -------------- next part --------------
> > A non-text attachment was scrubbed...
> > Name: Book4.csv
> > Type: text/csv
> > Size: 1320725 bytes
> > Desc: not available
> > URL: <
> > https://mailman.ucar.edu/pipermail/ncl-talk/
> attachments/20201203/a266e1eb/attachment.csv
> > >
> > -------------- next part --------------
> > A non-text attachment was scrubbed...
> > Name: test.ncl
> > Type: application/octet-stream
> > Size: 2564 bytes
> > Desc: not available
> > URL: <
> > https://mailman.ucar.edu/pipermail/ncl-talk/
> attachments/20201203/a266e1eb/attachment.obj
> > >
> >
> > ------------------------------
> >
> > Subject: Digest Footer
> >
> > _______________________________________________
> > ncl-talk mailing list
> > ncl-talk at mailman.ucar.edu
> > List instructions, subscriber options, unsubscribe:
> > https://mailman.ucar.edu/mailman/listinfo/ncl-talk
> >
> > ------------------------------
> >
> > End of ncl-talk Digest, Vol 205, Issue 3
> > ****************************************
> >
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: <https://mailman.ucar.edu/pipermail/ncl-talk/
> attachments/20201204/539c7862/attachment-0001.html>
>
> ------------------------------
>
> Message: 3
> Date: Fri, 4 Dec 2020 00:05:25 -0500
> From: Anahita Amiri Farahani <aamir003 at ucr.edu>
> To: Dennis Shea <shea at ucar.edu>
> Cc: "ncl-talk \[ncl-talk\]" <ncl-talk at ucar.edu>
> Subject: Re: [ncl-talk] Plot csv file
> Message-ID:
> <CAABdVrzKKXHmB--=s0q0XHdzDDLG7YJyBQkKZq4THzA=57
> 4VLg at mail.gmail.com>
> Content-Type: text/plain; charset="utf-8"
>
> Thank you so much!
>
> On Thu, Dec 3, 2020 at 7:12 PM Dennis Shea <shea at ucar.edu> wrote:
>
> > There are several *wrong* assignments. For example, the following is not
> > correct.
> >
> > total= new((/5,21488/),typeof(bc),bc at _FillValue)
> > [SNIP]
> > total!0 = "lat"
> > total!1 = "lon"
> > ======
> > The following dimension names would be more appropriate for the 'total'
> > variable. There are 5 (=ncol) variables.
> >
> > N = *dimsizes*(lat1)
> > total= *new*((/ncol,N/),*typeof*(bc),bc at _FillValue)
> > total!0 = "variable" ; 0->longitude, 1->latitude, 2->NO,
> > 3->NO2, 4->BC
> > total!1 = "value"
> > =======
> > The attached is a modification of the following example.
> >
> > See: *http://www.ncl.ucar.edu/Applications/station.shtml*
> > <http://www.ncl.ucar.edu/Applications/station.shtml>
> > Example 1
> >
> >
> >
> >
> > On Thu, Dec 3, 2020 at 11:22 AM Anahita Amiri Farahani via ncl-talk <
> > ncl-talk at mailman.ucar.edu> wrote:
> >
> >> Hello all,
> >>
> >> I tried to read and plot the csv data, but I faced some errors. I've
> >> attached both data and the code.
> >>
> >> Thank you,
> >> ana
> >> _______________________________________________
> >> ncl-talk mailing list
> >> ncl-talk at mailman.ucar.edu
> >> List instructions, subscriber options, unsubscribe:
> >> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
> >
> >
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: <https://mailman.ucar.edu/pipermail/ncl-talk/
> attachments/20201204/4f22cf88/attachment-0001.html>
>
> ------------------------------
>
> Message: 4
> Date: Fri, 4 Dec 2020 16:31:01 +0800
> From: ali mughal <mughalali655 at gmail.com>
> To: Ncl-talk <ncl-talk at ucar.edu>
> Subject: [ncl-talk] python transition
> Message-ID:
> <CAB5RZE2H_p3m0YOUDi6sqtPKF+ZnQ-sPAx-jt6FxgbwXP+6LRw at mail.
> gmail.com>
> Content-Type: text/plain; charset="utf-8"
>
> Dear NCL users
>
> Could you tell me the equivalent of the following NCL command in python ?
>
> TC2=TA_MEAN(:,0,:,:)
>
> where TA_MEAN is a 4d variable with the following dimensions
>
> Dimensions and sizes: [Time | 745] x [bottom_top | 50] x [south_north |
> 129] x [west_east | 210])
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: <https://mailman.ucar.edu/pipermail/ncl-talk/
> attachments/20201204/87497f4a/attachment-0001.html>
>
> ------------------------------
>
> Message: 5
> Date: Fri, 4 Dec 2020 09:27:36 -0500
> From: Anil Kumar - NOAA Affiliate <anil.kumar at noaa.gov>
> To: ncl-talk at mailman.ucar.edu
> Subject: [ncl-talk] wrfout file reading
> Message-ID:
> <CA+8gui_--oMeV9bJXmhjJoG-cRCMsVpNwAcu+aE+n=cx18rRhw@
> mail.gmail.com>
> Content-Type: text/plain; charset="utf-8"
>
> Hi,
> I am trying to read multiple wrfout files. I am reading okay..no problem
> with the code extracting variables.
> Problem comes, when each file has multiple times in it and then I don't
> know how to extract the time series of variables .
> Below is the code, please let me know how to modify code to read many times
> within each file.
> ************************
> fs1 = systemfunc("ls ../WRFV3/run/wrfout_d01*")
> nfs1= dimsizes(fs1)
> f1 = addfiles(fs1+".nc","r")
> ntimes1 = dimsizes(f1) ; number of times in the file(s)
>
> loc1 = wrf_user_ij_to_ll(f1, 593., 546., True)
> print("lon/lat locations are: " + loc)
> uu1 = new((/ntimes1/),float)
> vv1 = new((/ntimes1/),float)
> psfc1 = new((/ntimes1/),float)
> i1=0
> do while (i1.le.ntimes1-1)
> uu1(i1) = f1[i1]->U(:,0,546,593) ; (time, lat,lon) ; (2,64,128)
> vv1(i1) = f1[i1]->V(:,0,546,593) ; V10(:,29.301,-94.797)
> psfc1(i1) = f1[i1]->PSFC(:,546,593)
> i1 = i1 + 1
> end do
> model_wspd1 = sqrt(uu1^2+vv1^2) ;wind_speed(uu,vv)
> model_wdir1 = wind_direction(uu1,vv1,0)
> model_psfc1 = psfc1/100.
> ************************
> --
> Dr. Anil Kumar
> Environmental Modeling Center/NCEP/NWS/NOAA
> College Park, MD, USA
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: <https://mailman.ucar.edu/pipermail/ncl-talk/
> attachments/20201204/3deb1c85/attachment-0001.html>
>
> ------------------------------
>
> Message: 6
> Date: Fri, 4 Dec 2020 17:42:27 +0000 (UTC)
> From: Ehsan Taghizadeh <ehsantaghizadeh at yahoo.com>
> To: Anil Kumar - NOAA Affiliate <anil.kumar at noaa.gov>
> Cc: Ncl-talk <ncl-talk at ucar.edu>
> Subject: Re: [ncl-talk] wrfout file reading
> Message-ID: <2131133709.2884490.1607103747892 at mail.yahoo.com>
> Content-Type: text/plain; charset="utf-8"
>
> Dear Dr. AnilAre your wrfout files hourly?However, why you don't use
> wrf_user_list_times to extract the list of available times?
> fs1 = systemfunc("ls ../WRFV3/run/wrfout_d01*")
> nfs1= dimsizes(fs1)do i = 0, nfs1-1
> ?f1 = addfile(fs1(i)+".nc","r") times = wrf_user_list_times(f1) ; get
> times in the file ntimes = dimsizes(times) ; number of times in
> the file(s)
> end doSincerely
> Ehsan
> On Friday, December 4, 2020, 05:57:54 PM GMT+3:30, Anil Kumar - NOAA
> Affiliate via ncl-talk <ncl-talk at mailman.ucar.edu> wrote:
>
> Hi,I am trying to read multiple wrfout files. I am reading okay..no
> problem with the code extracting variables.?Problem comes, when each file
> has multiple times in it and then I don't know how to extract the time
> series of variables .?Below is the code, please let me know how to modify
> code to read many times within each file.************************fs1 =
> systemfunc("ls ../WRFV3/run/wrfout_d01*")
> ?nfs1= dimsizes(fs1)
> ?f1 = addfiles(fs1+".nc","r")
> ? ntimes1 =? dimsizes(f1) ? ? ? ? ; number of times in the file(s)
>
> loc1 = wrf_user_ij_to_ll(f1, 593., 546., True)
> ?print("lon/lat locations are: " + loc)
> uu1 ? ? ?= new((/ntimes1/),float)
> vv1 ? ? ?= new((/ntimes1/),float)
> psfc1 ? ? ?= new((/ntimes1/),float)
> i1=0
> do while (i1.le.ntimes1-1)
> uu1(i1) ? ?= f1[i1]->U(:,0,546,593) ? ; (time, lat,lon) ? ?; (2,64,128)
> vv1(i1) ? ?= f1[i1]->V(:,0,546,593) ? ? ? ?; V10(:,29.301,-94.797)
> psfc1(i1) ?= f1[i1]->PSFC(:,546,593)
> i1 = i1 + 1
> end do
> ? ?model_wspd1 = sqrt(uu1^2+vv1^2) ;wind_speed(uu,vv)
> ? ?model_wdir1 = wind_direction(uu1,vv1,0)
> ? ?model_psfc1 = psfc1/100.
> ************************--
> Dr. Anil KumarEnvironmental Modeling Center/NCEP/NWS/NOAA
> College Park, MD, USA??
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at mailman.ucar.edu
> List instructions, subscriber options, unsubscribe:
> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: <https://mailman.ucar.edu/pipermail/ncl-talk/
> attachments/20201204/fc1b9e52/attachment-0001.html>
>
> ------------------------------
>
> Message: 7
> Date: Fri, 4 Dec 2020 13:01:45 -0500
> From: Anil Kumar - NOAA Affiliate <anil.kumar at noaa.gov>
> To: Ehsan Taghizadeh <ehsantaghizadeh at yahoo.com>
> Cc: Ncl-talk <ncl-talk at ucar.edu>
> Subject: Re: [ncl-talk] wrfout file reading
> Message-ID:
> <CA+8gui_xL-OFWfeLrL8YsQ+x0=oZXYty6y=LHssjCvY21ijWhg at mail.
> gmail.com>
> Content-Type: text/plain; charset="utf-8"
>
> Yes, wrf outputs are hourly files. Each hourly file has 10 multiple times
> (6 mins interval).
>
> Thanks,
> Anil
>
>
> On Fri, Dec 4, 2020 at 12:42 PM Ehsan Taghizadeh <
> ehsantaghizadeh at yahoo.com>
> wrote:
>
> > Dear Dr. Anil
> > Are your wrfout files hourly?
> > However, why you don't use wrf_user_list_times to extract the list of
> > available times?
> >
> > fs1 = systemfunc("ls ../WRFV3/run/wrfout_d01*")
> > nfs1= dimsizes(fs1)
> > do i = 0, nfs1-1
> > f1 = addfile(fs1(i)+".nc","r")
> > times = *wrf_user_list_times*(f1) ; get times in the file
> >
> > ntimes = dimsizes <https://www.ncl.ucar.edu/
> Document/Functions/Built-in/dimsizes.shtml>(times) ; number of
> times in the file(s)
> > end do
> >
> > *Sincerely*
> > *Ehsan*
> >
> > On Friday, December 4, 2020, 05:57:54 PM GMT+3:30, Anil Kumar - NOAA
> > Affiliate via ncl-talk <ncl-talk at mailman.ucar.edu> wrote:
> >
> >
> > Hi,
> > I am trying to read multiple wrfout files. I am reading okay..no problem
> > with the code extracting variables.
> > Problem comes, when each file has multiple times in it and then I don't
> > know how to extract the time series of variables .
> > Below is the code, please let me know how to modify code to read many
> > times within each file.
> > ************************
> > fs1 = systemfunc("ls ../WRFV3/run/wrfout_d01*")
> > nfs1= dimsizes(fs1)
> > f1 = addfiles(fs1+".nc","r")
> > ntimes1 = dimsizes(f1) ; number of times in the file(s)
> >
> > loc1 = wrf_user_ij_to_ll(f1, 593., 546., True)
> > print("lon/lat locations are: " + loc)
> > uu1 = new((/ntimes1/),float)
> > vv1 = new((/ntimes1/),float)
> > psfc1 = new((/ntimes1/),float)
> > i1=0
> > do while (i1.le.ntimes1-1)
> > uu1(i1) = f1[i1]->U(:,0,546,593) ; (time, lat,lon) ; (2,64,128)
> > vv1(i1) = f1[i1]->V(:,0,546,593) ; V10(:,29.301,-94.797)
> > psfc1(i1) = f1[i1]->PSFC(:,546,593)
> > i1 = i1 + 1
> > end do
> > model_wspd1 = sqrt(uu1^2+vv1^2) ;wind_speed(uu,vv)
> > model_wdir1 = wind_direction(uu1,vv1,0)
> > model_psfc1 = psfc1/100.
> > ************************
> > --
> > Dr. Anil Kumar
> > Environmental Modeling Center/NCEP/NWS/NOAA
> > College Park, MD, USA
> >
> >
> > _______________________________________________
> > ncl-talk mailing list
> > ncl-talk at mailman.ucar.edu
> > List instructions, subscriber options, unsubscribe:
> > https://mailman.ucar.edu/mailman/listinfo/ncl-talk
> >
>
>
> --
> Dr. Anil Kumar
> Environmental Modeling Center/NCEP/NWS/NOAA
> College Park, MD, USA
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: <https://mailman.ucar.edu/pipermail/ncl-talk/
> attachments/20201204/f6a4a763/attachment-0001.html>
>
> ------------------------------
>
> Message: 8
> Date: Fri, 4 Dec 2020 13:05:47 -0500
> From: Anil Kumar - NOAA Affiliate <anil.kumar at noaa.gov>
> To: Ehsan Taghizadeh <ehsantaghizadeh at yahoo.com>
> Cc: Ncl-talk <ncl-talk at ucar.edu>
> Subject: Re: [ncl-talk] wrfout file reading
> Message-ID:
> <CA+8gui8kW-L=eBdVhSv6BK5fXgdanjfRg_118O+TRXUozYQ8UQ at mail.
> gmail.com>
> Content-Type: text/plain; charset="utf-8"
>
> Yes, Hourly wrf outputs.
>
> and yes when I extract times using *wrf_user_list_times*(f1) that works
> fine but only times.
> but how to extract variables like U10 time series?
>
> Thanks,
> Anil
>
>
>
> On Fri, Dec 4, 2020 at 12:42 PM Ehsan Taghizadeh <
> ehsantaghizadeh at yahoo.com>
> wrote:
>
> > Dear Dr. Anil
> > Are your wrfout files hourly?
> > However, why you don't use wrf_user_list_times to extract the list of
> > available times?
> >
> > fs1 = systemfunc("ls ../WRFV3/run/wrfout_d01*")
> > nfs1= dimsizes(fs1)
> > do i = 0, nfs1-1
> > f1 = addfile(fs1(i)+".nc","r")
> > times = *wrf_user_list_times*(f1) ; get times in the file
> >
> > ntimes = dimsizes <https://www.ncl.ucar.edu/
> Document/Functions/Built-in/dimsizes.shtml>(times) ; number of
> times in the file(s)
> > end do
> >
> > *Sincerely*
> > *Ehsan*
> >
> > On Friday, December 4, 2020, 05:57:54 PM GMT+3:30, Anil Kumar - NOAA
> > Affiliate via ncl-talk <ncl-talk at mailman.ucar.edu> wrote:
> >
> >
> > Hi,
> > I am trying to read multiple wrfout files. I am reading okay..no problem
> > with the code extracting variables.
> > Problem comes, when each file has multiple times in it and then I don't
> > know how to extract the time series of variables .
> > Below is the code, please let me know how to modify code to read many
> > times within each file.
> > ************************
> > fs1 = systemfunc("ls ../WRFV3/run/wrfout_d01*")
> > nfs1= dimsizes(fs1)
> > f1 = addfiles(fs1+".nc","r")
> > ntimes1 = dimsizes(f1) ; number of times in the file(s)
> >
> > loc1 = wrf_user_ij_to_ll(f1, 593., 546., True)
> > print("lon/lat locations are: " + loc)
> > uu1 = new((/ntimes1/),float)
> > vv1 = new((/ntimes1/),float)
> > psfc1 = new((/ntimes1/),float)
> > i1=0
> > do while (i1.le.ntimes1-1)
> > uu1(i1) = f1[i1]->U(:,0,546,593) ; (time, lat,lon) ; (2,64,128)
> > vv1(i1) = f1[i1]->V(:,0,546,593) ; V10(:,29.301,-94.797)
> > psfc1(i1) = f1[i1]->PSFC(:,546,593)
> > i1 = i1 + 1
> > end do
> > model_wspd1 = sqrt(uu1^2+vv1^2) ;wind_speed(uu,vv)
> > model_wdir1 = wind_direction(uu1,vv1,0)
> > model_psfc1 = psfc1/100.
> > ************************
> > --
> > Dr. Anil Kumar
> > Environmental Modeling Center/NCEP/NWS/NOAA
> > College Park, MD, USA
> >
> >
> > _______________________________________________
> > ncl-talk mailing list
> > ncl-talk at mailman.ucar.edu
> > List instructions, subscriber options, unsubscribe:
> > https://mailman.ucar.edu/mailman/listinfo/ncl-talk
> >
>
>
> --
> Dr. Anil Kumar
> Environmental Modeling Center/NCEP/NWS/NOAA
> College Park, MD, USA
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: <https://mailman.ucar.edu/pipermail/ncl-talk/
> attachments/20201204/4a7be921/attachment-0001.html>
>
> ------------------------------
>
> Subject: Digest Footer
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at mailman.ucar.edu
> List instructions, subscriber options, unsubscribe:
> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
> ------------------------------
>
> End of ncl-talk Digest, Vol 205, Issue 5
> ****************************************
>
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20201205/1548d8ce/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: 6B57029C99184E59976BE6606EA02A83.png
Type: image/png
Size: 827 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20201205/1548d8ce/attachment.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: 5070BC9563F448D3BEE04BF94C48567B.png
Type: image/png
Size: 805 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20201205/1548d8ce/attachment-0001.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: 9C93A62AE27A49CDA4E88C4EC7C77C5C.png
Type: image/png
Size: 920 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20201205/1548d8ce/attachment-0002.png>
More information about the ncl-talk
mailing list