Re: ncl-talk Digest, Vol 114, Issue 2

From: jayasankar C B <cbjayasankar_at_nyahnyahspammersnyahnyah>
Date: Fri May 03 2013 - 04:49:56 MDT

dear From: Sha Feng

Subject: Re:help for overly xy-plot and contour-plot using NCL
can you try this link
i think u dont have to use overlay..
u can do by using

gsn_csm_xy2 command.

http://www.ncl.ucar.edu/Applications/Scripts/xy_7.ncl

On Thu, May 2, 2013 at 11:30 PM, <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. help for overly xy-plot and contour-plot using NCL (Sha Feng)
> 2. Plotting skewt from RUC grib files?
> (Gerry Creager - NOAA Affiliate)
> 3. Re: help for overly xy-plot and contour-plot using NCL
> (Will Hobbs)
> 4. Contours overlay problem (Soumik Basu)
> 5. Re: Contours overlay problem (Will Hobbs)
> 6. How to set the distance among legend so that they do not
> overlap (juki juki)
> 7. Re: Correlation Question (Nkese Mc Shine)
> 8. Re: Correlation Question (Dennis Shea)
> 9. Re: How to set the distance among legend so that they do not
> overlap (Dennis Shea)
> 10. Re: Correlation Question (Nkese Mc Shine)
> 11. Re: How to set the distance among legend so that they do not
> overlap (Mary Haley)
> 12. Re: Plotting skewt from RUC grib files? (Dennis Shea)
> 13. Re: Plotting skewt from RUC grib files?
> (Gerry Creager - NOAA Affiliate)
> 14. How to fix plot size and shape independent of paper size
> (Douglas Lowe)
> 15. Re: Correlation Question (Dennis Shea)
> 16. How to control "correlation" text position in Taylor Diagram.
> (mmkamal@uwaterloo.ca)
>
>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Wed, 1 May 2013 12:16:06 -0700
> From: Sha Feng <sfeng@mail.usf.edu>
> Subject: help for overly xy-plot and contour-plot using NCL
> To: "ncl-talk@ucar.edu" <ncl-talk@ucar.edu>
> Message-ID:
> <CAO_gYL_6ZaM5N_1P23+QX+=
> iFVwVk624zUaHPWW5oPb2KWM2xA@mail.gmail.com>
> Content-Type: text/plain; charset="iso-8859-1"
>
> Hi ncl-talkers,
>
> I was trying to plot a time-height distribution of vertical velocity and
> overlay the precipitation on it in the same figure. I tried "overlay" and
> "polylines", but all failed.
>
> Data are omega(time,z) and rain(time). They both have the same time
> dimension.
>
> In the figure, the x-axis is time dimension, the left y-axis is height, and
> the right y-axis is precipitation amount.
>
> Any clues to make it happen?
>
> Thank you so much,
> Sha
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL:
> http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20130501/2db5d216/attachment.html
>
> ------------------------------
>
> Message: 2
> Date: Wed, 1 May 2013 17:19:05 -0500
> From: Gerry Creager - NOAA Affiliate <gerry.creager@noaa.gov>
> Subject: Plotting skewt from RUC grib files?
> To: ncl-talk <ncl-talk@ucar.edu>
> Message-ID:
> <
> CAFufYD7CzEi1KFfne2xWSra3u13P74VT-heySTM2cmUY_s1miw@mail.gmail.com>
> Content-Type: text/plain; charset="windows-1252"
>
> I'm trying to sort through this exercise. Has anyone else done this already
> before I reinvent the wheel?
>
> Thanks, Gerry
> --
> 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)
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL:
> http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20130501/3d7cbb04/attachment.html
>
> ------------------------------
>
> Message: 3
> Date: Wed, 01 May 2013 23:13:03 +0000
> From: Will Hobbs <Will.Hobbs@utas.edu.au>
> Subject: Re: help for overly xy-plot and contour-plot using
> NCL
> To: Sha Feng <sfeng@mail.usf.edu>, "ncl-talk@ucar.edu"
> <ncl-talk@ucar.edu>
> Message-ID:
> <7B230BB1F8E9BD449F56BA16FC664F793FB3E2BF@MBXSBYN5.utas.ad.internal
> >
> Content-Type: text/plain; charset="us-ascii"
>
> Sha
>
> Polylines are not the way to go for what you want to do.
>
> There are a number of ways of plotting what you want, all of which are
> shown in the 'contour-on-contour' example page; if you haven't looked at
> this then I strongly urge you to.
>
> http://www.ncl.ucar.edu/Applications/conOncon.shtml
>
> The overlay method should work, but without seeing your script it's
> impossible to see what the problem is. If you haven't already done so, take
> a look at the example page, if you're still having problems then attach or
> cut and paste a relatively clean version of your script (commenting is
> always appreciated).
>
> Will
>
> From: Sha Feng <sfeng@mail.usf.edu<mailto:sfeng@mail.usf.edu>>
> Date: Thursday, 2 May 2013 5:16 AM
> To: "ncl-talk@ucar.edu<mailto:ncl-talk@ucar.edu>" <ncl-talk@ucar.edu
> <mailto:ncl-talk@ucar.edu>>
> Subject: help for overly xy-plot and contour-plot using NCL
>
> Hi ncl-talkers,
>
> I was trying to plot a time-height distribution of vertical velocity and
> overlay the precipitation on it in the same figure. I tried "overlay" and
> "polylines", but all failed.
>
> Data are omega(time,z) and rain(time). They both have the same time
> dimension.
>
> In the figure, the x-axis is time dimension, the left y-axis is height,
> and the right y-axis is precipitation amount.
>
> Any clues to make it happen?
>
> Thank you so much,
> Sha
>
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL:
> http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20130501/08dc070e/attachment.html
>
> ------------------------------
>
> Message: 4
> Date: Wed, 01 May 2013 17:48:46 -0800
> From: Soumik Basu <sbasu@iarc.uaf.edu>
> Subject: Contours overlay problem
> To: Ncl Talk <ncl-talk@ucar.edu>
> Message-ID: <faef538dac4f01f9ed7240dc38e1e495@arthur.iarc.uaf.edu>
> Content-Type: text/plain; charset=UTF-8; format=flowed
>
> Hi,
>
> I am trying to overlay non shaded contours on shaded contours (this part
> is working fine) but I want my non shaded contour lines to be colored
> following the same color table as the shaded contours. Possibly I am
> missing some resources as a result of which my non shaded contours are
> black.
>
> Is there any example showing that? here is the part of the script
> describing the resources for shaded and non shaded contours.
>
> ;************************************************************************
> ; Resources for shaded
> ;************************************************************************
> res = True
>
> res@gsnDraw = False
> res@gsnFrame = False
> res@gsnMaximize = False
> res@gsnSpreadColors = True
> res@gsnSpreadColors = True ; use full
> colormap
> res@gsnAddCyclic = True
>
> res@cnFillOn = True
> res@cnLinesOn = False
> res@cnLevelSelectionMode = "ManualLevels"
> res@cnMinLevelValF = -1000
> res@cnMaxLevelValF = 1000
> res@cnLevelSpacingF = 50
> res@cnInfoLabelOn = False ; turn off
> cn info label
> res@cnLineLabelsOn = False
> res@cnLineLabelFontHeightF = 0.015
>
> res@tmLabelAutoStride = True
> res@tmXTOn = False
> res@tmYROn = False
> res@tmXBLabelFontHeightF = 0.030
> res@tmYLLabelFontHeightF = 0.030
>
> res@mpLimitMode = "LatLon"
> res@mpMinLatF = 40.
> res@mpMaxLatF = 60.
> res@mpMinLonF = 290.
> res@mpMaxLonF = 380.
> res@mpGeophysicalLineThicknessF = 1.75
> res@mpCenterLonF = 0
>
> res@lbOrientation = "vertical" ;vertical label bar
> res@lbLabelBarOn = False
> res@lbLabelAutoStride = True
>
> ;************************************************************************
> ; Resources for non shaded
> ;************************************************************************
>
> sres = True
>
> sres@gsnDraw = False
> sres@gsnFrame = False
> sres@gsnMaximize = False
> sres@gsnSpreadColors = True
> sres@gsnSpreadColors = True ; use full
> colormap
> sres@gsnAddCyclic = True
>
> sres@cnFillOn = False
> sres@cnLinesOn = True
> sres@cnLevelSelectionMode = "ManualLevels"
> sres@cnMinLevelValF = -1000
> sres@cnMaxLevelValF = 1000
> sres@cnLevelSpacingF = 50
> sres@cnInfoLabelOn = False ; turn off
> cn info label
> sres@cnLineLabelsOn = False
> sres@cnLineLabelFontHeightF = 0.015
>
> sres@tmLabelAutoStride = True
> sres@tmXTOn = False
> sres@tmYROn = False
> sres@tmXBLabelFontHeightF = 0.030
> sres@tmYLLabelFontHeightF = 0.030
>
> sres@lbOrientation = "vertical" ;vertical label bar
> sres@lbLabelBarOn = False
> sres@lbLabelAutoStride = True
>
>
> Thanks,
> Soumik
>
> --
> "Numbing the pain for a while will make it worse when you finally feel
> it.
> ~J.K. Rowling"
> *****************************************
> Soumik Basu
> Graduate Student, Research Assistant
> International Arctic Research Center
> University of Alaska Fairbanks
> Fairbanks,Alaska,USA
> *****************************************
>
>
> ------------------------------
>
> Message: 5
> Date: Thu, 02 May 2013 02:01:17 +0000
> From: Will Hobbs <Will.Hobbs@utas.edu.au>
> Subject: Re: Contours overlay problem
> To: Soumik Basu <sbasu@iarc.uaf.edu>, Ncl Talk <ncl-talk@ucar.edu>
> Message-ID:
> <7B230BB1F8E9BD449F56BA16FC664F793FB3F30E@MBXSBYN5.utas.ad.internal
> >
> Content-Type: text/plain; CHARSET=US-ASCII
>
> Soumik
>
>
> You should be able to colour your contour lines by setting
>
> sres@cnMonoLineColors = False
>
> getvalues plot@contour ;retrieve contour plot colormap
> "cnFiLLColors" : cols
> end getvalues
>
> sres@cnLineColors = cols ;set line colors
>
> The 'getvalues' statement is outlined in a previous post:
> http://www.ncl.ucar.edu/Support/talk_archives/2012/0731.html
>
> Obviously, you will have to have already plotted a filled contour base
> plot(named 'plot' in my example) or there won't be a colormap to retrieve.
> You can then draw the contour lines on a separate plot, and use overlay()
> to add it to the base plot.
>
> Will
>
> On 2/05/13 11:48 AM, "Soumik Basu" <sbasu@iarc.uaf.edu> wrote:
>
> >Hi,
> >
> >I am trying to overlay non shaded contours on shaded contours (this part
> >is working fine) but I want my non shaded contour lines to be colored
> >following the same color table as the shaded contours. Possibly I am
> >missing some resources as a result of which my non shaded contours are
> >black.
> >
> >Is there any example showing that? here is the part of the script
> >describing the resources for shaded and non shaded contours.
> >
> >;************************************************************************
> >; Resources for shaded
> >;************************************************************************
> >res = True
> >
> >res@gsnDraw = False
> >res@gsnFrame = False
> >res@gsnMaximize = False
> >res@gsnSpreadColors = True
> >res@gsnSpreadColors = True ; use full
> >colormap
> >res@gsnAddCyclic = True
> >
> >res@cnFillOn = True
> >res@cnLinesOn = False
> >res@cnLevelSelectionMode = "ManualLevels"
> >res@cnMinLevelValF = -1000
> >res@cnMaxLevelValF = 1000
> >res@cnLevelSpacingF = 50
> >res@cnInfoLabelOn = False ; turn off
> >cn info label
> >res@cnLineLabelsOn = False
> >res@cnLineLabelFontHeightF = 0.015
> >
> >res@tmLabelAutoStride = True
> >res@tmXTOn = False
> >res@tmYROn = False
> >res@tmXBLabelFontHeightF = 0.030
> >res@tmYLLabelFontHeightF = 0.030
> >
> >res@mpLimitMode = "LatLon"
> >res@mpMinLatF = 40.
> >res@mpMaxLatF = 60.
> >res@mpMinLonF = 290.
> >res@mpMaxLonF = 380.
> >res@mpGeophysicalLineThicknessF = 1.75
> >res@mpCenterLonF = 0
> >
> >res@lbOrientation = "vertical" ;vertical label bar
> >res@lbLabelBarOn = False
> >res@lbLabelAutoStride = True
> >
> >;************************************************************************
> >; Resources for non shaded
> >;************************************************************************
> >
> >sres = True
> >
> >sres@gsnDraw = False
> >sres@gsnFrame = False
> >sres@gsnMaximize = False
> >sres@gsnSpreadColors = True
> >sres@gsnSpreadColors = True ; use full
> >colormap
> >sres@gsnAddCyclic = True
> >
> >sres@cnFillOn = False
> >sres@cnLinesOn = True
> >sres@cnLevelSelectionMode = "ManualLevels"
> >sres@cnMinLevelValF = -1000
> >sres@cnMaxLevelValF = 1000
> >sres@cnLevelSpacingF = 50
> >sres@cnInfoLabelOn = False ; turn off
> >cn info label
> >sres@cnLineLabelsOn = False
> >sres@cnLineLabelFontHeightF = 0.015
> >
> >sres@tmLabelAutoStride = True
> >sres@tmXTOn = False
> >sres@tmYROn = False
> >sres@tmXBLabelFontHeightF = 0.030
> >sres@tmYLLabelFontHeightF = 0.030
> >
> >sres@lbOrientation = "vertical" ;vertical label bar
> >sres@lbLabelBarOn = False
> >sres@lbLabelAutoStride = True
> >
> >
> >Thanks,
> >Soumik
> >
> >--
> >"Numbing the pain for a while will make it worse when you finally feel
> >it.
> >~J.K. Rowling"
> >*****************************************
> >Soumik Basu
> >Graduate Student, Research Assistant
> >International Arctic Research Center
> >University of Alaska Fairbanks
> >Fairbanks,Alaska,USA
> >*****************************************
> >_______________________________________________
> >ncl-talk mailing list
> >List instructions, subscriber options, unsubscribe:
> >http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
>
> ------------------------------
>
> Message: 6
> Date: Wed, 1 May 2013 21:41:02 -0700 (PDT)
> From: juki juki <juky_emc2@yahoo.com>
> Subject: How to set the distance among legend so that they
> do not overlap
> To: ncl-talk <ncl-talk@ucar.edu>
> Message-ID:
> <1367469662.88461.YahooMailNeo@web162801.mail.bf1.yahoo.com>
> Content-Type: text/plain; charset="iso-8859-1"
>
> Hi all ncl users..
>
> Thanks for your help and nice discussion so far. I have plotted my data to
> XY plot with 16 legends. The problem is that the legends overlap each
> other. How to overcome ?such problem ? The following is simple ncl code to
> control my current legend:
>
> ?res@pmLegendDisplayMode ? ?= "Always" ? ? ? ? ? ?; turn on legend
> ? res@pmLegendSide ? ? ? ? ? = "Top" ? ? ? ? ? ? ? ; Change location of?
> ? res@pmLegendParallelPosF ? = .6 ? ? ? ? ? ? ? ? ?; move units right
> ? res@pmLegendOrthogonalPosF = -0.25 ? ? ? ? ? ? ? ?; move units down
> ? res@pmLegendWidthF ? ? ? ? = 0.06 ? ? ? ? ? ? ? ?; Change width and
> ? res@pmLegendHeightF ? ? ? ?= 0.08 ? ? ? ? ? ? ? ?; height of legend.
> ? res@lgPerimOn ? ? ? ? ? ? ?= False ? ? ? ? ? ? ? ; turn off box around
> ? res@lgLabelFontHeightF ? ? = .012 ? ? ? ? ? ? ? ? ; label font height
>
> Thank you.
> joeky
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL:
> http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20130501/480bdd5e/attachment.html
>
> ------------------------------
>
> Message: 7
> Date: Thu, 2 May 2013 09:22:52 +0000
> From: Nkese Mc Shine <Nkese.McShine@sta.uwi.edu>
> Subject: Re: Correlation Question
> To: Will Hobbs <Will.Hobbs@utas.edu.au>, "ncl-talk@ucar.edu"
> <ncl-talk@ucar.edu>
> Message-ID:
> <7368311F6A984441893E050C8422DDBA3126D68E@exMS01vm.sauwi.uwi.tt>
> Content-Type: text/plain; charset="iso-8859-1"
>
> Will,
>
> Thank you for the explanation and suggestion however I do not think that
> will work as pattern_cor the rightmost dimension must be lon, and my
> monthly data which is 1 value is in time.
>
> Thank you.
> Nkese.
> ________________________________
> From: Will Hobbs [Will.Hobbs@utas.edu.au]
> Sent: 30 April 2013 19:37
> To: Nkese Mc Shine; ncl-talk@ucar.edu
> Subject: Re: Correlation Question
>
> Nkese
>
> >I would like to make a month to month correlation. I would like to
> correlate 1 month of sst(lat,lon,time) with 1 month of precipitation which
> actually is 1 value. Is it possible?
>
> No, you can't correlate a single value (in this case a month). This is
> nothing to do with NCL's capability, but fundamental to what a correlation
> is, i.e. a comparison of two variances. A single realisation (in your case
> a month) does not have a variance.
>
> What you can do (in case I misunderstood your question) is correlate the
> spatial pattern with that of a a single month; the output would be a time
> series of correlations. This can be done in NCL using the function
> 'pattern_cor()'
>
> Will
>
> From: Nkese Mc Shine <Nkese.McShine@sta.uwi.edu<mailto:
> Nkese.McShine@sta.uwi.edu>>
> Date: Wednesday, 1 May 2013 12:11 AM
> To: "ncl-talk@ucar.edu<mailto:ncl-talk@ucar.edu>" <ncl-talk@ucar.edu
> <mailto:ncl-talk@ucar.edu>>
> Subject: Correlation Question
>
> Dear All,
>
> I would like to make a month to month correlation. I would like to
> correlate 1 month of sst(lat,lon,time) with 1 month of precipitation which
> actually is 1 value. Is it possible? Attached is my script and below is the
> output. Please let me know, what I can do if it is possible.
>
> Variable: sst_anom
> Type: float
> Total Size: 105408000 bytes
> 26352000 values
> Number of Dimensions: 3
> Dimensions and sizes: [lat | 100] x [lon | 360] x [time | 732]
> Coordinates:
> lat: [49.5..-49.5]
> lon: [-179.5..179.5]
> time: [194601..200612]
> Number Of Attributes: 12
> anomaly_op_ncl : Annual Cycle Removed:function
> rmMonAnnCycLLT:contributed.ncl
> description : HadISST 1.1 monthly average sea surface temperature
> actual_range : ( -1.8, 34.76296 )
> missing_value : -1e+30
> _FillValue : -1e+30
> scale_factor : 1
> cell_methods : time: lat: lon: mean
> add_offset : 0
> units : degC
> standard_name : sea_surface_temperature
> long_name : Monthly 1 degree resolution SST
> reference : function rmMonAnnCycLLT in contrib.ncl
> (0)
> (0) Monthly 1 degree resolution SST: min=-8.70077 max=6.98223
>
> Variable: sst_anom_jan
> Type: float
> Total Size: 144000 bytes
> 36000 values
> Number of Dimensions: 3
> Dimensions and sizes: [lat | 100] x [lon | 360] x [time | 1]
> Coordinates:
> lat: [49.5..-49.5]
> lon: [-179.5..179.5]
> time: [199301..199301]
> Number Of Attributes: 12
> reference : function rmMonAnnCycLLT in contrib.ncl
> long_name : Monthly 1 degree resolution SST
> standard_name : sea_surface_temperature
> units : degC
> add_offset : 0
> cell_methods : time: lat: lon: mean
> scale_factor : 1
> _FillValue : -1e+30
> missing_value : -1e+30
> actual_range : ( -1.8, 34.76296 )
> description : HadISST 1.1 monthly average sea surface temperature
> anomaly_op_ncl : Annual Cycle Removed:function
> rmMonAnnCycLLT:contributed.ncl
> (0)
> (0) Monthly 1 degree resolution SST: min=-2.4631 max=2.59599
>
> Variable: precipanom_data_ts
> Type: float
> Total Size: 2928 bytes
> 732 values
> Number of Dimensions: 1
> Dimensions and sizes: [time | 732]
> Coordinates:
> time: [194601..200612]
> Number Of Attributes: 1
> _FillValue : 9.96921e+36
> (0)
> (0) min=-1.73324 max=3.01711
>
> Variable: precipJJA1993
> Type: float
> Total Size: 4 bytes
> 1 values
> Number of Dimensions: 1
> Dimensions and sizes: [time | 1]
> Coordinates:
> time: [199307..199307]
> Number Of Attributes: 2
> _FillValue : 9.96921e+36
> NMO : 6
> (0)
> (0) min=0.267831 max=0.267831
> warning:escorc: Non-fatal conditions encountered in series or xstd equals
> zero.
> Possibly, all values of a series are constant.
>
> warning:escorc: Most likely, one or more series consisted of all constant
> values
>
> Variable: corr
> Type: float
> Total Size: 144000 bytes
> 36000 values
> Number of Dimensions: 2
> Dimensions and sizes: [100] x [360]
> Coordinates:
> Number Of Attributes: 1
> _FillValue : -1e+30
> (0)
> (0) min=-1e+30 max=-1e+30
>
> Regards,
> Nkese.
>
> CONFIDENTIALITY: This email (including any attachments) may contain
> confidential, proprietary and/or privileged information. Any duplication,
> copying, distribution, dissemination, transmission, disclosure or use in
> any manner of this email (including any attachments) without the
> authorisation of the sender is strictly prohibited. If you receive this
> email (including any attachments) in error, please notify the sender and
> delete this email (including any attachments) from your system. Thank you.
>
>
> CONFIDENTIALITY: This email (including any attachments) may contain
> confidential, proprietary and/or privileged information. Any duplication,
> copying, distribution, dissemination, transmission, disclosure or use in
> any manner of this email (including any attachments) without the
> authorisation of the sender is strictly prohibited. If you receive this
> email (including any attachments) in error, please notify the sender and
> delete this email (including any attachments) from your system. Thank you.
>
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL:
> http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20130502/089eaaa2/attachment.html
>
> ------------------------------
>
> Message: 8
> Date: Thu, 02 May 2013 06:30:50 -0600
> From: Dennis Shea <shea@ucar.edu>
> Subject: Re: Correlation Question
> To: Nkese Mc Shine <Nkese.McShine@sta.uwi.edu>
> Cc: "ncl-talk@ucar.edu" <ncl-talk@ucar.edu>
> Message-ID: <51825C7A.9070805@ucar.edu>
> Content-Type: text/plain; charset=ISO-8859-1; format=flowed
>
> Ho Nkese
>
> If you have an array: x(lat,lon,time)
>
> You can use the named dimensions to reorder
>
> xNew = x(time|:,lat|:,lon|:)
> printVarSummary(xNew)
>
> Then you can compute the pattern correaltion(s)
>
> Good Luck
>
> On 5/2/13 3:22 AM, Nkese Mc Shine wrote:
> > Will,
> >
> > Thank you for the explanation and suggestion however I do not think that
> will work as pattern_cor the rightmost dimension must be lon, and my
> monthly data which is 1 value is in time.
> >
> > Thank you.
> > Nkese.
> > ________________________________
> > From: Will Hobbs [Will.Hobbs@utas.edu.au]
> > Sent: 30 April 2013 19:37
> > To: Nkese Mc Shine; ncl-talk@ucar.edu
> > Subject: Re: Correlation Question
> >
> > Nkese
> >
> >> I would like to make a month to month correlation. I would like to
> correlate 1 month of sst(lat,lon,time) with 1 month of precipitation which
> actually is 1 value. Is it possible?
> >
> > No, you can't correlate a single value (in this case a month). This is
> nothing to do with NCL's capability, but fundamental to what a correlation
> is, i.e. a comparison of two variances. A single realisation (in your case
> a month) does not have a variance.
> >
> > What you can do (in case I misunderstood your question) is correlate the
> spatial pattern with that of a a single month; the output would be a time
> series of correlations. This can be done in NCL using the function
> 'pattern_cor()'
> >
> > Will
> >
> > From: Nkese Mc Shine <Nkese.McShine@sta.uwi.edu<mailto:
> Nkese.McShine@sta.uwi.edu>>
> > Date: Wednesday, 1 May 2013 12:11 AM
> > To: "ncl-talk@ucar.edu<mailto:ncl-talk@ucar.edu>" <ncl-talk@ucar.edu
> <mailto:ncl-talk@ucar.edu>>
> > Subject: Correlation Question
> >
> > Dear All,
> >
> > I would like to make a month to month correlation. I would like to
> correlate 1 month of sst(lat,lon,time) with 1 month of precipitation which
> actually is 1 value. Is it possible? Attached is my script and below is the
> output. Please let me know, what I can do if it is possible.
> >
> > Variable: sst_anom
> > Type: float
> > Total Size: 105408000 bytes
> > 26352000 values
> > Number of Dimensions: 3
> > Dimensions and sizes: [lat | 100] x [lon | 360] x [time | 732]
> > Coordinates:
> > lat: [49.5..-49.5]
> > lon: [-179.5..179.5]
> > time: [194601..200612]
> > Number Of Attributes: 12
> > anomaly_op_ncl : Annual Cycle Removed:function
> rmMonAnnCycLLT:contributed.ncl
> > description : HadISST 1.1 monthly average sea surface temperature
> > actual_range : ( -1.8, 34.76296 )
> > missing_value : -1e+30
> > _FillValue : -1e+30
> > scale_factor : 1
> > cell_methods : time: lat: lon: mean
> > add_offset : 0
> > units : degC
> > standard_name : sea_surface_temperature
> > long_name : Monthly 1 degree resolution SST
> > reference : function rmMonAnnCycLLT in contrib.ncl
> > (0)
> > (0) Monthly 1 degree resolution SST: min=-8.70077 max=6.98223
> >
> > Variable: sst_anom_jan
> > Type: float
> > Total Size: 144000 bytes
> > 36000 values
> > Number of Dimensions: 3
> > Dimensions and sizes: [lat | 100] x [lon | 360] x [time | 1]
> > Coordinates:
> > lat: [49.5..-49.5]
> > lon: [-179.5..179.5]
> > time: [199301..199301]
> > Number Of Attributes: 12
> > reference : function rmMonAnnCycLLT in contrib.ncl
> > long_name : Monthly 1 degree resolution SST
> > standard_name : sea_surface_temperature
> > units : degC
> > add_offset : 0
> > cell_methods : time: lat: lon: mean
> > scale_factor : 1
> > _FillValue : -1e+30
> > missing_value : -1e+30
> > actual_range : ( -1.8, 34.76296 )
> > description : HadISST 1.1 monthly average sea surface temperature
> > anomaly_op_ncl : Annual Cycle Removed:function
> rmMonAnnCycLLT:contributed.ncl
> > (0)
> > (0) Monthly 1 degree resolution SST: min=-2.4631 max=2.59599
> >
> > Variable: precipanom_data_ts
> > Type: float
> > Total Size: 2928 bytes
> > 732 values
> > Number of Dimensions: 1
> > Dimensions and sizes: [time | 732]
> > Coordinates:
> > time: [194601..200612]
> > Number Of Attributes: 1
> > _FillValue : 9.96921e+36
> > (0)
> > (0) min=-1.73324 max=3.01711
> >
> > Variable: precipJJA1993
> > Type: float
> > Total Size: 4 bytes
> > 1 values
> > Number of Dimensions: 1
> > Dimensions and sizes: [time | 1]
> > Coordinates:
> > time: [199307..199307]
> > Number Of Attributes: 2
> > _FillValue : 9.96921e+36
> > NMO : 6
> > (0)
> > (0) min=0.267831 max=0.267831
> > warning:escorc: Non-fatal conditions encountered in series or xstd
> equals zero.
> > Possibly, all values of a series are constant.
> >
> > warning:escorc: Most likely, one or more series consisted of all
> constant values
> >
> > Variable: corr
> > Type: float
> > Total Size: 144000 bytes
> > 36000 values
> > Number of Dimensions: 2
> > Dimensions and sizes: [100] x [360]
> > Coordinates:
> > Number Of Attributes: 1
> > _FillValue : -1e+30
> > (0)
> > (0) min=-1e+30 max=-1e+30
> >
> > Regards,
> > Nkese.
> >
> > CONFIDENTIALITY: This email (including any attachments) may contain
> confidential, proprietary and/or privileged information. Any duplication,
> copying, distribution, dissemination, transmission, disclosure or use in
> any manner of this email (including any attachments) without the
> authorisation of the sender is strictly prohibited. If you receive this
> email (including any attachments) in error, please notify the sender and
> delete this email (including any attachments) from your system. Thank you.
> >
> >
> > CONFIDENTIALITY: This email (including any attachments) may contain
> confidential, proprietary and/or privileged information. Any duplication,
> copying, distribution, dissemination, transmission, disclosure or use in
> any manner of this email (including any attachments) without the
> authorisation of the sender is strictly prohibited. If you receive this
> email (including any attachments) in error, please notify the sender and
> delete this email (including any attachments) from your system. Thank you.
> >
> >
> >
> >
> > _______________________________________________
> > ncl-talk mailing list
> > List instructions, subscriber options, unsubscribe:
> > http://mailman.ucar.edu/mailman/listinfo/ncl-talk
> >
>
>
> ------------------------------
>
> Message: 9
> Date: Thu, 02 May 2013 06:35:49 -0600
> From: Dennis Shea <shea@ucar.edu>
> Subject: Re: How to set the distance among legend so that
> they do not overlap
> To: juki juki <juky_emc2@yahoo.com>
> Cc: ncl-talk <ncl-talk@ucar.edu>
> Message-ID: <51825DA5.7060702@ucar.edu>
> Content-Type: text/plain; charset=ISO-8859-1; format=flowed
>
> Change
> res@pmLegendHeightF = 0.08
> To a larger number
>
> and/or res@lgLabelFontHeightF = .012
>
> to a smaller number.
>
> You will have to experiment
>
>
> On 5/1/13 10:41 PM, juki juki wrote:
> > Hi all ncl users..
> >
> > Thanks for your help and nice discussion so far. I have plotted my data
> to XY plot with 16 legends. The problem is that the legends overlap each
> other. How to overcome such problem ? The following is simple ncl code to
> control my current legend:
> >
> > res@pmLegendDisplayMode = "Always" ; turn on legend
> > res@pmLegendSide = "Top" ; Change location of
> > res@pmLegendParallelPosF = .6 ; move units right
> > res@pmLegendOrthogonalPosF = -0.25 ; move units down
> > res@pmLegendWidthF = 0.06 ; Change width and
> > res@pmLegendHeightF = 0.08 ; height of legend.
> > res@lgPerimOn = False ; turn off box
> around
> > res@lgLabelFontHeightF = .012 ; label font height
> >
> > Thank you.
> > joeky
> >
> >
> >
> > _______________________________________________
> > ncl-talk mailing list
> > List instructions, subscriber options, unsubscribe:
> > http://mailman.ucar.edu/mailman/listinfo/ncl-talk
> >
>
>
> ------------------------------
>
> Message: 10
> Date: Thu, 2 May 2013 14:48:23 +0000
> From: Nkese Mc Shine <Nkese.McShine@sta.uwi.edu>
> Subject: Re: Correlation Question
> To: Dennis Shea <shea@ucar.edu>
> Cc: "ncl-talk@ucar.edu" <ncl-talk@ucar.edu>
> Message-ID:
> <7368311F6A984441893E050C8422DDBA3126D6AB@exMS01vm.sauwi.uwi.tt>
> Content-Type: text/plain; charset="us-ascii"
>
> Dear All,
>
> I probably was not clear on what I wanted to do. I would like to correlate
> one(1) month sst(lat,lon,time) with 1 month of precipitation which is
> actually one(1) value precip(time) and I wanted the correlation to be
> (lat,lon) so that it can be plotted on a map. So what I understand is this
> can not be done. Since my precipitation data has only the time dimension
> then using the pattern_cor() can not be used either because the output is
> (time) only and my precipitation data would have to also contain lat and
> lon.
>
> If there any other suggestions to what I may be able to do, please let me
> know. Thank you for all your help.
>
> Regards,
> Nkese.
> ________________________________________
> From: Dennis Shea [shea@ucar.edu]
> Sent: 02 May 2013 08:30
> To: Nkese Mc Shine
> Cc: ncl-talk@ucar.edu
> Subject: Re: Correlation Question
>
> Ho Nkese
>
> If you have an array: x(lat,lon,time)
>
> You can use the named dimensions to reorder
>
> xNew = x(time|:,lat|:,lon|:)
> printVarSummary(xNew)
>
> Then you can compute the pattern correaltion(s)
>
> Good Luck
>
> On 5/2/13 3:22 AM, Nkese Mc Shine wrote:
> > Will,
> >
> > Thank you for the explanation and suggestion however I do not think that
> will work as pattern_cor the rightmost dimension must be lon, and my
> monthly data which is 1 value is in time.
> >
> > Thank you.
> > Nkese.
> > ________________________________
> > From: Will Hobbs [Will.Hobbs@utas.edu.au]
> > Sent: 30 April 2013 19:37
> > To: Nkese Mc Shine; ncl-talk@ucar.edu
> > Subject: Re: Correlation Question
> >
> > Nkese
> >
> >> I would like to make a month to month correlation. I would like to
> correlate 1 month of sst(lat,lon,time) with 1 month of precipitation which
> actually is 1 value. Is it possible?
> >
> > No, you can't correlate a single value (in this case a month). This is
> nothing to do with NCL's capability, but fundamental to what a correlation
> is, i.e. a comparison of two variances. A single realisation (in your case
> a month) does not have a variance.
> >
> > What you can do (in case I misunderstood your question) is correlate the
> spatial pattern with that of a a single month; the output would be a time
> series of correlations. This can be done in NCL using the function
> 'pattern_cor()'
> >
> > Will
> >
> > From: Nkese Mc Shine <Nkese.McShine@sta.uwi.edu<mailto:
> Nkese.McShine@sta.uwi.edu>>
> > Date: Wednesday, 1 May 2013 12:11 AM
> > To: "ncl-talk@ucar.edu<mailto:ncl-talk@ucar.edu>" <ncl-talk@ucar.edu
> <mailto:ncl-talk@ucar.edu>>
> > Subject: Correlation Question
> >
> > Dear All,
> >
> > I would like to make a month to month correlation. I would like to
> correlate 1 month of sst(lat,lon,time) with 1 month of precipitation which
> actually is 1 value. Is it possible? Attached is my script and below is the
> output. Please let me know, what I can do if it is possible.
> >
> > Variable: sst_anom
> > Type: float
> > Total Size: 105408000 bytes
> > 26352000 values
> > Number of Dimensions: 3
> > Dimensions and sizes: [lat | 100] x [lon | 360] x [time | 732]
> > Coordinates:
> > lat: [49.5..-49.5]
> > lon: [-179.5..179.5]
> > time: [194601..200612]
> > Number Of Attributes: 12
> > anomaly_op_ncl : Annual Cycle Removed:function
> rmMonAnnCycLLT:contributed.ncl
> > description : HadISST 1.1 monthly average sea surface temperature
> > actual_range : ( -1.8, 34.76296 )
> > missing_value : -1e+30
> > _FillValue : -1e+30
> > scale_factor : 1
> > cell_methods : time: lat: lon: mean
> > add_offset : 0
> > units : degC
> > standard_name : sea_surface_temperature
> > long_name : Monthly 1 degree resolution SST
> > reference : function rmMonAnnCycLLT in contrib.ncl
> > (0)
> > (0) Monthly 1 degree resolution SST: min=-8.70077 max=6.98223
> >
> > Variable: sst_anom_jan
> > Type: float
> > Total Size: 144000 bytes
> > 36000 values
> > Number of Dimensions: 3
> > Dimensions and sizes: [lat | 100] x [lon | 360] x [time | 1]
> > Coordinates:
> > lat: [49.5..-49.5]
> > lon: [-179.5..179.5]
> > time: [199301..199301]
> > Number Of Attributes: 12
> > reference : function rmMonAnnCycLLT in contrib.ncl
> > long_name : Monthly 1 degree resolution SST
> > standard_name : sea_surface_temperature
> > units : degC
> > add_offset : 0
> > cell_methods : time: lat: lon: mean
> > scale_factor : 1
> > _FillValue : -1e+30
> > missing_value : -1e+30
> > actual_range : ( -1.8, 34.76296 )
> > description : HadISST 1.1 monthly average sea surface temperature
> > anomaly_op_ncl : Annual Cycle Removed:function
> rmMonAnnCycLLT:contributed.ncl
> > (0)
> > (0) Monthly 1 degree resolution SST: min=-2.4631 max=2.59599
> >
> > Variable: precipanom_data_ts
> > Type: float
> > Total Size: 2928 bytes
> > 732 values
> > Number of Dimensions: 1
> > Dimensions and sizes: [time | 732]
> > Coordinates:
> > time: [194601..200612]
> > Number Of Attributes: 1
> > _FillValue : 9.96921e+36
> > (0)
> > (0) min=-1.73324 max=3.01711
> >
> > Variable: precipJJA1993
> > Type: float
> > Total Size: 4 bytes
> > 1 values
> > Number of Dimensions: 1
> > Dimensions and sizes: [time | 1]
> > Coordinates:
> > time: [199307..199307]
> > Number Of Attributes: 2
> > _FillValue : 9.96921e+36
> > NMO : 6
> > (0)
> > (0) min=0.267831 max=0.267831
> > warning:escorc: Non-fatal conditions encountered in series or xstd
> equals zero.
> > Possibly, all values of a series are constant.
> >
> > warning:escorc: Most likely, one or more series consisted of all
> constant values
> >
> > Variable: corr
> > Type: float
> > Total Size: 144000 bytes
> > 36000 values
> > Number of Dimensions: 2
> > Dimensions and sizes: [100] x [360]
> > Coordinates:
> > Number Of Attributes: 1
> > _FillValue : -1e+30
> > (0)
> > (0) min=-1e+30 max=-1e+30
> >
> > Regards,
> > Nkese.
> >
> > CONFIDENTIALITY: This email (including any attachments) may contain
> confidential, proprietary and/or privileged information. Any duplication,
> copying, distribution, dissemination, transmission, disclosure or use in
> any manner of this email (including any attachments) without the
> authorisation of the sender is strictly prohibited. If you receive this
> email (including any attachments) in error, please notify the sender and
> delete this email (including any attachments) from your system. Thank you.
> >
> >
> > CONFIDENTIALITY: This email (including any attachments) may contain
> confidential, proprietary and/or privileged information. Any duplication,
> copying, distribution, dissemination, transmission, disclosure or use in
> any manner of this email (including any attachments) without the
> authorisation of the sender is strictly prohibited. If you receive this
> email (including any attachments) in error, please notify the sender and
> delete this email (including any attachments) from your system. Thank you.
> >
> >
> >
> >
> > _______________________________________________
> > ncl-talk mailing list
> > List instructions, subscriber options, unsubscribe:
> > http://mailman.ucar.edu/mailman/listinfo/ncl-talk
> >
>
>
> CONFIDENTIALITY: This email (including any attachments) may contain
> confidential, proprietary and/or privileged information. Any duplication,
> copying, distribution, dissemination, transmission, disclosure or use in
> any manner of this email (including any attachments) without the
> authorisation of the sender is strictly prohibited. If you receive this
> email (including any attachments) in error, please notify the sender and
> delete this email (including any attachments) from your system. Thank you.
>
>
>
>
> ------------------------------
>
> Message: 11
> Date: Thu, 2 May 2013 09:03:45 -0600
> From: Mary Haley <haley@ucar.edu>
> Subject: Re: How to set the distance among legend so that
> they do not overlap
> To: juki juki <juky_emc2@yahoo.com>
> Cc: ncl-talk <ncl-talk@ucar.edu>
> Message-ID: <D21BAE03-A7EB-4A01-8DD0-CBD7CB10D7BE@ucar.edu>
> Content-Type: text/plain; charset=us-ascii
>
> You can visit our legend page for more help:
>
> http://www.ncl.ucar.edu/Applications/legend.shtml
>
> As Dennis said, you need to experiment with values, depending on what you
> want.
>
> If you continue to have problems, it will help to see what your image
> looks like.
>
>
> --Mary
>
>
> On May 2, 2013, at 6:35 AM, Dennis Shea wrote:
>
> > Change
> > res@pmLegendHeightF = 0.08
> > To a larger number
> >
> > and/or res@lgLabelFontHeightF = .012
> >
> > to a smaller number.
> >
> > You will have to experiment
> >
> >
> > On 5/1/13 10:41 PM, juki juki wrote:
> >> Hi all ncl users..
> >>
> >> Thanks for your help and nice discussion so far. I have plotted my data
> to XY plot with 16 legends. The problem is that the legends overlap each
> other. How to overcome such problem ? The following is simple ncl code to
> control my current legend:
> >>
> >> res@pmLegendDisplayMode = "Always" ; turn on legend
> >> res@pmLegendSide = "Top" ; Change location of
> >> res@pmLegendParallelPosF = .6 ; move units right
> >> res@pmLegendOrthogonalPosF = -0.25 ; move units down
> >> res@pmLegendWidthF = 0.06 ; Change width and
> >> res@pmLegendHeightF = 0.08 ; height of legend.
> >> res@lgPerimOn = False ; turn off box
> around
> >> res@lgLabelFontHeightF = .012 ; label font height
> >>
> >> Thank you.
> >> joeky
> >>
> >>
> >>
> >> _______________________________________________
> >> 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: 12
> Date: Thu, 02 May 2013 09:46:42 -0600
> From: Dennis Shea <shea@ucar.edu>
> Subject: Re: Plotting skewt from RUC grib files?
> To: Gerry Creager - NOAA Affiliate <gerry.creager@noaa.gov>
> Cc: ncl-talk <ncl-talk@ucar.edu>
> Message-ID: <51828A62.3020003@ucar.edu>
> Content-Type: text/plain; charset=ISO-8859-1; format=flowed
>
> Gerry, see Example 6 at
> http://www.ncl.ucar.edu/Applications/skewt.shtml
> Cheers!
>
> On 5/1/13 4:19 PM, Gerry Creager - NOAA Affiliate wrote:
> > I'm trying to sort through this exercise. Has anyone else done this
> already
> > before I reinvent the wheel?
> >
> > Thanks, Gerry
> >
> >
> >
> > _______________________________________________
> > ncl-talk mailing list
> > List instructions, subscriber options, unsubscribe:
> > http://mailman.ucar.edu/mailman/listinfo/ncl-talk
> >
>
>
> ------------------------------
>
> Message: 13
> Date: Thu, 2 May 2013 10:53:20 -0500
> From: Gerry Creager - NOAA Affiliate <gerry.creager@noaa.gov>
> Subject: Re: Plotting skewt from RUC grib files?
> To: Dennis Shea <shea@ucar.edu>
> Cc: ncl-talk <ncl-talk@ucar.edu>
> Message-ID:
> <
> CAFufYD7OXQORx17smaLYEJO+KOTnP_Ko9-_OFHdmnPy1j1dqgQ@mail.gmail.com>
> Content-Type: text/plain; charset="windows-1252"
>
> Thanks, Dennis. I'll try that now. I was certain there was an example but I
> overlooked the obvious one. And I really did look!
>
> gerry
>
>
> On Thu, May 2, 2013 at 10:46 AM, Dennis Shea <shea@ucar.edu> wrote:
>
> > Gerry, see Example 6 at
> > http://www.ncl.ucar.edu/**Applications/skewt.shtml<
> http://www.ncl.ucar.edu/Applications/skewt.shtml>
> > Cheers!
> >
> >
> > On 5/1/13 4:19 PM, Gerry Creager - NOAA Affiliate wrote:
> >
> >> I'm trying to sort through this exercise. Has anyone else done this
> >> already
> >> before I reinvent the wheel?
> >>
> >> Thanks, Gerry
> >>
> >>
> >>
> >> ______________________________**_________________
> >> ncl-talk mailing list
> >> List instructions, subscriber options, unsubscribe:
> >> http://mailman.ucar.edu/**mailman/listinfo/ncl-talk<
> 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)
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL:
> http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20130502/a1b72a32/attachment.html
>
> ------------------------------
>
> Message: 14
> Date: Thu, 2 May 2013 16:54:00 +0000
> From: Douglas Lowe <Douglas.Lowe@manchester.ac.uk>
> Subject: How to fix plot size and shape independent of
> paper size
> To: "ncl-talk@ucar.edu" <ncl-talk@ucar.edu>
> Message-ID:
> <43E6B083008E774B87C4283E0FFA4E70010117F6EF@MBXP02.ds.man.ac.uk>
> Content-Type: text/plain; charset="us-ascii"
>
> Hi,
>
> I'm working on plotting wrf data using ncl with the purpose of creating
> semi-transparent png images for layering one on top of another. For this we
> need the plot size to always be the same size and shape - which is
> difficult to do if it is set using the paper size (and can change if the
> title string is too long, or the colorbar labels are different lengths). Is
> there anyway with ncl to set the plot size independently of the paper size?
>
> Thanks,
> Doug
>
> ------------------------------
>
> Message: 15
> Date: Thu, 02 May 2013 10:57:28 -0600
> From: Dennis Shea <shea@ucar.edu>
> Subject: Re: Correlation Question
> To: Nkese Mc Shine <Nkese.McShine@sta.uwi.edu>
> Cc: "ncl-talk@ucar.edu" <ncl-talk@ucar.edu>
> Message-ID: <51829AF8.4070800@ucar.edu>
> Content-Type: text/plain; charset=ISO-8859-1; format=flowed
>
> Maybe 'we' are misunderstanding you but, by definition,
> a correlation is calculated via:
>
> cor = SUM [(X(t)-Xave)*(Y(t)-Yave)}]/(Xstd*Ystd)
>
> so if X=>sst(lat,lon,time)
> and Y=>prc(time)
>
> where time is size one, Y(t)-Yave is zero because
> the Yave=prc. The avergae of one value is the mean.
>
> Further, the Ystd is 0 since there is only one value.
>
> same for X if time is size one.
>
> On 5/2/13 8:48 AM, Nkese Mc Shine wrote:
> > Dear All,
> >
> > I probably was not clear on what I wanted to do. I would like to
> correlate one(1) month sst(lat,lon,time) with 1 month of precipitation
> which is actually one(1) value precip(time) and I wanted the correlation to
> be (lat,lon) so that it can be plotted on a map. So what I understand is
> this can not be done. Since my precipitation data has only the time
> dimension then using the pattern_cor() can not be used either because the
> output is (time) only and my precipitation data would have to also contain
> lat and lon.
> >
> > If there any other suggestions to what I may be able to do, please let
> me know. Thank you for all your help.
> >
> > Regards,
> > Nkese.
> > ________________________________________
> > From: Dennis Shea [shea@ucar.edu]
> > Sent: 02 May 2013 08:30
> > To: Nkese Mc Shine
> > Cc: ncl-talk@ucar.edu
> > Subject: Re: Correlation Question
> >
> > Ho Nkese
> >
> > If you have an array: x(lat,lon,time)
> >
> > You can use the named dimensions to reorder
> >
> > xNew = x(time|:,lat|:,lon|:)
> > printVarSummary(xNew)
> >
> > Then you can compute the pattern correaltion(s)
> >
> > Good Luck
> >
> > On 5/2/13 3:22 AM, Nkese Mc Shine wrote:
> >> Will,
> >>
> >> Thank you for the explanation and suggestion however I do not think
> that will work as pattern_cor the rightmost dimension must be lon, and my
> monthly data which is 1 value is in time.
> >>
> >> Thank you.
> >> Nkese.
> >> ________________________________
> >> From: Will Hobbs [Will.Hobbs@utas.edu.au]
> >> Sent: 30 April 2013 19:37
> >> To: Nkese Mc Shine; ncl-talk@ucar.edu
> >> Subject: Re: Correlation Question
> >>
> >> Nkese
> >>
> >>> I would like to make a month to month correlation. I would like to
> correlate 1 month of sst(lat,lon,time) with 1 month of precipitation which
> actually is 1 value. Is it possible?
> >>
> >> No, you can't correlate a single value (in this case a month). This is
> nothing to do with NCL's capability, but fundamental to what a correlation
> is, i.e. a comparison of two variances. A single realisation (in your case
> a month) does not have a variance.
> >>
> >> What you can do (in case I misunderstood your question) is correlate
> the spatial pattern with that of a a single month; the output would be a
> time series of correlations. This can be done in NCL using the function
> 'pattern_cor()'
> >>
> >> Will
> >>
> >> From: Nkese Mc Shine <Nkese.McShine@sta.uwi.edu<mailto:
> Nkese.McShine@sta.uwi.edu>>
> >> Date: Wednesday, 1 May 2013 12:11 AM
> >> To: "ncl-talk@ucar.edu<mailto:ncl-talk@ucar.edu>" <ncl-talk@ucar.edu
> <mailto:ncl-talk@ucar.edu>>
> >> Subject: Correlation Question
> >>
> >> Dear All,
> >>
> >> I would like to make a month to month correlation. I would like to
> correlate 1 month of sst(lat,lon,time) with 1 month of precipitation which
> actually is 1 value. Is it possible? Attached is my script and below is the
> output. Please let me know, what I can do if it is possible.
> >>
> >> Variable: sst_anom
> >> Type: float
> >> Total Size: 105408000 bytes
> >> 26352000 values
> >> Number of Dimensions: 3
> >> Dimensions and sizes: [lat | 100] x [lon | 360] x [time | 732]
> >> Coordinates:
> >> lat: [49.5..-49.5]
> >> lon: [-179.5..179.5]
> >> time: [194601..200612]
> >> Number Of Attributes: 12
> >> anomaly_op_ncl : Annual Cycle Removed:function
> rmMonAnnCycLLT:contributed.ncl
> >> description : HadISST 1.1 monthly average sea surface temperature
> >> actual_range : ( -1.8, 34.76296 )
> >> missing_value : -1e+30
> >> _FillValue : -1e+30
> >> scale_factor : 1
> >> cell_methods : time: lat: lon: mean
> >> add_offset : 0
> >> units : degC
> >> standard_name : sea_surface_temperature
> >> long_name : Monthly 1 degree resolution SST
> >> reference : function rmMonAnnCycLLT in contrib.ncl
> >> (0)
> >> (0) Monthly 1 degree resolution SST: min=-8.70077 max=6.98223
> >>
> >> Variable: sst_anom_jan
> >> Type: float
> >> Total Size: 144000 bytes
> >> 36000 values
> >> Number of Dimensions: 3
> >> Dimensions and sizes: [lat | 100] x [lon | 360] x [time | 1]
> >> Coordinates:
> >> lat: [49.5..-49.5]
> >> lon: [-179.5..179.5]
> >> time: [199301..199301]
> >> Number Of Attributes: 12
> >> reference : function rmMonAnnCycLLT in contrib.ncl
> >> long_name : Monthly 1 degree resolution SST
> >> standard_name : sea_surface_temperature
> >> units : degC
> >> add_offset : 0
> >> cell_methods : time: lat: lon: mean
> >> scale_factor : 1
> >> _FillValue : -1e+30
> >> missing_value : -1e+30
> >> actual_range : ( -1.8, 34.76296 )
> >> description : HadISST 1.1 monthly average sea surface temperature
> >> anomaly_op_ncl : Annual Cycle Removed:function
> rmMonAnnCycLLT:contributed.ncl
> >> (0)
> >> (0) Monthly 1 degree resolution SST: min=-2.4631 max=2.59599
> >>
> >> Variable: precipanom_data_ts
> >> Type: float
> >> Total Size: 2928 bytes
> >> 732 values
> >> Number of Dimensions: 1
> >> Dimensions and sizes: [time | 732]
> >> Coordinates:
> >> time: [194601..200612]
> >> Number Of Attributes: 1
> >> _FillValue : 9.96921e+36
> >> (0)
> >> (0) min=-1.73324 max=3.01711
> >>
> >> Variable: precipJJA1993
> >> Type: float
> >> Total Size: 4 bytes
> >> 1 values
> >> Number of Dimensions: 1
> >> Dimensions and sizes: [time | 1]
> >> Coordinates:
> >> time: [199307..199307]
> >> Number Of Attributes: 2
> >> _FillValue : 9.96921e+36
> >> NMO : 6
> >> (0)
> >> (0) min=0.267831 max=0.267831
> >> warning:escorc: Non-fatal conditions encountered in series or xstd
> equals zero.
> >> Possibly, all values of a series are constant.
> >>
> >> warning:escorc: Most likely, one or more series consisted of all
> constant values
> >>
> >> Variable: corr
> >> Type: float
> >> Total Size: 144000 bytes
> >> 36000 values
> >> Number of Dimensions: 2
> >> Dimensions and sizes: [100] x [360]
> >> Coordinates:
> >> Number Of Attributes: 1
> >> _FillValue : -1e+30
> >> (0)
> >> (0) min=-1e+30 max=-1e+30
> >>
> >> Regards,
> >> Nkese.
> >>
> >> CONFIDENTIALITY: This email (including any attachments) may contain
> confidential, proprietary and/or privileged information. Any duplication,
> copying, distribution, dissemination, transmission, disclosure or use in
> any manner of this email (including any attachments) without the
> authorisation of the sender is strictly prohibited. If you receive this
> email (including any attachments) in error, please notify the sender and
> delete this email (including any attachments) from your system. Thank you.
> >>
> >>
> >> CONFIDENTIALITY: This email (including any attachments) may contain
> confidential, proprietary and/or privileged information. Any duplication,
> copying, distribution, dissemination, transmission, disclosure or use in
> any manner of this email (including any attachments) without the
> authorisation of the sender is strictly prohibited. If you receive this
> email (including any attachments) in error, please notify the sender and
> delete this email (including any attachments) from your system. Thank you.
> >>
> >>
> >>
> >>
> >> _______________________________________________
> >> ncl-talk mailing list
> >> List instructions, subscriber options, unsubscribe:
> >> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
> >>
> >
> >
> > CONFIDENTIALITY: This email (including any attachments) may contain
> confidential, proprietary and/or privileged information. Any duplication,
> copying, distribution, dissemination, transmission, disclosure or use in
> any manner of this email (including any attachments) without the
> authorisation of the sender is strictly prohibited. If you receive this
> email (including any attachments) in error, please notify the sender and
> delete this email (including any attachments) from your system. Thank you.
> >
> >
>
>
> ------------------------------
>
> Message: 16
> Date: Thu, 02 May 2013 13:07:02 -0400
> From: mmkamal@uwaterloo.ca
> Subject: How to control "correlation" text position in
> Taylor Diagram.
> To: "ncl-talk@ucar.edu" <ncl-talk@ucar.edu>
> Message-ID:
> <20130502130702.43044c0d7sd28zs4@www.nexusmail.uwaterloo.ca>
> Content-Type: text/plain; charset="iso-8859-1"
>
> Hi,
>
> I have increased the values of rxy@tmXBValues rxy@tmXBLabels
> rxy@tmYLValues rxy@tmYLLabels from default 1.65 to 3.00 but failed to
> change the position of the text "correlation". Now the "correlation"
> appeared inside the circle. I am attaching the plot for your
> convenience. Another issue is the RMSD circle does not end at the
> X-axis.Could anyone please suggest me how to control the position of
> the text.
>
>
>
> Thanks
> Kamal
>
>
>
>
>
>
> -------------- next part --------------
> A non-text attachment was scrubbed...
> Name: Prcp_July_Taylor.4.png
> Type: image/png
> Size: 136041 bytes
> Desc: not available
> Url :
> http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20130502/addc3675/attachment.png
>
> ------------------------------
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk@ucar.edu
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
> End of ncl-talk Digest, Vol 114, Issue 2
> ****************************************
>

-- 
*Jayasankar C B*
*
*
*Project Assistant*
*CSIR-CMMACS*
*NAL Bhelur campus*
*Bangalore-37.*

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Fri May 3 04:50:12 2013

This archive was generated by hypermail 2.1.8 : Tue May 07 2013 - 08:38:33 MDT