Taylor diagrams for model/data deep Pacific water 14C intercomparison


At a glance

You can download a PDF version of each diagram below by clicking on it.

See the commands that produce these diagrams

 Main Taylor diagram 

 Multi-component Taylor diagram 

  Total space-time field
  Multi-component diagram

The commands that produce these diagrams:

1st step: The command that computes the statistics NetCDF file (C14/comp_stat.nc) 

>>> python model_vs_data_comp_stat.py -oC14/comp_stat.nc -tC14/comp_stat.txt data/KEY_C14.nc C14 "model_output/%(model).nc" DC14

The four last parameters are, respectively:
data/KEY_C14.nc:
NetCDF file  containing estimates or observation data. This file is derived from data-based annual mean 14C in Pacific Ocean provided by Bob Key.
C14:
Variable id within this file
model_output/%(model).nc:
Template for NetCDF files containing model output. Note that this parameter is quoted to protect it from the shell since it contains the character '%' (with special meaning in shell syntax)
DC14:
Variable id common to all model output files

Due to the use of option '-t', this command also produces a text file containing all the statistics (C14/comp_stat.txt).

2nd step: The commands that read the statistics file and produce the Taylor diagrams, above

Main Taylor diagram

>>> python plot_taylor_diagrams.py -i38 -C1,8,2 -t "Deep Pacific water C14 Taylor diagram" -pbrd -o C14/comp_stat/ C14/comp_stat/comp_stat.nc


This command results in 3 diagrams (in 3 files): the main one, as above, plus 2 for different statistical components.

Option '-i' sets the size of marker or symbol for each models on the diagram (allowed size ranges from 0 to 100)

Option -C selects a subset of statistical components. Since observation data and model output are annual means, only spatial components are selected.

Multi-component Taylor diagram

>>> python plot_taylor_diagrams.py -i32 -t "Deep Pacific water C14 Taylor diagram" -sp -C1,8,2 -o C14/comp_stat/multi_component.ps C14/comp_stat/comp_stat.nc


The same set of statistical components as on the left is selected.

Option '-o' sets the output Postscript file.

Note that, though we specify the diagram title, some other meta-information (units) is extracted from the data or model output files.
The model list and their code letter are extracted from a miscellaneous file: model_dictionary.py.