Commit b4e536dd authored by Kennedy, Joseph H's avatar Kennedy, Joseph H
Browse files

Descriptive plot titles for all extensions

Fixes #20
parent 2cd6d2b4
......@@ -14,6 +14,8 @@ Kay and Gettelman, 2009). The black solid lines denote the 0, 1000, 2000, and
3000 meter Greenland ice sheet elevation contours as seen by CESM.
"""
title = "High cloud cover over Greenland"
def make_plot(config, out_path='.'):
# ---------------- Data source in TITAN ------------------------
......@@ -189,7 +191,7 @@ def make_plot(config, out_path='.'):
img_link = os.path.join(os.path.basename(out_path),
os.path.basename(wks_img + '.' + wks_type))
img_elem = el.image('CESM_cldhgh_ANN',
img_elem = el.image(title,
' '.join(describe.split()),
img_link)
img_elem['Height'] = config['image_height']
......
......@@ -14,6 +14,8 @@ Kay and Gettelman, 2009). The black solid lines denote the 0, 1000, 2000, and
3000 meter Greenland ice sheet elevation contours as seen by CESM.
"""
title = "Low cloud cover over Greenland"
def make_plot(config, out_path='.'):
# ---------------- Data source in TITAN ------------------------
......@@ -189,7 +191,7 @@ def make_plot(config, out_path='.'):
img_link = os.path.join(os.path.basename(out_path),
os.path.basename(wks_img + '.' + wks_type))
img_elem = el.image('CESM_cldlow_ANN',
img_elem = el.image(title,
' '.join(describe.split()),
img_link)
img_elem['Height'] = config['image_height']
......
......@@ -14,6 +14,8 @@ Kay and Gettelman, 2009). The black solid lines denote the 0, 1000, 2000, and
3000 meter Greenland ice sheet elevation contours as seen by CESM.
"""
title = "Cloud cover over Greenland"
def make_plot(config, out_path='.'):
# ---------------- Data source in TITAN ------------------------
......@@ -189,7 +191,7 @@ def make_plot(config, out_path='.'):
img_link = os.path.join(os.path.basename(out_path),
os.path.basename(wks_img + '.' + wks_type))
img_elem = el.image('CESM_cldtot_ANN',
img_elem = el.image(title,
' '.join(describe.split()),
img_link)
img_elem['Height'] = config['image_height']
......
......@@ -13,6 +13,8 @@ for CESM (red), ISCCP (green; Rossow and Schiffer, 1999), and CLOUDSAT (cyan;
Kay and Gettelman, 2009).
"""
title = "Monthly total high cloud cover over Greenland"
def make_plot(config, out_path='.'):
......@@ -68,7 +70,7 @@ def make_plot(config, out_path='.'):
img_link = os.path.join(os.path.basename(out_path),
os.path.basename(img_path))
img_elem = el.image('yearly_cycle_cldhgh',
img_elem = el.image(title,
' '.join(describe.split()),
img_link)
if config:
......
......@@ -13,6 +13,8 @@ for CESM (red), ISCCP (green; Rossow and Schiffer, 1999), and CLOUDSAT (cyan;
Kay and Gettelman, 2009).
"""
title = "Monthly total low cloud cover over Greenland"
def make_plot(config, out_path='.'):
......@@ -69,7 +71,7 @@ def make_plot(config, out_path='.'):
img_link = os.path.join(os.path.basename(out_path),
os.path.basename(img_path))
img_elem = el.image('yearly_cycle_cldlow',
img_elem = el.image(title,
' '.join(describe.split()),
img_link)
if config:
......
......@@ -13,6 +13,8 @@ for CESM (red), ISCCP (green; Rossow and Schiffer, 1999), and CLOUDSAT (cyan;
Kay and Gettelman, 2009).
"""
title = "Monthly total cloud cover over Greenland"
def make_plot(config, out_path='.'):
......@@ -69,7 +71,7 @@ def make_plot(config, out_path='.'):
img_link = os.path.join(os.path.basename(out_path),
os.path.basename(img_path))
img_elem = el.image('yearly_cycle_cldtot',
img_elem = el.image(title,
' '.join(describe.split()),
img_link)
if config:
......
......@@ -17,6 +17,8 @@ difference (CISM-A - init) over Jakobshavn glacier (West Greenland; Zwally et al
"""
# plots velnorm, zooms in on Jakobshavn 69.2N 51.1W
title = "Surface velocity norm over Jakobshavn glacier"
def make_plot(config, out_path='.'):
# ---------------- Read Data -----------------------------------
......@@ -168,7 +170,7 @@ def make_plot(config, out_path='.'):
img_link = os.path.join(os.path.basename(out_path),
os.path.basename(wks_img + '.' + wks_type))
img_elem = el.image('CISMA_js',
img_elem = el.image(title,
' '.join(describe.split()),
img_link)
img_elem['Height'] = config['image_height']
......
......@@ -17,6 +17,8 @@ difference (CISM-A - init) over Petermann glacier (North-West Greenland; Zwally
"""
# Plots velnorm, zooms in on Petermann glacier, 80.5N -59.5
title = "Surface velocity norm over Petermann glacier"
def make_plot(config, out_path='.'):
# ---------------- Read Data -----------------------------------
......@@ -168,7 +170,7 @@ def make_plot(config, out_path='.'):
img_link = os.path.join(os.path.basename(out_path),
os.path.basename(wks_img + '.' + wks_type))
img_elem = el.image('CISMA_pm',
img_elem = el.image(title,
' '.join(describe.split()),
img_link)
img_elem['Height'] = config['image_height']
......
......@@ -16,6 +16,8 @@ for RACMO 2.3 and CISM-Albany (after spin-up forced by RACMO 2.0) interpolated t
the coarser RACMO grid.
"""
title = "Surface mass balance comparison"
def make_plot(config, out_path='.'):
# --------------------------------------------------------------
......@@ -89,7 +91,7 @@ def make_plot(config, out_path='.'):
img_link = os.path.join(os.path.basename(out_path),
os.path.basename(img_path))
img_elem = el.image('scatter_CISMA_RACMO23',
img_elem = el.image(title,
' '.join(describe.split()),
img_link)
img_elem['Height'] = config['image_height']
......
......@@ -15,6 +15,8 @@ forced with RACMO 2.0 (left), RACMO 2.3 (right), and their difference
(CISM-Albany - RACMO 2.3) with CISM-Albany interpolated to the coarser RACMO 2.3 grid.
"""
title = "Surface mass balance over Greenland"
def make_plot(config, out_path='.'):
......@@ -179,7 +181,7 @@ def make_plot(config, out_path='.'):
# Creat multiple figures and draw, which now contains the elevation and temperature
# "[1,3]" indicates 1 row, 3 columns.
map_title = ["CISM-A","RACMO23", "CISM-A--RACMO23"]
map_title = ["CISM-A", "RACMO23", "CISM-A--RACMO23"]
nmap = 3
plot = []
......@@ -196,7 +198,7 @@ def make_plot(config, out_path='.'):
img_link = os.path.join(os.path.basename(out_path),
os.path.basename(wks_img + '.' + wks_type))
img_elem = el.image('CISMA_RACMO23',
img_elem = el.image(title,
' '.join(describe.split()),
img_link)
img_elem['Height'] = config['image_height']
......
......@@ -9,7 +9,12 @@ from netCDF4 import Dataset
from livvkit.util import elements as el
describe = """Surface temperature (°C) of the Greenland ice-sheet after CISM-Albany spin-up."""
describe = """
Surface temperature (°C) of the Greenland ice sheet after CISM-Albany spin-up.
"""
title = "Surface temperature of the Greenland ice sheet"
def make_plot(config, out_path='.'):
......@@ -152,7 +157,7 @@ def make_plot(config, out_path='.'):
img_link = os.path.join(os.path.basename(out_path),
os.path.basename(wks_img + '.' + wks_type))
img_elem = el.image('CISMA_temp',
img_elem = el.image(title,
' '.join(describe.split()),
img_link)
img_elem['Height'] = config['image_height']
......
......@@ -16,6 +16,8 @@ The norm of the surface velocity (m s^-1) from CISM initialization dataset
difference (CISM-A - init) over Greenland.
"""
title = "Surface velocity norm over Greenland"
def make_plot(config, out_path='.'):
# ---------------- Read Data -----------------------------------
......@@ -167,7 +169,7 @@ def make_plot(config, out_path='.'):
img_link = os.path.join(os.path.basename(out_path),
os.path.basename(wks_img + '.' + wks_type))
img_elem = el.image('CISMA_velnorm',
img_elem = el.image(title,
' '.join(describe.split()),
img_link)
img_elem['Height'] = config['image_height']
......
......@@ -17,6 +17,8 @@ difference (CISM-A - init) over North-East Greenland (Zwally et al. (2012) basin
"""
# zoomed in velnorm over Zachariae glacier at 78N, -30E
title = "Surface velocity norm over North-East Greenland"
def make_plot(config, out_path='.'):
# ---------------- Read Data -----------------------------------
......@@ -168,7 +170,7 @@ def make_plot(config, out_path='.'):
img_link = os.path.join(os.path.basename(out_path),
os.path.basename(wks_img + '.' + wks_type))
img_elem = el.image('CISMA_zach',
img_elem = el.image(title,
' '.join(describe.split()),
img_link)
img_elem['Height'] = config['image_height']
......
......@@ -18,6 +18,8 @@ Greenland ice sheet elevation contours as seen by the models (CESM's contours
shown in the difference plot).
"""
title = "Summer surface albedo"
def make_plot(config, out_path='.'):
# ---------------- Data source in TITAN ------------------------
......@@ -229,7 +231,7 @@ def make_plot(config, out_path='.'):
img_link = os.path.join(os.path.basename(out_path),
os.path.basename(wks_img + '.' + wks_type))
img_elem = el.image('CESM_RACMO23_albedo',
img_elem = el.image(title,
' '.join(describe.split()),
img_link)
img_elem['Height'] = config['image_height']
......
......@@ -19,6 +19,8 @@ Greenland ice sheet elevation contours as seen by the models (CESM's contours
shown in the difference plot).
"""
title = "Summer latent heat flux"
def make_plot(config, out_path='.'):
# ---------------- Data source in TITAN ------------------------
......@@ -228,7 +230,7 @@ def make_plot(config, out_path='.'):
img_link = os.path.join(os.path.basename(out_path),
os.path.basename(wks_img + '.' + wks_type))
img_elem = el.image('CESM_RACMO23_laft',
img_elem = el.image(title,
' '.join(describe.split()),
img_link)
img_elem['Height'] = config['image_height']
......
......@@ -19,6 +19,8 @@ Greenland ice sheet elevation contours as seen by the models (CESM's contours
shown in the difference plot).
"""
title = "Summer downward longwave radiation"
def make_plot(config, out_path='.'):
# ---------------- Data source in TITAN ------------------------
......@@ -230,7 +232,7 @@ def make_plot(config, out_path='.'):
img_link = os.path.join(os.path.basename(out_path),
os.path.basename(wks_img + '.' + wks_type))
img_elem = el.image('CESM_RACMO23_lwsd',
img_elem = el.image(title,
' '.join(describe.split()),
img_link)
img_elem['Height'] = config['image_height']
......
......@@ -19,6 +19,8 @@ Greenland ice sheet elevation contours as seen by the models (CESM's contours
shown in the difference plot).
"""
title = "Net summer longwave radiation"
def make_plot(config, out_path='.'):
# ---------------- Data source in TITAN ------------------------
......@@ -238,7 +240,7 @@ def make_plot(config, out_path='.'):
img_link = os.path.join(os.path.basename(out_path),
os.path.basename(wks_img + '.' + wks_type))
img_elem = el.image('CESM_RACMO23_lwsn',
img_elem = el.image(title,
' '.join(describe.split()),
img_link)
img_elem['Height'] = config['image_height']
......
......@@ -19,6 +19,8 @@ Greenland ice sheet elevation contours as seen by the models (CESM's contours
shown in the difference plot).
"""
title = "Net summer radiation"
def make_plot(config, out_path='.'):
# ---------------- Data source in TITAN ------------------------
......@@ -247,7 +249,7 @@ def make_plot(config, out_path='.'):
img_link = os.path.join(os.path.basename(out_path),
os.path.basename(wks_img + '.' + wks_type))
img_elem = el.image('CESM_RACMO23_rnet',
img_elem = el.image(title,
' '.join(describe.split()),
img_link)
img_elem['Height'] = config['image_height']
......
......@@ -19,6 +19,8 @@ elevation contours as seen by the models (CESM's contours shown in the
difference plot).
"""
title = "Net shortwave radiation"
def make_plot(config, out_path='.'):
# ---------------- Data source in TITAN ------------------------
......@@ -231,7 +233,7 @@ def make_plot(config, out_path='.'):
img_link = os.path.join(os.path.basename(out_path),
os.path.basename(wks_img + '.' + wks_type))
img_elem = el.image('CESM_RACMO23_senf',
img_elem = el.image(title,
' '.join(describe.split()),
img_link)
img_elem['Height'] = config['image_height']
......
......@@ -19,6 +19,8 @@ the 0, 1000, 2000, and 3000 meter Greenland ice sheet elevation contours as seen
the models (CESM's contours shown in the difference plot).
"""
title = "Summer surface mass balance"
def make_plot(config, out_path='.'):
# ---------------- Data source in TITAN ------------------------
......@@ -252,7 +254,7 @@ def make_plot(config, out_path='.'):
img_link = os.path.join(os.path.basename(out_path),
os.path.basename(wks_img + '.' + wks_type))
img_elem = el.image('CESM_RACMO23_smb',
img_elem = el.image(title,
' '.join(describe.split()),
img_link)
img_elem['Height'] = config['image_height']
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment