statistics for Taylor Diagrams

From: FERMEPIN Solange <solange.fermepin_at_nyahnyahspammersnyahnyah>
Date: Thu May 16 2013 - 08:50:06 MDT

Hello everybody,

I am trying to make some Taylor Diagrams. I was following the example
taylor_3.ncl while using the statistics defined in taylor_7.ncl.

I was getting ~0 correlation coeficient fot only two variables in my
timeseries and reasonable values for all the others.

I wrote again the pattern correlation by following "Taylor Diagram
Primer" document posted in the site.

I re run the script and now I get reasonable values (where before I got
~0) and exactly the same values for all the other variabales.

Does anyone have any idea why? Is there any statistical reason that I am

I am using ncl version 6.1.2. Below the relevant parts of my code:

Thanks a lot!

;;; Statistics part of the script

  ratio = new((/ncase, nvar/),float)
  cc = new((/ncase, nvar/),float)

  gcm = new((/ntim,nlat,nlon/),float)
  obs = new((/ntim,nlat,nlon/),float)
  do nc=0,ncase-1
   do nv=0,nvar-1

    gcm = all_gcm(nc,nv,:,:,:)
    obs = all_obs(nc,nv,:,:,:)

    tot_area = sum(aire)
    tot_gcm = sum(aire*gcm)
    tot_obs = sum(aire*obs)
    mean_gcm = tot_gcm/tot_area
    mean_obs = tot_obs/tot_area

    wvar_gcm = sum(aire*(gcm-mean_gcm)^2)/tot_area
    wvar_obs = sum(aire*(obs-mean_obs)^2)/tot_area

    wcc_gcm = (sum(aire*gcm*obs)-tot_gcm*tot_obs/tot_area)/ \

    ;wcc_gcm = sum(aire*(gcm-mean_gcm)*(obs-mean_obs))/tot_area/ \
    ; (wvar_gcm*wvar_obs)^0.5

    ratio(nc,nv) = (wvar_gcm/wvar_obs)^0.5
    cc(nc,nv) = wcc_gcm

   end do
  end do

Solange Fermepin, PhD student
Laboratoire de Météorologie Dynamique - CNRS/IPSL
Tour 45-55, 3ème étage
4 place Jussieu
75252 Paris cedex 05
Tel: +33 (0) 1 44 27 52 77
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
Received on Thu May 16 08:51:10 2013

This archive was generated by hypermail 2.1.8 : Thu May 30 2013 - 11:38:10 MDT