NCL_taylor_3.py

This script illustrates the following concepts:

• Creating a basic Taylor diagram using geocat-viz Taylor diagram function.

• Adding labels to a Taylor diagram

See following URLs to see the reproduced NCL plot & script:

• Original NCL script: https://www.ncl.ucar.edu/Applications/Scripts/taylor_3.ncl

• Original NCL plot: https://www.ncl.ucar.edu/Applications/Images/taylor_3_lg.png

Note: Due to limitations of matplotlib’s axisartist toolkit, we cannot include minor tick marks

between 0.9 and 0.99, as seen in the original NCL plot.

Import packages:

import matplotlib.pyplot as plt

import geocat.viz as gv

Create dummy data:

# Case A
CA_ratio = [
    1.230,
    0.988,
    1.092,
    1.172,
    1.064,
    0.966,
    1.079,
    0.781,
]  # standard deviation
CA_cc = [
    0.958,
    0.973,
    0.740,
    0.743,
    0.922,
    0.982,
    0.952,
    0.433,
]  # correlation coefficient

# Case B
CB_ratio = [
    1.129,
    0.996,
    1.016,
    1.134,
    1.023,
    0.962,
    1.048,
    0.852,
]  # standard deviation
CB_cc = [
    0.963,
    0.975,
    0.801,
    0.814,
    0.946,
    0.984,
    0.968,
    0.647,
]  # correlation coefficient

Plot:

# Create figure and TaylorDiagram instance
fig = plt.figure(figsize=(10, 10))
dia = gv.TaylorDiagram(fig=fig, label='REF')

# Add models to Taylor diagram
dia.add_model_set(CA_ratio, CA_cc, color='red', marker='o', label='Case A', fontsize=16)

dia.add_model_set(
    CB_ratio, CB_cc, color='blue', marker='o', label='Case B', fontsize=16
)

# Create model name list
namearr = ['SLP', 'Tsfc', 'Prc', 'Prc 30S-30N', 'LW', 'SW', 'U300', 'Guess']

# Add model name
dia.add_model_name(namearr, fontsize=16)

# Add figure legend
dia.add_legend(fontsize=16)

# Show the plot
plt.show()