[ncl-talk] RATIO in Taylor Diagram

[Dennis Shea]

See attached: Hopefully that will answer your question.

---

Please do not send repeat questions to ncl-talk
All questions are flagged and will be answered. Many of us are on
travel or vacation.

Good Luck
D,


On Mon, Sep 14, 2015 at 8:14 PM, Vanúcia Schumacher
<vanucia-schumacher at hotmail.com> wrote:
> Dear users,
>
> I'm confused with the calculation of RATIO (normalized root-mean-saquare
> (RMS) differences), relation to Taylor Diagram,  the comments that I have
> read here in the group, it was suggested calculated in 3 ways:
>
> 1) Use of the function:   dim_rmsd
>
> or
>
> 2) RATIO
> ; temporal variance at each grid point [local]
>
> ;   vref_var_T = dim_variance_n(rdata, 0 ) ; (lat,lon)
> ;   vcase_var_T = dim_variance_n(cdata, 0 )
>
> ;                 ; wgted areal *local* temporal variance
>
> ;   wvar_ref_T = sum(wgt_S*vref_var_T)/sumwgt_S
> ;   wvar_case_T = sum(wgt_S*vcase_var_T)/sumwgt_S
>
>                  ; sqrt of ratio of spatially weighted variances
>
> ;   wvar_ratio_T = (wvar_case_T/wvar_ref_T)^0.5
>
> I did not understand what would that sumwgt_S ?
>
> or 3) more explicitly, for xc and xo on the same grid: wgtc=wgto
>
> ;         wgtc  = conform_dims(dimsizes(xc), gw, 0)    ; make 2d for gw[*]
> ;         xavgc = sum(wgtc*xc)/sum(wgtc)               ; control centered
> mean
> ;         xavgo = sum(wgtc*xo)/sum(wgtc)
>
> ;     (b) compute the sum of the centered area weighted variances.
>
> ;         dc2   = sum(wgtc*(xc-xavgc)^2)               ; control ; centered
> about xavgc
> ;         do2   = sum(wgtc*(xo-xavgo)^2)
> ;         rat   = sqrt(do2/dc2)
>
> xc and xo is the variance or datasets?
>
>
> I would like to know the step by step to correct this calculation RATIO, if
> someone can help me, please.
>
> Thanks
>
>
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: taylor_diag_stats.ncl
Type: application/octet-stream
Size: 3153 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20150917/2fe59582/attachment.obj