Add figure captions and references to the dynamics extension

For #20

dynamics/dynamics/cisma_js_standalone.py

+# -*- coding: utf-8 -*-
+
 import os
 import Ngl
 
@@ -7,7 +9,12 @@ from netCDF4 import Dataset
 from livvkit.util import elements as el
 
 
-describe = """CISMA_js plot"""
+describe = """
+The norm of the surface velocity (m s^-1) from CISM initialization dataset 
+(left; Joughin et al., 2010a,b), CISM-Albany after spin-up (middle), and their
+difference (CISM-A - init) over Jakobshavn glacier (West Greenland; Zwally et al. 
+(2012) basins 7).   
+"""
 # plots velnorm, zooms in on Jakobshavn 69.2N 51.1W
 
 

dynamics/dynamics/cisma_pm_standalone.py

+# -*- coding: utf-8 -*-
+
 import os
 import Ngl
 
@@ -7,7 +9,12 @@ from netCDF4 import Dataset
 from livvkit.util import elements as el
 
 
-describe = """CISMA_pm plot"""
+describe = """
+The norm of the surface velocity (m s^-1) from CISM initialization dataset 
+(left; Joughin et al., 2010a,b), CISM-Albany after spin-up (middle), and their
+difference (CISM-A - init) over Petermann glacier (North-West Greenland; Zwally et al. 
+(2012) basins 1).   
+"""
 # Plots velnorm, zooms in on Petermann glacier, 80.5N -59.5
 
 

dynamics/dynamics/cisma_racmo23_scatter_standalone.py

+# -*- coding: utf-8 -*-
+
 import os
 
 import numpy as np
@@ -8,7 +10,11 @@ from netCDF4 import Dataset
 from livvkit.util import elements as el
 
 
-describe = """scatter_CISMA_RACMO23 plot."""
+describe = """
+Comparison of the surface mass balance (kg m^-2 a^-1) at grid point over Greenland 
+for RACMO 2.3 and CISM-Albany (after spin-up forced by RACMO 2.0) interpolated to 
+the coarser RACMO grid.
+"""
 
 
 def make_plot(config, out_path='.'):

dynamics/dynamics/cisma_racmo23_standalone.py

+# -*- coding: utf-8 -*-
+
 import os
 import Ngl
 
@@ -7,7 +9,11 @@ from netCDF4 import Dataset
 from livvkit.util import elements as el
 
 
-describe = """CISMA_RACMO23 plot"""
+describe = """
+Surface mass balance (kg m^-2 a^-1) over Greenland for CISM-Albany after spin-up 
+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.
+"""
 
 
 def make_plot(config, out_path='.'):

dynamics/dynamics/cisma_temp_standalone.py

+# -*- coding: utf-8 -*-
+
 import os
 import Ngl
 
@@ -7,7 +9,7 @@ from netCDF4 import Dataset
 from livvkit.util import elements as el
 
 
-describe = """CISMA_temp plot"""
+describe = """Surface temperature (°C) of the Greenland ice-sheet after CISM-Albany spin-up."""
 
 
 def make_plot(config, out_path='.'):

dynamics/dynamics/cisma_velnorm_standalone.py

+# -*- coding: utf-8 -*-
+
 import os
 import Ngl
 
@@ -8,7 +10,11 @@ from netCDF4 import Dataset
 from livvkit.util import elements as el
 
 
-describe = """CISMA_velnorm plot"""
+describe = """
+The norm of the surface velocity (m s^-1) from CISM initialization dataset 
+(left; Joughin et al., 2010a,b), CISM-Albany after spin-up (middle), and their
+difference (CISM-A - init) over Greenland.   
+"""
 
 
 def make_plot(config, out_path='.'):

dynamics/dynamics/cisma_zach_standalone.py

+# -*- coding: utf-8 -*-
+
 import os
 import Ngl
 
@@ -8,7 +10,11 @@ from netCDF4 import Dataset
 from livvkit.util import elements as el
 
 
-describe = """CISMA_zach plot"""
+describe = """
+The norm of the surface velocity (m s^-1) from CISM initialization dataset 
+(left; Joughin et al., 2010a,b), CISM-Albany after spin-up (middle), and their
+difference (CISM-A - init) over North-East Greenland (Zwally et al. (2012) basins 1 and 2).   
+"""
 # zoomed in velnorm over Zachariae glacier at 78N, -30E
 
 

dynamics/dynamics/utils.py

+# -*- coding: utf-8 -*-
+
+from __future__ import unicode_literals
+
 import pybtex.database
 import pybtex.io
 

dynamics/dynamics_cisma.bib

@@ -9,3 +9,53 @@
     PAGES = {854--869},
     YEAR = {2017},
 }
