[ncl-talk] error running NCL : UTF8 Locale not supported
Mary Haley
haley at ucar.edu
Thu Jun 4 10:02:12 MDT 2015
Heri,
In the future, please post installation issues to ncl-install and not
ncl-talk.
Also, please don't include a bunch of digest email from other postings,
because it becomes hard to read.
You said you couldn't run any NCL commands, but you are running nedit"
which is not an NCL command. This is a UNIX editor, and is developed and
supported by somebody else.
If you google:
"nedit UTF8 locale not supported"
You will get lots of hits with potential solutions.
If you are having problems running NCL, then please include the command
that you are typing and the error message in a new email post, so we can
read it better.
Thanks,
--Mary
On Thu, Jun 4, 2015 at 2:48 AM, Heri Kuswanto <kuswanto.its82 at yahoo.com>
wrote:
> Dear all,
>
> I tried to install ncl in my WIndows system. It works but when I try to
> run it by writing "nedit" it give the following message " UTF8 locale not
> supported". Consequently, I cannot run any *.ncl command. The detail is as
> follows:
>
> $ startxwin
>
> --------------------------------------------------------------
> Could not load crashreporter dll
> Welcome to the XWin X Server
> Vendor: The Cygwin/X Project
> Release: 1.17.1.0
> OS: CYGWIN_NT-6.1 ROSYIDI-PC 2.0.1(0.287/5/3) 2015-04-30 18:13 i686
> OS: Windows 7 [Windows NT 6.1 build 7600] (Win32)
> Package: version 1.17.1-4 built 2015-05-01
>
> XWin was started with the following command line:
>
> /usr/bin/XWin :0 -multiwindow -auth
> /home/ROSYIDI/.serverauth.5928
>
> Could not load crashreporter dll
> (II) xorg.conf is not supported
> (II) See http://x.cygwin.com/docs/faq/cygwin-x-faq.html for more
> information
> LoadPreferences: /home/ROSYIDI/.XWinrc not found
> LoadPreferences: Loading /etc/X11/system.XWinrc
> LoadPreferences: Done parsing the configuration file...
> winDetectSupportedEngines - DirectDraw4 installed, allowing ShadowDDNL
> winDetectSupportedEngines - Returning, supported engines 00000005
> winSetEngine - Multi Window or Rootless => ShadowGDI
> winScreenInit - Using Windows display depth of 32 bits per pixel
> winAllocateFBShadowGDI - Creating DIB with width: 1366 height: 768 depth:
> 32
> winFinishScreenInitFB - Masks: 00ff0000 0000ff00 000000ff
> winInitVisualsShadowGDI - Masks 00ff0000 0000ff00 000000ff BPRGB 8 d 24
> bpp 32
> MIT-SHM extension disabled due to lack of kernel support
> XFree86-Bigfont extension local-client optimization disabled due to lack
> of shared memory support in the kernel
> glWinSelectGLimplementation: Loaded 'cygnativeGLthunk.dll'
> (II) AIGLX: Testing pixelFormatIndex 1
> GL_VERSION: 4.1.10825 Compatibility Profile Context
> GL_VENDOR: ATI Technologies Inc.
> GL_RENDERER: AMD Radeon HD 6320 Graphics
> (II) AIGLX: enabled GLX_SGI_make_current_read
> (II) AIGLX: enabled GLX_MESA_copy_sub_buffer
> (II) AIGLX: enabled GLX_SGI_swap_control and GLX_MESA_swap_control
> (II) AIGLX: enabled GLX_SGIX_pbuffer
> (II) AIGLX: enabled GLX_ARB_multisample and GLX_SGIS_multisample
> (II) 95 pixel formats reported by wglGetPixelFormatAttribivARB
> (II) AIGLX: Set GLX version to 1.4
> (II) 24 fbConfigs
> (II) ignored pixel formats: 0 not OpenGL, 12 RBGA float, 24 RGBA unsigned
> float, 0 unknown pixel type, 35 unaccelerated
> (II) GLX: Initialized Win32 native WGL GL provider for screen 0
> winPointerWarpCursor - Discarding first warp: 683 384
> (--) 5 mouse buttons found
> (--) Setting autorepeat to delay=500, rate=31
> (--) Windows keyboard layout: "00000409" (00000409) "US", type 7
> (--) Found matching XKB configuration "English (USA)"
> (--) Model = "pc105" Layout = "us" Variant = "none" Options = "none"
> Rules = "base" Model = "pc105" Layout = "us" Variant = "none" Options =
> "none"
> winInitMultiWindowWM - DISPLAY=:0.0
> winMultiWindowXMsgProc - DISPLAY=:0.0
> winProcEstablishConnection - winInitClipboard returned.
> winClipboardThreadProc - DISPLAY=:0.0
> OS maintains clipboard viewer chain: yes
> winInitMultiWindowWM - XOpenDisplay () returned and successfully opened
> the display.
> winMultiWindowXMsgProc - XOpenDisplay () returned and successfully opened
> the display.
> winClipboardProc - XOpenDisplay () returned and successfully opened the
> display.
> winMultiWindowXMsgProcErrorHandler - ERROR: BadMatch (invalid parameter
> attributes)
>
> ---------------------------------------------------------------------------------------------------------------------------------------
> $nedit
>
> UTF8 locale not supported.
>
> --------------
> Anyone know the solution, please?
>
> Thanks for your help.
>
> Heri
>
>
>
>
> On Thursday, June 4, 2015 1:00 AM, "ncl-talk-request at ucar.edu" <
> 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. Re: exporting variables and dimensions in excel sheet
> (Rick Brownrigg)
> 2. Re: Editing a netcdf4 file (David Brown)
> 3. Re: fileattdef seg fault on CESM LE files (David Brown)
> 4. Re: drawing RMS semi-circles in Taylor diagram (Alan Brammer)
> 5. Re: functions handling input of unknown dimension (Dennis Shea)
> 6. Re: Panel Plot (Dennis Shea)
>
>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Wed, 3 Jun 2015 08:07:07 -0600
> From: Rick Brownrigg <brownrig at ucar.edu>
> Subject: Re: [ncl-talk] exporting variables and dimensions in excel
> sheet
> To: amit timilsina <timilsinaamit87 at huskers.unl.edu>
> Cc: "ncl-talk at ucar.edu" <ncl-talk at ucar.edu>
> Message-ID:
> <CAGKRhbEKr8MqfHvsaNvxCWWy8FaqA2EyTUqCUXu3Q2Jh7oOB7A at mail.gmail.com>
> Content-Type: text/plain; charset="utf-8"
>
> Hi,
>
> Are you wanting to import these data into an Excel spreadsheet from the
> NetCDF file? If using NCL for this purpose, you likely want to write a
> comma-separated-file (csv), which Excel could directly import. Check out
> the ascii-output examples at:
>
> http://ncl.ucar.edu/Applications/o-ascii.shtml
>
> particularly the write_table function.
>
> Rick
>
> On Wed, Jun 3, 2015 at 7:58 AM, amit timilsina <
> timilsinaamit87 at huskers.unl.edu> wrote:
>
> >
> > Hello all,
> > The attached file contains tasmin_diff variable, and month,lat and lon
> > dimensions. Is it possible to export these all in combination in excel
> > sheet? Please share link any procedure described. Thanks in advance
> >
> >
> > Regards,
> > Amit
> > _______________________________________________
> > 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/20150603/c6682ba7/attachment-0001.html
>
>
> ------------------------------
>
> Message: 2
> Date: Wed, 3 Jun 2015 08:50:48 -0600
> From: David Brown <dbrown at ucar.edu>
> Subject: Re: [ncl-talk] Editing a netcdf4 file
> To: Gerry Creager - NOAA Affiliate <gerry.creager at noaa.gov>
> Cc: ncl-talk <ncl-talk at ucar.edu>
> Message-ID:
> <CAC92F7fVj6XAh_MEMGiS24wMMK9v13Pw-3piOcz6rvaVAofn6g at mail.gmail.com>
> Content-Type: text/plain; charset=UTF-8
>
> Hi Dustin and Gerry,
> I have confirmed there is a problem. It has been partially but not
> completely fixed in the development version of the code. I am looking
> into the remaining issue and will let you know when it has been fixed.
> -dave
>
> On Tue, Jun 2, 2015 at 10:34 AM, Gerry Creager - NOAA Affiliate
> <gerry.creager at noaa.gov> wrote:
> > Additional information: If we do an addfile("<filename>","r"), we can
> open
> > and read the file. If we use "w", well, we can open the file but
> assigning
> > an existing variable to an internal variable fails on the assign (var =
> > a->QRAIN). Of course, because it fails there, we can't see if it'd also
> fail
> > on write.
> >
> > Thanks
> > Gerry
> >
> > On Mon, Jun 1, 2015 at 7:15 PM, Gerry Creager - NOAA Affiliate
> > <gerry.creager at noaa.gov> wrote:
> >>
> >> I've started looking at this with Dusty. Appreciate any thought you
> might
> >> have as to why we're failing.
> >>
> >> I've uploaded the input file (wrfinput_d02) to the FTP incoming site.
> I'm
> >> testing with the script:
> >>
> >> var = "QRAIN"
> >> print(var)
> >> a=addfile("wrfinput_d02.nc", "w")
> >> var_orig=a->QRAIN
> >> print(var_orig)
> >>
> >>
> >> The result looks like:gerry at loki:~/NEWSe/NCL> ncl test.ncl
> >> Copyright (C) 1995-2015 - All Rights Reserved
> >> University Corporation for Atmospheric Research
> >> NCAR Command Language Version 6.3.0
> >> The use of this software is governed by a License Agreement.
> >> See http://www.ncl.ucar.edu/ for more details.
> >>
> >>
> >> Variable: var
> >> Type: string
> >> Total Size: 8 bytes
> >> 1 values
> >> Number of Dimensions: 1
> >> Dimensions and sizes: [1]
> >> Coordinates:
> >> (0) QRAIN
> >> fatal:["NclNetCDF4.c":3073]:NclNetCDF4.c: Error in nc_def_var_chunking
> in
> >> file (wrfinput_d02.nc) for writing, at line: 3073
> >>
> >> line 3074 of NclNetCDF4.c: NetCDF: Attempt to define var properties,
> like
> >> deflate, after enddef.
> >>
> >>
> >> On Fri, May 29, 2015 at 8:54 PM, Dustan Wheatley - NOAA Affiliate
> >> <dustan.wheatley at noaa.gov> wrote:
> >>>
> >>> I am trying to overwrite a netcdf4 file analogous to what is shown at
> >>> http://www.ncl.ucar.edu/Applications/change_netCDF.shtml. In the
> attached
> >>> ncl script, I open the file to be modified at line 35 and attempt the
> >>> modification at line 73 (with a lot of calculations in between).
> While I
> >>> have double-checked that the modified field (var_mod) has the same
> >>> dimensionality as the original variable (QRAIN), I am still get a
> >>> segmentation fault/core dump. Never had this problem when I tried to
> >>> manipulate netcdf3 files. Could it possibly be related to the netcdf4
> >>> compression? It?s possible I may be trying to do something that?s no
> longer
> >>> possible.
> >>>
> >>> Thanks in advance for any suggestions.
> >>>
> >>> -Dusty Wheatley
> >>>
> >>> _______________________________________________
> >>> ncl-talk mailing list
> >>> ncl-talk at ucar.edu
> >>> List instructions, subscriber options, unsubscribe:
> >>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
> >>>
> >>
> >>
> >>
> >> --
> >> Gerry Creager
> >> NSSL/CIMMS
> >> 405.325.6371
> >> ++++++++++++++++++++++
> >> ?Big whorls have little whorls,
> >> That feed on their velocity;
> >> And little whorls have lesser whorls,
> >> And so on to viscosity.?
> >> Lewis Fry Richardson (1881-1953)
> >
> >
> >
> >
> > --
> > Gerry Creager
> > NSSL/CIMMS
> > 405.325.6371
> > ++++++++++++++++++++++
> > ?Big whorls have little whorls,
> > That feed on their velocity;
> > And little whorls have lesser whorls,
> > And so on to viscosity.?
> > Lewis Fry Richardson (1881-1953)
> >
> > _______________________________________________
> > ncl-talk mailing list
> > ncl-talk at ucar.edu
> > List instructions, subscriber options, unsubscribe:
> > http://mailman.ucar.edu/mailman/listinfo/ncl-talk
> >
>
>
> ------------------------------
>
> Message: 3
> Date: Wed, 3 Jun 2015 09:36:55 -0600
> From: David Brown <dbrown at ucar.edu>
> Subject: Re: [ncl-talk] fileattdef seg fault on CESM LE files
> To: "David B. Reusch" <dreusch at ees.nmt.edu>
> Cc: NCL Talk <ncl-talk at ucar.edu>
> Message-ID:
> <CAC92F7dnPLD=J4_Odjydw7hkALzbn1vGoLXUHGGxhRwNpMs4CA at mail.gmail.com>
> Content-Type: text/plain; charset=UTF-8
>
> Hi David,
> I want to let you know that this issue has been fixed in the
> development version of NCL. Let us know if you'd like to try it out on
> yellowstone.
> -dave
>
> On Tue, May 26, 2015 at 2:11 PM, Mary Haley <haley at ucar.edu> wrote:
> > Hi Dave,
> >
> > Thanks for providing the script.
> >
> > I was able to reproduce your results with NCL 6.2.1, NCL 6.3.0 (both
> intel
> > and gcc versions) and a test version of NCL, which is our development
> > version. I've created a ticket, NCL-2223.
> >
> > If/when we get this working, I can point you to a development version
> that
> > you can use on yellowstone.
> >
> > --Mary
> >
> >
> >
> > On Tue, May 26, 2015 at 11:37 AM, David B. Reusch <dreusch at ees.nmt.edu>
> > wrote:
> >>
> >> I seem to have found an NCL seg fault when copying global attributes
> with
> >> fileattdef from CESM Large Ensemble (LE) files. This was part of a
> larger
> >> script I've been using to regrid for multiple years. The attached script
> >> just does the fileattdef step and it works on at least one of my other,
> >> non-CESM LE files. I've tried it on multiple CESM LE files with the
> same
> >> results. Copying the attributes individually works so I'm assuming it's
> >> related to doing all the attributes together.
> >>
> >> This is on yellowstone with NCL 6.3.
> >>
> >> My temporary workaround is to just comment out this step so that I can
> >> continue my file processing. Since that means losing global metadata, a
> >> better solution/fix would be preferred.
> >>
> >> Thanks,
> >> Dave
> >>
> >> --
> >> Associate Research Professor of Climatology
> >> Dept of Earth and Environmental Science
> >> MSEC 304; 801 Leroy Place
> >> New Mexico Tech
> >> Socorro, NM 87801
> >>
> >>
> >> _______________________________________________
> >> ncl-talk mailing list
> >> List instructions, subscriber options, unsubscribe:
> >> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
> >>
> >
> >
> > _______________________________________________
> > ncl-talk mailing list
> > List instructions, subscriber options, unsubscribe:
> > http://mailman.ucar.edu/mailman/listinfo/ncl-talk
> >
>
>
> ------------------------------
>
> Message: 4
> Date: Wed, 3 Jun 2015 11:40:14 -0400
> From: Alan Brammer <abrammer at albany.edu>
> Subject: Re: [ncl-talk] drawing RMS semi-circles in Taylor diagram
> To: LI Qi <liqi123sh at qq.com>
> Cc: ncl-talk <ncl-talk at ucar.edu>
> Message-ID:
> <CA+V766SnYssWRYVuoSJeQcyZz7MBYivwjTcKzBwO_=tz67qnCg at mail.gmail.com>
> Content-Type: text/plain; charset="utf-8"
>
> I would say you want to edit the constants on the below lines. I'm not
> sure how these were obtained though, trial and error might be the simplest
> way.
> These are just restricting how long the line is, with your base plot I
> would think that even your 3rd circle can span the full 180-360. Increase
> 0.77 and 0.61 until the circles reach a respective edge.
> n3 = floattointeger( 0.77*npts )
> n4 = floattointeger( 0.61*npts )
>
>
> Good luck.
>
> On 3 June 2015 at 09:43, LI Qi <liqi123sh at qq.com> wrote:
>
> > Dear NCL developers,
> >
> > I use the script "taylor_diagram.ncl" to draw Taylor diagram. I extend
> the
> > outer curve from the original 1.65 to 1.75 to fit my data, but it seems
> > that the parameters w.r.t. RMS semi-circles need also be changed
> > accordingly (as is shown below, the semi-circles are incomplete).
> > The issue is, I have no idea how to change it.
> >
> >
> > Below is a snippet of the script, which I guess is related to the RMS
> > semi-circles
> > Any hints about how to change the parameters would be sincerely
> > appreciated!
> >
> > Best,
> > Qi
> >
> > dx = 0.25
> > ncon = 4 ; 0.75, 0.50, 0.25,
> > 0.0
> > npts = 100 ; arbitrary
> > ang = fspan(180,360,npts)*rad
> >
> > dum9 = new(ncon,graphic)
> >
> > do n=1,ncon
> > rr = n*dx ; radius from 1.0 [OBS] abscissa
> > xx = 1. + rr*cos(ang)
> > yy = fabs( rr*sin(ang) )
> > if (n.le.2) then
> > dum9(n-1) = gsn_add_polyline(wks,taylor,xx,yy,respl)
> > end if
> > if (n.eq.3) then
> > n3 = floattointeger( 0.77*npts )
> > dum9(n-1) =
> > gsn_add_polyline(wks,taylor,xx(0:n3),yy(0:n3),respl)
> > end if
> > if (n.eq.4) then
> > n4 = floattointeger( 0.61*npts )
> > dum9(n-1) =
> > gsn_add_polyline(wks,taylor,xx(0:n4),yy(0:n4),respl)
> > end if
> > end do
> >
> > _______________________________________________
> > 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/20150603/4b3ede87/attachment-0001.html
>
> -------------- next part --------------
> A non-text attachment was scrubbed...
> Name: 2E05F4FE at 0FC9496C.9F046F55
> Type: application/octet-stream
> Size: 56453 bytes
> Desc: not available
> Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20150603/4b3ede87/attachment-0001.obj
>
>
> ------------------------------
>
> Message: 5
> Date: Wed, 3 Jun 2015 10:12:41 -0600
> From: Dennis Shea <shea at ucar.edu>
> Subject: Re: [ncl-talk] functions handling input of unknown dimension
> To: Walter Kolczynski <walter.kolczynski at noaa.gov>
> Cc: "ncl-talk at ucar.edu" <ncl-talk at ucar.edu>
> Message-ID:
> <CAOF1d_6MzNWNh_-eycynEZvbehBF8Z2SjJjsOAv6ND-kRbx_Lw at mail.gmail.com>
> Content-Type: text/plain; charset="utf-8"
>
> Will, Walt
>
> To my knowledge there is no 'simple' way to do what Will originally posted.
>
> Walt's 'flattening' suggestion would be 'clean' if it is general and hidden
> with a simple interface. However, I am not sure that it would be more
> efficient. More memory would be required for the 'flat' array.
> --
>
> re: WH's comment: "particularly 'clean' or efficient code"
>
> The nested 'if' approach is certainly cumbersome
>
> if (rankx.eq.1) then
> ...
> end if
> if (rankx.eq.2) then
> ...
> end if
> if (rankx.eq.3) then
> ...
> end if
> if (rankx.eq.4) then
> ...
> end if
>
> but I don't think it is particularly inefficient.
>
>
> On Tue, Jun 2, 2015 at 11:09 AM, Walter Kolczynski <
> walter.kolczynski at noaa.gov> wrote:
>
> > Will,
> >
> > Depending on what you are doing, you might be able to flatten into a
> lower
> > dimensional array, do whatever operation, then convert back. I've
> attached
> > an example to get you started: it returns a random sample in the
> left-most
> > dimension of any 2+ dimensional array. To run it just go into an
> > interactive NCL session, load the file, and run the test procedure:
> >
> > [user at machine ~]$ ncl
> > Copyright (C) 1995-2014 - All Rights Reserved
> > University Corporation for Atmospheric Research
> > NCAR Command Language Version 6.2.0
> > The use of this software is governed by a License Agreement.
> > See http://www.ncl.ucar.edu/ for more details.
> > ncl 0> load "takeSample.ncl"
> > ncl 1> test_takeSample()
> >
> > - Walter
> >
> >
> > On 02-Jun-15 02:43, Will Hobbs wrote:
> >
> > Hi all
> >
> > I wonder if anyone has any suggestions for a generic way of writing
> > functions that deal with input/output arrays which may have different
> > numbers of dimensions, e.g.
> >
> > function generic_example(input:numeric)
> >
> > begin
> > output = some_subset_of_(input)
> >
> > return(output)
> > end
> >
> > where input/output might be 1, 2, 3, 4 or more dimensions, and so a
> > different line is needed to do the subsetting step.
> >
> > At the moment I just use nested 'if' statements based on the rank of the
> > array (where rank is the number of dimensions, i.e. dimsizes(dim
> > sizes(input)), and this works fine. But, it does not make for
> particularly
> > 'clean' or efficient code, so I wonder if anyone has any nicer ways of
> > doing this.
> >
> > Many thanks
> >
> > Will
> >
> >
> >
> > University of Tasmania Electronic Communications Policy (December, 2014).
> > This email is confidential, and is for the intended recipient only.
> > Access, disclosure, copying, distribution, or reliance on any of it by
> > anyone outside the intended recipient organisation is prohibited and may
> be
> > a criminal offence. Please delete if obtained in error and email
> > confirmation to the sender. The views expressed in this email are not
> > necessarily the views of the University of Tasmania, unless clearly
> > intended otherwise.
> >
> >
> > _______________________________________________
> > ncl-talk mailing listncl-talk at ucar.edu
> > List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
> >
> >
> > --
> > Walter Kolczynski, Jr.
> > Global Ensemble Team
> > NOAA/NWS/NCEP/EMC (via I.M. Systems Group)(301) 683-3781
> >
> >
> > _______________________________________________
> > 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/20150603/92a38701/attachment-0001.html
>
>
> ------------------------------
>
> Message: 6
> Date: Wed, 3 Jun 2015 11:33:33 -0600
> From: Dennis Shea <shea at ucar.edu>
> Subject: Re: [ncl-talk] Panel Plot
> To: Amadou Coulibaly <mpapin24 at gmail.com>
> Cc: "ncl-talk at ucar.edu" <ncl-talk at ucar.edu>
> Message-ID:
> <CAOF1d_62pZAVQ=MtHYNicnQdtdatfet_kKH7yee363+wkTM_+w at mail.gmail.com>
> Content-Type: text/plain; charset="utf-8"
>
> Looking at scripts can be a tedious and time consuming task.
> --
>
> Your code
>
> infile = addfile("/home/amadou/NCL/all_new.nc", "r")
> temp = infile->air_temperature
> Temp = temp(256392:258575,7) ; for 5 years
> (2012_JFMAMJ)
> time = infile->time
> Time = time(256392:258575)
>
> index = ind(Temp(:) .ne. Temp at _FillValue) ; indices of
> Temp without missing values
> printVarSummary(index)
> Time0 = time(256444:258575)
> Temp1 = new((/2132/), float)
> Temp1 = Temp(index)
>
> [snip]
>
> [1] You can **not** use the following to check for missing values. These
> can *only* be checked via the 'ismissing' function.
>
> index = ind(Temp(:) .ne. Temp at _FillValue) ; indices
> of Temp without missing values
>
> Use
>
> index = ind(.not.ismissing(Temp)) ; indices of Temp without missing
> values
>
> https://www.ncl.ucar.edu/Document/Functions/Built-in/ismissing.shtml
> https://www.ncl.ucar.edu/Document/Functions/Built-in/ind.shtml
>
> [2]
> As I mentioned to you in a previous ncl-talk response, you can not just
> eliminate missing values and perform a wavelet analysis on the 'good'
> values. Omitting values will result in 'numbers' but they are not the
> 'correct numbers'. Omitting values results in spectral gaps. If you have
> (say) hourly values and very short spans with missing values you could
> perform a linear interpolation. This will affect the spectra but likely not
> in any substantive manner.
>
> [3]
> You have two different sets of index values: 256392:258575 and
> 256444:258575 ... this is likely the source of your plotting issue.
>
> Good programming practice suggest that you use variables rather than hard
> wired constants.
>
> tStrt = 256392
> tLast = 258575
>
> and use these
>
> temp = infile->air_temperature
> Temp = temp(tStrt:tLast,7) ; for 5 years
> (2012_JFMAMJ)
> time = infile->time
> Time = time(tStrt:tLast)
>
> Actually .... better
>
> Temp = infile->air_temperature(tStrt:tLast,7)
> printVarSummary(Temp)
>
>
>
>
>
>
>
>
>
> On Wed, Jun 3, 2015 at 6:00 AM, Amadou Coulibaly <mpapin24 at gmail.com>
> wrote:
>
> > Dear NCL users,
> >
> > I tried to do a panel plot between "Wavelet Power Spectrum" and
> "Timesires"
> > from hourly temperature data of April-May-June. But this data contained
> 61
> > hours of missing values that I cancelled. When I plotted my data, I
> > remarked that the x-axis on the wavelet power spectrum does not match the
> > timeserie one.
> >
> > What is wrong with my script?
> >
> > Please find attached my script and my plot.
> >
> > Best
> > --
> > * COULIBALY AMADOU *
> > PhD Student on West African Climate System (WACS)
> > FUTA - Federal University of Technology of Akure, Nigeria
> >
> > *Visiting Student - *
> > *University of Cologne, Germany**Institute of Geophysics and Meteorology*
> > Pohligstr. 3 / Office 3.102
> > D-50969 K?ln
> >
> > *Project*: WASCAL (West African Science Service Centre on Climate Change
> > and Adapted Land Use)
> > Phone:(+234) 810 795 2836 /(+223) 65 67 27 57 /(+226) 64 57 37 27 /+49
> > 15218352574
> >
> > E-mail: mpapin24 at gmail.com / coulibalya68 at yahoo.com
> >
> > "*The time is always right to do right": Nelson Mandela*
> >
> > *"Character is like a tree and reputation like a shadow. The shadow is
> > what we think of it; the tree is the real thing" : Abraham Lincoln*
> >
> > *"Do what you can, with what you have, where you are" Theodore Roosevelt*
> >
> >
> >
> >
> >
> >
> > _______________________________________________
> > 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/20150603/2ea068de/attachment-0001.html
>
>
> ------------------------------
>
> _______________________________________________
> 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 139, Issue 5
> ****************************************
>
>
>
> _______________________________________________
> 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/20150604/ff7850a4/attachment.html
More information about the ncl-talk
mailing list