Sabeer,
If you want to take a linear axis and make it irregular, you need to
provide a new coordinate array that defines your new Y axis values.
NCL can't make any assumptions about recalibrating the Y axis based on
tick mark values.
You currently have a "freq" coordinate array that is linearly spaced.
One thing you can do is provide a new coordinate array with 15
elements that defines the new spacing that you want to see on the Y
axis. You have a start with the values you provided so far; can you
fill this in to get 15 values?
--Mary
On Jul 4, 2010, at 10:59 PM, Sabeerali(sebi) wrote:
> Thanks Mary
> Yes, I want the plot similar to example1 (axes_1.ncl)..But when I
> use res@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
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Jul 6 10:47:18 2010
This archive was generated by hypermail 2.1.8 : Fri Jul 09 2010 - 16:06:22 MDT