Re: Make the spacing between the tick marks equal without changing the explicit tick mark values

Thanks Mary
Yes, I want the plot similar to example1
(axes_1.ncl<http://www.ncl.ucar.edu/Applications/Scripts/axes_1.ncl>)..But
when I use *res_at_gsnYAxisIrregular2Linear*<http://www.ncl.ucar.edu/Document/Graphics/Resources/gsn.shtml#gsnYAxisIrregular2Linear>=True
it doesn't make any changes to my plot. I think this resourse will work only
if my Y axis coordinate arrays are not linear. But i my case I think the Y
axis coordinate arrays are linear and I am explicitly defining the tick mark
values i.e (/0.01, 0.0125, 0.016666667, 0.025,0.05, 0.1/)..Note that
interval between the tick mark increases as we go from left to right in the
given values.(bottom to top of Y axis). The print(WORK2&freq) gives the
following values.

Variable: freq (coordinate)
Type: float
Total Size: 60 bytes
            15 values
Number of Dimensions: 1
Dimensions and sizes: [freq | 15]
Coordinates:
(0) 0
(1) 0.008196721
(2) 0.01639344
(3) 0.02459016
(4) 0.03278688
(5) 0.04098361
(6) 0.04918033
(7) 0.05737705
(8) 0.06557377
(9) 0.07377049
(10) 0.08196721
(11) 0.09016393
(12) 0.09836066
(13) 0.1065574
(14) 0.1147541

Any help would be highly appreciated
sabeer

On 2 July 2010 23:30, <ncl-talk-request@ucar.edu> wrote:

> Send ncl-talk mailing list submissions to
> ncl-talk@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@ucar.edu
>
> You can reach the person managing the list at
> ncl-talk-owner@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: RAMS HDF5 format? (Dreher, Joseph Civ USAF AFTAC AFTAC/TMNV)
> 2. Re: How can I use equivalent latitudes? (Dennis Shea)
> 3. Re: Adjusting tickmark labels with gsn_csm_contour_map
> (Mary Haley)
> 4. overplot (WuDraco)
> 5. Re: overplot (debasish@sca.uqam.ca)
> 6. how to calculate POP (Principal Oscillation Pattern) with NCL
> (???)
> 7. noise after 2D interpolation (Li, Yongzuo)
> 8. Make the spacing between the tick marks equal without
> changing the explicit tick mark values. (Sabeerali(sebi))
> 9. Re: RAMS HDF5 format? (Mary Haley)
> 10. Re: Make the spacing between the tick marks equal without
> changing the explicit tick mark values. (Mary Haley)
> 11. Re: noise after 2D interpolation (Mary Haley)
> 12. lag composites (Straub, Katherine)
> 13. Re: noise after 2D interpolation (Dennis Shea)
> 14. Re: lag composites (Dennis Shea)
> 15. invalid pointer error message (Ioana Colfescu)
> 16. v5.2.1 for rcm2points to solve interpolation noise (Li, Yongzuo)
> 17. Re: invalid pointer error message (Dennis Shea)
>
>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Thu, 1 Jul 2010 15:14:54 -0400
> From: "Dreher, Joseph Civ USAF AFTAC AFTAC/TMNV"
> <Joseph.Dreher@patrick.af.mil>
> Subject: Re: RAMS HDF5 format?
> To: "Dennis Shea" <shea@ucar.edu>
> Cc: ncl-talk@ucar.edu
> Message-ID:
> <
> A57E54C24C2801428660457002D1109307022CA2@ENPAEX04.aftac.patrick.afspc.ds.af.mil
> >
>
> Content-Type: text/plain; charset="us-ascii"
>
> Thanks for the information. I will try uploading several sample files
> early next week.
>
> Joe
>
> -----Original Message-----
> From: Dennis Shea [mailto:shea@ucar.edu]
> Sent: Thursday, July 01, 2010 12:34 PM
> To: Dreher, Joseph Civ USAF AFTAC AFTAC/TMNV
> Cc: ncl-talk@ucar.edu
> Subject: Re: RAMS HDF5 format?
>
> Hello,
>
> I believe that HDF5 will be in v6.0.0.
>
> There will be a 5.2.1 release before that.
> I am not sure how much of HDF5 will be supported in that release.
>
> Can you please send a sample RAMS HDF5 file?
> The NCL team can use it for testing. If NCL's development branch
> reads the file successfully, I am pretty sure a beta version could
> be created for your use.
>
> [1] Please send the output from:
>
> %> uname -a
>
> [2]
>
> ftp ftp.cgd.ucar.edu
> anonymous
> Use your email address for the password
> cd incoming
> put ...
> put ...
> :
> quit
>
> Send me an email with the following information:
>
> * The names of the files after you have successfully completed the
> transfer.
> We cannot look at the ftp/incoming directory for security
> reasons.
>
>
> Cheers
>
>
> On 7/1/10 8:57 AM, Dreher, Joseph Civ USAF AFTAC AFTAC/TMNV wrote:
> > Hello all,
> >
> > I recently attended the WRF user's workshop and saw the overview on
> NCL.
> > I believe I saw a note that the next release is supposed to include
> > readers for HDF5 format? In my organization we currently are running
> > RAMS and now WRF and would like to have the capability to read both
> the
> > RAMS HDF5 and WRF netCDF file formats and compare with NCL. Is it
> still
> > true that the next scheduled release will include HDF5?
> >
> > Thanks in advance,
> >
> > Joe
> >
> >
> >
> > _______________________________________________
> > ncl-talk mailing list
> > List instructions, subscriber options, unsubscribe:
> > http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
> ------------------------------
>
> Message: 2
> Date: Thu, 01 Jul 2010 16:55:29 -0600
> From: Dennis Shea <shea@ucar.edu>
> Subject: Re: How can I use equivalent latitudes?
> To: "Runde, Theresa" <Theresa.Runde@dlr.de>
> Cc: ncl-talk@ucar.edu
> Message-ID: <4C2D1CE1.2000107@ucar.edu>
> Content-Type: text/plain; charset="iso-8859-1"
>
> Hi Theresa,
>
> Attached is a sample script that will plot the data
> using the conventional latitudes and the equivalent latitudes.
> The key is the explicitly create the sfXArray and sfYArray
> resources and convert the corresponding data array to one-dimensional.
>
> diri = "./"
> fili = "PV_Eqlat_600K_22503_2403.nc"
> f = addfile(diri+fili, "r")
>
> x = f->vortpot_0600K(nt,:,:) ; (48,64)
> elat = f->eqlat_0600K(nt,:,:) ; (48,64)
> elon = conform_dims(dimsizes(elat), x&lon, 1)
>
> res@sfXArray = ndtooned(elon)
> res@sfYArray = ndtooned(elat)
> x1d = ndtooned(x)
>
> res@gsnLeftString = "Equivalent Latitudes"
> plot = gsn_csm_contour_map_polar(wks,x1d,res)
>
> Cheers
> D
>
> On 6/29/10 8:16 AM, Runde, Theresa wrote:
> > Hello,
> >
> > I have a big problem with handling equivalent latitudes and using them
> > on a regular grid for plotting. I calculated equivalent latitudes and
> > now I have an equivalent latitude array of equlat=equlat(time,lat,lon)
> > for each potential temperature level. That means: my new (equivalent)
> > latitudes are 3 dimensional.
> >
> > Now I want to plot (for example) potential vorticity (PV) dependent on
> > these new equivalent latitudes (to get a zonally symmetric polar vortex).
> > At the moment my PV field is also 3 dimensional (time x lat x lon), but
> > regular (latitudes and longitudes are always the same). How can I
> > transform my PV field (time x lat(1d) x lon) to a PV field (time x
> > equivalent lats(3d) x lon) and plot this? The time dependence is
> necessary.
> >
> > Thank you for your help,
> > Theresa Runde
> > _______________________________________________
> > ncl-talk mailing list
> > List instructions, subscriber options, unsubscribe:
> > http://mailman.ucar.edu/mailman/listinfo/ncl-talk
> -------------- next part --------------
> An embedded and charset-unspecified text was scrubbed...
> Name: equivLat.ncl
> Url:
> http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20100701/280001b9/attachment.pl
>
> ------------------------------
>
> Message: 3
> Date: Thu, 1 Jul 2010 17:24:17 -0600
> From: Mary Haley <haley@ucar.edu>
> Subject: Re: Adjusting tickmark labels with
> gsn_csm_contour_map
> To: Emilie Vanvyve <evanvyve@ucar.edu>
> Cc: ncl-talk@ucar.edu
> Message-ID: <BAE1F542-1704-4E05-B6C6-4748884E8C3F@ucar.edu>
> Content-Type: text/plain; charset=us-ascii
>
> Hi Emilie,
>
> It's not possible to control your own map tickmarks using the method you
> tried,
> unfortunately. This has been on our list of things to improve for awhile
> now.
>
> I created an example that shows how to create your own map tickmarks so
> you can put the labels where you want them.
>
> See example 9 at:
>
> http://www.ncl.ucar.edu/Applications/mptick.shtml#ex9
>
> You will want to copy the "add_map_tickmarks" function to your own script,
> and then call it as shown in this example. I used the same map resources
> that you set for the first frame of this example.
>
> You can control the size of the tickmark labels by setting:
>
> res@tmXBLabelFontHeightF
>
> This will affect the other axes labels as well.
>
> To control the length of the tickmarks, set:
>
> res@tmXBMajorOutwardLengthF
> res@tmXBOutwardLengthF
>
> to the same value. Try a small value, like 0.01 and then work from there.
> This will propagate to the other tickmark lengths as well.
>
>
> --Mary
>
>
> On Jul 1, 2010, at 9:01 AM, Emilie Vanvyve wrote:
>
> > Hello,
> >
> > I am new to NCL and struggling to customise a simple plot. I am not doing
> > anything fancy and don't understand why it doesn't work, hence this
> > e-mail.
> >
> > The plot is simple, it shows the topography of the US. I am trying to
> make
> > the lon/lat labels match the lon/lat grid lines overlaid on the filled
> > contours. I can modify the grid lines to be every 10 degrees, but I can't
> > modify the axis labels to match this (they are every 15 degrees in lon
> and
> > 5 degrees in lat). Note that I could as well modify the grid lines
> spacing
> > to match the labels, but this offers less flexibility.
> >
> > The topography is stored in the 'terrain' variable and is on a regular
> > longitude-latitude grid. I plot these data with:
> >
> > plot = gsn_csm_contour_map(wks,terrain,res)
> > draw(plot)
> > frame(wks)
> >
> > where 'res' has previously been set as:
> >
> > res = True
> > res@gsnDraw = False
> > res@gsnFrame = False
> >
> > res@mpProjection = "LambertConformal"
> > res@mpLambertParallel1F = 10
> > res@mpLambertParallel2F = 70
> > res@mpLambertMeridianF = -100
> > res@mpLimitMode = "LatLon"
> > res@mpMinLatF = 20
> > res@mpMaxLatF = 55
> > res@mpMinLonF = -140
> > res@mpMaxLonF = -60
> > res@mpGridAndLimbOn = True
> > res@mpGridSpacingF = 10.
> > res@pmTickMarkDisplayMode = "Always"
> >
> > res@tmXBMode = "Manual" ; I also tried "Explicit"
> > res@tmXBTickStartF = -130.
> > res@tmXBTickSpacingF = 10.
> > res@tmYLMode = "Manual"
> > res@tmYLTickStartF = 20.
> > res@tmYLTickSpacingF = 10.
> >
> > res@cnFillOn = True
> > res@cnLinesOn = False
> >
> > At execution of the script, I get:
> >
> > warning:MapPlotSetValues: tm[XB|XT|YL|YR]Mode resources are not
> > currently enabled for MapPlot tick marks
> >
> > I guess that's why it doesn't work, but why does this message appear and
> > how can I solve my problem?
> >
> > Sorry if this has already been asked, I did look through the on-line
> > documentation and previous e-mail messages, but couldn't find any
> > explanation (that I understand).
> >
> > Thanks in advance for any help,
> > Emilie
> >
> > _______________________________________________
> > ncl-talk mailing list
> > List instructions, subscriber options, unsubscribe:
> > http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
>
> ------------------------------
>
> Message: 4
> Date: Thu, 1 Jul 2010 23:52:33 +0000
> From: WuDraco <aos_rain@hotmail.com>
> Subject: overplot
> To: <ncl-talk@ucar.edu>
> Message-ID: <SNT134-w26234F8E88A6980BEC2EB0F9CD0@phx.gbl>
> Content-Type: text/plain; charset="gb2312"
>
>
> Hello,
>
> I would like to plot a contour. Then overlay a line box on it. However, the
> output file has the contour and the line box in two different figures. Do
> you have any idea what is wrong in my code? My code is enclosed. Thanks!
>
> -Longtao
>
> _________________________________________________________________
> ????????????????
>
> http://cn.bing.com/search?q=%E5%A4%A9%E6%B0%94%E9%A2%84%E6%8A%A5&form=MICHJ2
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL:
> http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20100701/f88579a8/attachment.html
> -------------- next part --------------
> A non-text attachment was scrubbed...
> Name: wrfout_domains.ncl
> Type: application/octet-stream
> Size: 2097 bytes
> Desc: not available
> Url :
> http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20100701/f88579a8/attachment.obj
>
> ------------------------------
>
> Message: 5
> Date: Thu, 1 Jul 2010 21:50:04 -0400 (EDT)
> From: debasish@sca.uqam.ca
> Subject: Re: overplot
> To: "WuDraco" <aos_rain@hotmail.com>
> Cc: ncl-talk@ucar.edu
> Message-ID:
> <50500.174.91.227.163.1278035404.squirrel@webmail.sca.uqam.ca>
> Content-Type: text/plain;charset=iso-8859-1
>
> Longtao,
>
> Please try it by adding
>
> mpres@gsnDraw = False ; don't draw yet
> mpres@gsnFrame = False ; don't advance ye
>
> Also check the example 4 from Polygons, Polymarkers, and Polylines webpage
> http://www.ncl.ucar.edu/Applications/polyg.shtml
>
> Thanks
> debasish
>
> >
> > Hello,
> >
> > I would like to plot a contour. Then overlay a line box on it. However,
> > the output file has the contour and the line box in two different
> figures.
> > Do you have any idea what is wrong in my code? My code is enclosed.
> > Thanks!
> >
> > -Longtao
> >
> > _________________________________________________________________
> > ????????????????????????????????
> >
> http://cn.bing.com/search?q=%E5%A4%A9%E6%B0%94%E9%A2%84%E6%8A%A5&form=MICHJ2_______________________________________________
> > ncl-talk mailing list
> > List instructions, subscriber options, unsubscribe:
> > http://mailman.ucar.edu/mailman/listinfo/ncl-talk
> >
>
>
>
>
> ------------------------------
>
> Message: 6
> Date: Mon, 7 Jun 2010 17:24:45 +0800
> From: ??? <annddyshow@gmail.com>
> Subject: how to calculate POP (Principal Oscillation
> Pattern) with NCL
> To: <ncl-talk@ucar.edu>
> Message-ID: <06BF5F11180A4B07B5B5CC6EC1C08002@XiaoCLPC>
> Content-Type: text/plain; charset="gb2312"
>
> Dear Madam/Sir,
> How to use NCL to do POP (Principal Oscillation Pattern) analysis? Is there
> any example?
>
> In the mailing list, I've got a similar thread
> http://www.ncl.ucar.edu/Support/talk_archives/2009/0265.html
>
> Thanks a lot!
> Xiao
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL:
> http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20100607/28fbb0f0/attachment.html
>
> ------------------------------
>
> Message: 7
> Date: Fri, 2 Jul 2010 06:44:26 +0000
> From: "Li, Yongzuo" <yongzuo.li@ou.edu>
> Subject: noise after 2D interpolation
> To: "ncl-talk@ucar.edu" <ncl-talk@ucar.edu>
> Cc: "karan@ngi.msstate.edu" <karan@ngi.msstate.edu>
> Message-ID:
> <
> 3D70D3147AC59C42A3C346D8F9FA2A0A2996505A@it-monad.sooner.net.ou.edu>
> Content-Type: text/plain; charset="us-ascii"
>
> Hi,
>
>
> I use
>
> ftp://ftp.erc.msstate.edu/outgoing/yli/ABL/wrf-minus-ruc.ncl
>
> to make
>
> ftp://ftp.erc.msstate.edu/outgoing/yli/ABL/wrf-minus-ruc-final.pdf
>
> There are a lot of noise on page 3 of the pdf file after
>
> u1d = rcm2points(lat2d,lon2d,uzoom,y1d,x1d,2)
> v1d = rcm2points(lat2d,lon2d,vzoom,y1d,x1d,2)
>
> Is any smooth function available in NCL to remove noise?
>
>
> Yongzuo
>
>
> ------------------------------
>
> Message: 8
> Date: Fri, 2 Jul 2010 15:54:05 +0530
> From: "Sabeerali(sebi)" <sabeerl@gmail.com>
> Subject: Make the spacing between the tick marks equal
> without changing the explicit tick mark values.
> To: ncl-talk@ucar.edu
> Message-ID:
> <AANLkTildedbjZvvqIueTysdvMCwnAc20Aj4KnITL1xE9@mail.gmail.com>
> Content-Type: text/plain; charset="iso-8859-1"
>
> Hi everyone,
>
> I have a frequency-wave number spectrum. In which the X axis is wave
> number(-4 to 4) and Y axis frequecy. The Y axis frequency values are
> explicitly defined.
> My Y axis values are (/0.01, 0.0125, 0.016666667, 0.025, 0.05, 0.1/). So
> when we go from bottom to top, obviously the spacing between the tick marks
> increases. That is quite clear. In the plot, I want to make the spacing
> between each tick mark to be equal without changing the above tick mark
> values( i.e. extent the bottom and compress the top of plot)
> Is there any resources or method in NCL so that I can make the spacing
> between the tick marks equal (By using the above values)
>
> Any help is highly appreciated..
> sabeer
>
>
> the plotting part is given below.
>
>
> wks = gsn_open_wks("eps", "freq_spec_new2")
> gsn_define_colormap(wks,"prcp_2")
> res=True
> res@cnFillOn=True
> res@gsnSpreadColors=True
> NW=4
> fMin=-0.5
> fMax=0.5
> WORK2=pow2({-NW:NW},{0:fMax})
> izero=ind(WORK2&freq.eq.0.0)
> WORK2(:,izero) = min(WORK2)
> res@tmYLMode = "Explicit" ; explicit labels
> pltPeriods = (/10,20,40,60,80,100/)
> res@tmYLValues = 1./pltPeriods ; freq location labels (/0.01,
> 0.0125, 0.016666667, 0.025, 0.05, 0.1/).
> res@tmYLLabels = pltPeriods ; labels themselves
> res@tmLabelAutoStride = True ; nice stride on labels
> res@tiXAxisOffsetYF = 0.2
>
> plot = gsn_csm_contour(wks,WORK2(freq|:,wave|:),res)
> end
>
>
> --
> **********************************
> Sabeerali
> Climate and Global Modeling Division
> Indian Institute of Tropical Meteorology
> Pashan, Pune, 411 008
> ****************************************
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL:
> http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20100702/6c1dc5d4/attachment.html
>
> ------------------------------
>
> Message: 9
> Date: Fri, 2 Jul 2010 07:36:43 -0600
> From: Mary Haley <haley@ucar.edu>
> Subject: Re: RAMS HDF5 format?
> To: "Dreher, Joseph Civ USAF AFTAC AFTAC/TMNV"
> <Joseph.Dreher@patrick.af.mil>
> Cc: ncl-talk Talk <ncl-talk@ucar.edu>
> Message-ID: <58231F56-9AD8-4949-B7E4-E0C2FEC94721@ucar.edu>
> Content-Type: text/plain; charset=windows-1252
>
>
> On Jul 1, 2010, at 10:33 AM, Dennis Shea wrote:
>
> > Hello,
> >
> > I believe that HDF5 will be in v6.0.0.
> >
> > There will be a 5.2.1 release before that.
> > I am not sure how much of HDF5 will be supported in that release.
>
> The 5.2.1 release will be mostly a bug fix release. We don't plan to have
> HDF5 support in it.
>
> However, in the near future, we might be able to give a 6.0.0 alpha version
> to users for testing
> HDF5.
>
> --Mary
>
> >
> > Can you please send a sample RAMS HDF5 file?
> > The NCL team can use it for testing. If NCL's development branch
> > reads the file successfully, I am pretty sure a beta version could
> > be created for your use.
> >
> > [1] Please send the output from:
> >
> > %> uname -a
> >
> > [2]
> >
> > ftp ftp.cgd.ucar.edu
> > anonymous
> > Use your email address for the password
> > cd incoming
> > put ...
> > put ...
> > :
> > quit
> >
> > Send me an email with the following information:
> >
> > * The names of the files after you have successfully completed the
> > transfer.
> > We cannot look at the ftp/incoming directory for security reasons.
> >
> >
> > Cheers
> >
> >
> > On 7/1/10 8:57 AM, Dreher, Joseph Civ USAF AFTAC AFTAC/TMNV wrote:
> >> Hello all,
> >>
> >> I recently attended the WRF user?s workshop and saw the overview on NCL.
> >> I believe I saw a note that the next release is supposed to include
> >> readers for HDF5 format? In my organization we currently are running
> >> RAMS and now WRF and would like to have the capability to read both the
> >> RAMS HDF5 and WRF netCDF file formats and compare with NCL. Is it still
> >> true that the next scheduled release will include HDF5?
> >>
> >> Thanks in advance,
> >>
> >> Joe
> >>
> >>
> >>
> >> _______________________________________________
> >> 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: 10
> Date: Fri, 2 Jul 2010 07:40:08 -0600
> From: Mary Haley <haley@ucar.edu>
> Subject: Re: Make the spacing between the tick marks equal
> without changing the explicit tick mark values.
> To: "Sabeerali(sebi)" <sabeerl@gmail.com>
> Cc: ncl-talk@ucar.edu
> Message-ID: <69727805-03ED-4C8A-978F-9BDE221AD477@ucar.edu>
> Content-Type: text/plain; charset=us-ascii
>
> Hi Sabeer,
>
> I think what you want is to set the gsnYAxisIrregular2Linear resource to
> True. See example 1 at:
>
> http://www.ncl.ucar.edu/Applications/axes.shtml#ex1
>
> --Mary
>
> On Jul 2, 2010, at 4:24 AM, Sabeerali(sebi) wrote:
>
> > Hi everyone,
> >
> > I have a frequency-wave number spectrum. In which the X axis is wave
> number(-4 to 4) and Y axis frequecy. The Y axis frequency values are
> explicitly defined.
> > My Y axis values are (/0.01, 0.0125, 0.016666667, 0.025, 0.05, 0.1/). So
> when we go from bottom to top, obviously the spacing between the tick marks
> increases. That is quite clear. In the plot, I want to make the spacing
> between each tick mark to be equal without changing the above tick mark
> values( i.e. extent the bottom and compress the top of plot)
> > Is there any resources or method in NCL so that I can make the spacing
> between the tick marks equal (By using the above values)
> >
> > Any help is highly appreciated..
> > sabeer
> >
> >
> > the plotting part is given below.
> >
> >
> > wks = gsn_open_wks("eps", "freq_spec_new2")
> > gsn_define_colormap(wks,"prcp_2")
> > res=True
> > res@cnFillOn=True
> > res@gsnSpreadColors=True
> > NW=4
> > fMin=-0.5
> > fMax=0.5
> > WORK2=pow2({-NW:NW},{0:fMax})
> > izero=ind(WORK2&freq.eq.0.0)
> > WORK2(:,izero) = min(WORK2)
> > res@tmYLMode = "Explicit" ; explicit labels
> > pltPeriods = (/10,20,40,60,80,100/)
> > res@tmYLValues = 1./pltPeriods ; freq location labels (/0.01,
> 0.0125, 0.016666667, 0.025, 0.05, 0.1/).
> > res@tmYLLabels = pltPeriods ; labels themselves
> > res@tmLabelAutoStride = True ; nice stride on
> labels
> > res@tiXAxisOffsetYF = 0.2
> >
> > plot = gsn_csm_contour(wks,WORK2(freq|:,wave|:),res)
> > end
> >
> >
> > --
> > **********************************
> > Sabeerali
> > Climate and Global Modeling Division
> > Indian Institute of Tropical Meteorology
> > Pashan, Pune, 411 008
> > ****************************************
> > _______________________________________________
> > ncl-talk mailing list
> > List instructions, subscriber options, unsubscribe:
> > http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
>
> ------------------------------
>
> Message: 11
> Date: Fri, 2 Jul 2010 07:43:25 -0600
> From: Mary Haley <haley@ucar.edu>
> Subject: Re: noise after 2D interpolation
> To: "Li, Yongzuo" <yongzuo.li@ou.edu>
> Cc: "ncl-talk@ucar.edu" <ncl-talk@ucar.edu>, "karan@ngi.msstate.edu"
> <karan@ngi.msstate.edu>
> Message-ID: <53D71B73-C15B-4DB7-9969-1A99C34E424D@ucar.edu>
> Content-Type: text/plain; charset=us-ascii
>
> Hi Yongzuo,
>
> What about setting the contour resource cnSmoothingOn to True? Or,
> if you are generating a raster plot (which doesn't seem to be the case
> here), you can set cnRasterSmoothingOn to True.
>
> --Mary
>
> On Jul 2, 2010, at 12:44 AM, Li, Yongzuo wrote:
>
> > Hi,
> >
> >
> > I use
> >
> > ftp://ftp.erc.msstate.edu/outgoing/yli/ABL/wrf-minus-ruc.ncl
> >
> > to make
> >
> > ftp://ftp.erc.msstate.edu/outgoing/yli/ABL/wrf-minus-ruc-final.pdf
> >
> > There are a lot of noise on page 3 of the pdf file after
> >
> > u1d = rcm2points(lat2d,lon2d,uzoom,y1d,x1d,2)
> > v1d = rcm2points(lat2d,lon2d,vzoom,y1d,x1d,2)
> >
> > Is any smooth function available in NCL to remove noise?
> >
> >
> > Yongzuo
> > _______________________________________________
> > ncl-talk mailing list
> > List instructions, subscriber options, unsubscribe:
> > http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
>
> ------------------------------
>
> Message: 12
> Date: Fri, 2 Jul 2010 15:04:34 +0000
> From: "Straub, Katherine" <straubk@susqu.edu>
> Subject: lag composites
> To: "ncl-talk@ucar.edu" <ncl-talk@ucar.edu>
> Message-ID:
> <92518EBDD479794BBF96B657409A42CE01A5A0@su-mail1.susqu.edu>
> Content-Type: text/plain; charset="us-ascii"
>
> I have a list of specific dates (YYYYMMDD) I'm using in a composite. I can
> make the composite for this list just fine, but am wondering if there is an
> easy way to extend the composite backward and forward in time, like from day
> -5 to day +5, while still compositing around actual dates rather than
> indices.
>
> Thanks,
> Kathy Straub
>
>
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL:
> http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20100702/97b5df9b/attachment.html
>
> ------------------------------
>
> Message: 13
> Date: Fri, 02 Jul 2010 10:09:11 -0600
> From: Dennis Shea <shea@ucar.edu>
> Subject: Re: noise after 2D interpolation
> To: "Li, Yongzuo" <yongzuo.li@ou.edu>
> Cc: "ncl-talk@ucar.edu" <ncl-talk@ucar.edu>, "karan@ngi.msstate.edu"
> <karan@ngi.msstate.edu>
> Message-ID: <4C2E0F27.6080700@ucar.edu>
> Content-Type: text/plain; charset=ISO-8859-1; format=flowed
>
>
> The v5.2.1 release has another option [opt=2]
>
> http://www.ncl.ucar.edu/Document/Functions/Built-in/rcm2points.shtml
>
> You can provide the following information
>
> %> gcc --version
>
> %> uname -a
>
> a v5.2.1 beta version can be made available.
>
>
> On 7/2/10 12:44 AM, Li, Yongzuo wrote:
> > Hi,
> >
> >
> > I use
> >
> > ftp://ftp.erc.msstate.edu/outgoing/yli/ABL/wrf-minus-ruc.ncl
> >
> > to make
> >
> > ftp://ftp.erc.msstate.edu/outgoing/yli/ABL/wrf-minus-ruc-final.pdf
> >
> > There are a lot of noise on page 3 of the pdf file after
> >
> > u1d = rcm2points(lat2d,lon2d,uzoom,y1d,x1d,2)
> > v1d = rcm2points(lat2d,lon2d,vzoom,y1d,x1d,2)
> >
> > Is any smooth function available in NCL to remove noise?
> >
> >
> > Yongzuo
> > _______________________________________________
> > ncl-talk mailing list
> > List instructions, subscriber options, unsubscribe:
> > http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
> ------------------------------
>
> Message: 14
> Date: Fri, 02 Jul 2010 10:32:59 -0600
> From: Dennis Shea <shea@ucar.edu>
> Subject: Re: lag composites
> To: "Straub, Katherine" <straubk@susqu.edu>
> Cc: "ncl-talk@ucar.edu" <ncl-talk@ucar.edu>
> Message-ID: <4C2E14BB.5030400@ucar.edu>
> Content-Type: text/plain; charset=windows-1252; format=flowed
>
> Not enough info.
>
>
>
> x = f->X ; (time,lat,lon)
> ymd = f->date ; (time) ... units=>YYYYMMDD
>
> say, 'time' is (say) 'hours/days since .."
>
>
> then, ideally, you could
>
> x&time = ymd
> printVarSummary(x) ; note the 'time' coordinate values
>
> This would allow u to access a specific YYYYMMDD via
> coordinate subscripting
>
> q = x({20030412},:,:)
>
> If the +/-5 were in the same month, then
>
> q = x({20030412-5:20030412+5},:,:)
>
> However, spanning months would not be possible
>
> Say, the oriinal units were 'days since ..'
>
> dayLag = 5
> centerDate = 20030412
> ntStrt = ind(ymd.eq.(centerDate-dayLag))
> ntLast = ind(ymd.eq.(centerDate+dayLag))
>
> Good luck
>
>
>
>
>
> On 7/2/10 9:04 AM, Straub, Katherine wrote:
> > I have a list of specific dates (YYYYMMDD) I?m using in a composite. I
> > can make the composite for this list just fine, but am wondering if
> > there is an easy way to extend the composite backward and forward in
> > time, like from day -5 to day +5, while still compositing around actual
> > dates rather than indices.
> >
> > Thanks,
> >
> > Kathy Straub
> >
> >
> >
> > _______________________________________________
> > ncl-talk mailing list
> > List instructions, subscriber options, unsubscribe:
> > http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
> ------------------------------
>
> Message: 15
> Date: Fri, 2 Jul 2010 12:39:01 -0400 (EDT)
> From: Ioana Colfescu <colfescu@cola.iges.org>
> Subject: invalid pointer error message
> To: ncl-talk <ncl-talk@ucar.edu>
> Message-ID:
> <143391364.11611.1278088741065.JavaMail.root@mail.iges.org>
> Content-Type: text/plain; charset=utf-8
>
> Hi,
>
> Could someone please help me with the following -I use the following code :
>
> FILES=systemfunc("ls /data/oldfarm2/colfi/SODA/*.cdf")
> f1=addfiles(FILES,"r")
> ListSetType(f1,"cat");
> t=f1[:]->temp
>
> which works fine for other data but for this I get this error :
>
> *** glibc detected *** ncl: munmap_chunk(): invalid pointer:
> 0x00002b1bbc0e3010 ***
> ======= Backtrace: =========
> /lib64/libc.so.6(cfree+0x166)[0x37bca72886]
> ncl(NclFree+0x25)[0x5681de]
> ncl[0x5495bb]
> ncl(_NclDestroyObj+0x56)[0x553e09]
> ncl(CallLIST_READ_FILEVAR_OP+0x2900)[0x5e0398]
> ncl(_NclExecute+0x442)[0x5e09c5]
> ncl(yyparse+0x792)[0x51fc69]
> ncl(main+0xa3e)[0x51b4ce]
> /lib64/libc.so.6(__libc_start_main+0xf4)[0x37bca1d994]
> ncl(__gxx_personality_v0+0x14a)[0x51a9fa]
> ======= Memory map: ========
> 00400000-01437000 r-xp 00000000 00:17 2147824247
> /homes/colfescu/ncl/bin/nc
> 01537000-01722000 rw-p 01037000 00:17 2147824247
> /homes/colfescu/ncl/bin/ncl
>
>
> I don't know why the error appears.
>
> Ioana
>
>
> ------------------------------
>
> Message: 16
> Date: Fri, 2 Jul 2010 17:17:29 +0000
> From: "Li, Yongzuo" <yongzuo.li@ou.edu>
> Subject: v5.2.1 for rcm2points to solve interpolation noise
> To: "ncl-talk@ucar.edu" <ncl-talk@ucar.edu>
> Cc: "karan@ngi.msstate.edu" <karan@ngi.msstate.edu>
> Message-ID:
> <
> 3D70D3147AC59C42A3C346D8F9FA2A0A299650AF@it-monad.sooner.net.ou.edu>
> Content-Type: text/plain; charset="us-ascii"
>
> Hi ,
>
> According to Dennis's e-mail below, I provide
>
> Bazooka[33] yli$ gcc --version
> gcc (GCC) 4.1.2 20070115 (SUSE Linux)
> Copyright (C) 2006 Free Software Foundation, Inc.
> This is free software; see the source for copying conditions. There is NO
> warranty; not even for MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
>
> Bazooka[34] yli$ uname -a
> Linux Bazooka.HPC.MsState.Edu 2.6.16.60-0.42.8-smp #1 SMP Tue Dec 15
> 17:28:00 UTC 2009 x86_64 x86_64 x86_64 GNU/Linux
> Bazooka[35] yli$
>
> What is next step?
>
> Yongzuo
>
> ________________________________________
> From: Dennis Shea [shea@ucar.edu]
> Sent: Friday, July 02, 2010 11:09 AM
> To: Li, Yongzuo
> Cc: ncl-talk@ucar.edu; karan@ngi.msstate.edu
> Subject: Re: noise after 2D interpolation
>
> The v5.2.1 release has another option [opt=2]
>
> http://www.ncl.ucar.edu/Document/Functions/Built-in/rcm2points.shtml
>
> You can provide the following information
>
> %> gcc --version
>
> %> uname -a
>
> a v5.2.1 beta version can be made available.
>
>
> On 7/2/10 12:44 AM, Li, Yongzuo wrote:
> > Hi,
> >
> >
> > I use
> >
> > ftp://ftp.erc.msstate.edu/outgoing/yli/ABL/wrf-minus-ruc.ncl
> >
> > to make
> >
> > ftp://ftp.erc.msstate.edu/outgoing/yli/ABL/wrf-minus-ruc-final.pdf
> >
> > There are a lot of noise on page 3 of the pdf file after
> >
> > u1d = rcm2points(lat2d,lon2d,uzoom,y1d,x1d,2)
> > v1d = rcm2points(lat2d,lon2d,vzoom,y1d,x1d,2)
> >
> > Is any smooth function available in NCL to remove noise?
> >
> >
> > Yongzuo
> > _______________________________________________
> > ncl-talk mailing list
> > List instructions, subscriber options, unsubscribe:
> > http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
> ------------------------------
>
> Message: 17
> Date: Fri, 02 Jul 2010 11:30:05 -0600
> From: Dennis Shea <shea@ucar.edu>
> Subject: Re: invalid pointer error message
> To: Ioana Colfescu <colfescu@cola.iges.org>
> Cc: ncl-talk <ncl-talk@ucar.edu>
> Message-ID: <4C2E221D.5050409@ucar.edu>
> Content-Type: text/plain; charset=ISO-8859-1; format=flowed
>
> Not enough information.
>
> I speculate that you have run out of available memory or, perhaps,
> the collective size of 't' exceeded 2 GB.
>
> Good luck
>
> On 7/2/10 10:39 AM, Ioana Colfescu wrote:
> > Hi,
> >
> > Could someone please help me with the following -I use the following code
> :
> >
> > FILES=systemfunc("ls /data/oldfarm2/colfi/SODA/*.cdf")
> > f1=addfiles(FILES,"r")
> > ListSetType(f1,"cat");
> > t=f1[:]->temp
> >
> > which works fine for other data but for this I get this error :
> >
> > *** glibc detected *** ncl: munmap_chunk(): invalid pointer:
> 0x00002b1bbc0e3010 ***
> > ======= Backtrace: =========
> > /lib64/libc.so.6(cfree+0x166)[0x37bca72886]
> > ncl(NclFree+0x25)[0x5681de]
> > ncl[0x5495bb]
> > ncl(_NclDestroyObj+0x56)[0x553e09]
> > ncl(CallLIST_READ_FILEVAR_OP+0x2900)[0x5e0398]
> > ncl(_NclExecute+0x442)[0x5e09c5]
> > ncl(yyparse+0x792)[0x51fc69]
> > ncl(main+0xa3e)[0x51b4ce]
> > /lib64/libc.so.6(__libc_start_main+0xf4)[0x37bca1d994]
> > ncl(__gxx_personality_v0+0x14a)[0x51a9fa]
> > ======= Memory map: ========
> > 00400000-01437000 r-xp 00000000 00:17 2147824247
> /homes/colfescu/ncl/bin/nc
> > 01537000-01722000 rw-p 01037000 00:17 2147824247
> /homes/colfescu/ncl/bin/ncl
> >
> >
> > I don't know why the error appears.
> >
> > Ioana
> > _______________________________________________
> > ncl-talk mailing list
> > List instructions, subscriber options, unsubscribe:
> > http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
> ------------------------------
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk@ucar.edu
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
> End of ncl-talk Digest, Vol 80, Issue 2
> ***************************************
>

-- 
**********************************
Sabeerali
Climate and Global Modeling Division
Indian Institute of Tropical Meteorology
Pashan, Pune, 411 008
****************************************

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Sun Jul 4 22:59:15 2010