+
+
+@misc{Zwally2012,
+    AUTHOR = {H. Jay Zwally and Mario B. Giovinetto and Matthew A. Beckley and Jack L. Saba},
+    TITLE = {Antarctic and Greenland Drainage Systems},
+    NOTE = {GSFC Cryospheric Sciences Laboratory},
+    HOWPUBLISHED = {\url{http://icesat4.gsfc.nasa.gov/cryo_data/ant_grn_drainage_systems.php}},
+    YEAR = {2012}
+}
+
+@techreport{insar2010b,
+    AUTHOR = {Joughin, I. and B. Smith and I. Howat and T. Scambos},
+    TITLE = {{MEaSUREs Greenland Ice Velocity Map from InSAR Data}},
+    INSTITUTION = {National Snow and Ice Data Center},
+    ADDRESS = {Boulder, Colorado USA},
+    URL = {http://nsidc.org/data/docs/measures/nsidc0478_joughin/},
+    YEAR = {2010},
+    MONTH = {April},
+    DOI = {10.5067/MEASURES/CRYOSPHERE/nsidc-0478.001}
+}
+
+@Article{insar2010a,
+    AUTHOR = {Joughin, I. and B. Smith and I. M. Howat and T. Scambos and T. Moon},
+    TITLE = {Greenland Flow Variability from Ice-Sheet-Wide Velocity Mapping},
+    JOURNAL = {Journal of Glaciology},
+    VOLUME = {56},
+    YEAR = {2010},
+    PAGES = {415--430},
+    DOI = {10.3189/002214310792447734}
+}
+
+@article{racmo20,
+    AUTHOR = {van Angelen, J. H. and M.R. Van Den Broeke and B. Wouters and J. T. M. Lenaerts},
+    TITLE = {Contemporary (1960–2012) Evolution of the Climate and Surface Mass Balance of the Greenland Ice Sheet},
+    JOURNAL = {Surv. Geophys.},
+    DOI = {10.1007/s10712-013-9261-z},
+    VOLUME = {35},
+    YEAR = {2013},
+    PAGES = {1155--1174}
+ }
+
+@article{noel15 ,
+	AUTHOR = {B. N{\"o}el and W. J. van de Berg and E. van Meijgaard and P. Kuipers Munneke and R. S. W. van de Wal and M. R. van den Broeke},
+	TITLE = {Evaluation of the updated regional climate model RACMO2.3: summer snowfall impact on the Greenland ice sheet},
+	JOURNAL = {The Cryosphere},
+	YEAR = {2015},
+	DOI = {10.5194/tc-9-1831-2015},
+	VOLUME = {9},
+	PAGES = {1831--1844}
+}

dynamics/dynamics_cisma.py

+# -*- coding: utf-8 -*-
+#
 # Copyright (c) 2015,2016, UT-BATTELLE, LLC
 # All rights reserved.
 #
@@ -39,11 +41,12 @@ from livvkit.util import functions as fn
 with fn.temp_sys_path(os.path.dirname(__file__)):
     import dynamics.utils as utils
 
+    import dynamics.cisma_velnorm_standalone as cisma_velnorm
     import dynamics.cisma_js_standalone as cisma_jakobshavn
     import dynamics.cisma_pm_standalone as cisma_petermann
     import dynamics.cisma_zach_standalone as cisma_zachariae
-    import dynamics.cisma_temp_standalone as cisma_temp
-    import dynamics.cisma_velnorm_standalone as cisma_velnorm
+    # FIXME not clear what purpose this serves
+    # import dynamics.cisma_temp_standalone as cisma_temp
     import dynamics.cisma_racmo23_standalone as cisma_racmo_contour
     import dynamics.cisma_racmo23_scatter_standalone as cisma_racmo_scatter
 
@@ -65,11 +68,11 @@ def run(name, config):
 
     # PLOTS
     img_list = []
+    img_list.extend(cisma_velnorm.make_plot(config, out_path=img_dir))
     img_list.extend(cisma_jakobshavn.make_plot(config, out_path=img_dir))
     img_list.extend(cisma_petermann.make_plot(config, out_path=img_dir))
     img_list.extend(cisma_zachariae.make_plot(config, out_path=img_dir))
-    img_list.extend(cisma_temp.make_plot(config, out_path=img_dir))
-    img_list.extend(cisma_velnorm.make_plot(config, out_path=img_dir))
+    # img_list.extend(cisma_temp.make_plot(config, out_path=img_dir))
     img_list.extend(cisma_racmo_contour.make_plot(config, out_path=img_dir))
     img_list.extend(cisma_racmo_scatter.make_plot(config, out_path=img_dir))