Hi Dennis,
can you show me how does NCL caculate the spec@coher_probability?
Thanks.
2013-09-18
朱学明
发件人:Dennis Shea <shea@ucar.edu>
发送时间:2013-09-18 07:29
主题:Re: [ncl-talk] confidence interval/level
收件人:"朱学明"<557xiaozhu@163.com>
抄送:"ncl-talk@ucar.edu USERS"<ncl-talk@ucar.edu>
spec = specxy_anal(x,y,d,sm,pct)
PHASE = spec@phase
PHASE@_FillValue = 1e20
PHASE = where(spec@coher_probability(1).le.0.95,PHASE@_FillValue,PHASE)
=========
See: http://www.ncl.ucar.edu/Applications/xy.shtml#ex25
The right (R) plot should
resR = True ; plot mods desired
resR@xyMarkLineModes = "Markers" ; choose which have markers
resR@xyMarkers = 16 ; choose type of marker
resR@xyMarkerColor = "red" ; Marker color
resR@xyMarkerSizeF = 0.005 ; Marker size (default 0.01)
plot = gsn_csm_xy2(wks,spec@frq,spec@coher,PHASE,resL,resR)
On 9/17/13 4:04 PM, 朱学明 wrote:
> Hi all,
> I have tried to use the function specxy_anal get the coherance and phase of two timeseries data based on the example script spec_2.ncl. But I just want to get the coherence and phase satisfy the 0.95 confidence level, the picture likes that in the attachment. How can I do that?
> It looks like that the output of function specx_ci is not what I want.
>
> Any help would be appreciated. Thanks.
>
> 2013-09-18
>
>
>
> 朱学明
>
>
>
> _______________________________________________
> 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 Wed Sep 18 08:54:43 2013
This archive was generated by hypermail 2.1.8 : Wed Sep 18 2013 - 17:07:30 MDT