In [14]:

    

#!/usr/bin/python
# _*_ coding: latin-1 -*-
# Taylor Diagram - Based on Taylor (2001) - Journal Geophysical Research
# author: Mabel Calim Costa
# GMAO - INPE
# 20/02/2018

import numpy as np
from numpy import ma
import mpl_toolkits.axisartist.grid_finder as GF
import mpl_toolkits.axisartist.floating_axes as FA
import matplotlib.pyplot as plt
import netCDF4


    

In [15]:

    

import taylor_diag


    

In [16]:

    

import numpy as np
from pylab import *
z = np.linspace(0,2048,2048)
x = np.sin(50*np.pi*z)
y = np.cos(50*np.pi*z)
x2 = 2*np.sin(50*np.pi*z)

plot(x, 'k--')
plot(x2)
plot(y)


    

    
Out[16]:





[<matplotlib.lines.Line2D at 0x112d25ef0>]






    

In [17]:

    

series ={}
series[0] = x  # the first term will be always the reference signal 
series[1]= x2
series[2]= y


    

In [21]:

    

from matplotlib.projections import PolarAxes

taylor_diag.Taylor_diag(series,['x','x2','y'])


    

    




/Users/calim/Taylor_diagram/taylor_diag.py:40: RuntimeWarning: invalid value encountered in arccos
  tlocs = np.arccos(rlocs)        # Conversion to polar angles






    




---------------------------------------------------------------------------
NameError                                 Traceback (most recent call last)
<ipython-input-21-9740b3606668> in <module>()
      1 from matplotlib.projections import PolarAxes
      2 
----> 3 taylor_diag.Taylor_diag(series,['x','x2','y'])

~/Taylor_diagram/taylor_diag.py in Taylor_diag(series, names)
     52 
     53 
---> 54     tr = PolarAxes.PolarTransform()
     55 
     56     smin = 0

NameError: name 'PolarAxes' is not defined

In [20]:

    

print (taylor_diag.Taylor_diag.__doc__)


    

    




 Taylor Diagram : obs is reference data sample
        in a full diagram (0 --> npi)
        --------------------------------------------------------------------------
        Input: series     - dict with all time series (lists) to analyze  
               series[0]  - is the observation, the reference by default.
    

In [ ]: