Add sea ice production/melting analysis tasks

mpas_analysis/__main__.py

@@ -238,6 +238,18 @@ def build_analysis_list(config, controlConfig):
         config=config, mpasClimatologyTask=seaIceClimatolgyTask,
         hemisphere='SH', controlConfig=controlConfig))
     analyses.append(seaIceTimeSeriesTask)
+    analyses.append(sea_ice.ClimatologyMapSeaIceProduction(
+        config=config, mpas_climatology_task=seaIceClimatolgyTask,
+        hemisphere='NH', control_config=controlConfig))
+    analyses.append(sea_ice.ClimatologyMapSeaIceProduction(
+        config=config, mpas_climatology_task=seaIceClimatolgyTask,
+        hemisphere='SH', control_config=controlConfig))
+    analyses.append(sea_ice.ClimatologyMapSeaIceMelting(
+        config=config, mpas_climatology_task=seaIceClimatolgyTask,
+        hemisphere='NH', control_config=controlConfig))
+    analyses.append(sea_ice.ClimatologyMapSeaIceMelting(
+        config=config, mpas_climatology_task=seaIceClimatolgyTask,
+        hemisphere='SH', control_config=controlConfig))
 
     analyses.append(sea_ice.TimeSeriesSeaIce(config, seaIceTimeSeriesTask,
                                              controlConfig))

mpas_analysis/default.cfg

@@ -3736,3 +3736,157 @@ areaNH = IceArea_timeseries/iceAreaNH_climo_20180710.nc
 areaSH = IceArea_timeseries/iceAreaSH_climo_20180710.nc
 volNH = PIOMAS/PIOMASvolume_monthly_climo_20180710.nc
 volSH = none
+
+
+[climatologyMapSeaIceProductionNH]
+# options related to plotting horizontally remapped climatologies of
+# sea ice production against control model results and observations
+# in the northern hemisphere (NH)
+
+# colormap for model/observations
+colormapNameResult = dense
+# whether the colormap is indexed or continuous
+colormapTypeResult = indexed
+# color indices into colormapName for filled contours
+colormapIndicesResult = [20, 80, 110, 140, 170, 200, 230, 255]
+# colormap levels/values for contour boundaries
+colorbarLevelsResult = [0, 0.1, 0.25, 0.5, 1., 2., 3., 5., 10.]
+
+# colormap for differences
+colormapNameDifference = balance
+# whether the colormap is indexed or continuous
+colormapTypeDifference = indexed
+# color indices into colormapName for filled contours
+colormapIndicesDifference = [0, 32, 64, 96, 128, 128, 160, 192, 224, 255]
+# colormap levels/values for contour boundaries
+colorbarLevelsDifference = [-3., -2.5, -2, -0.5, -0.1, 0, 0.1, 0.5, 2, 2.5, 3.]
+
+# Times for comparison times
+seasons =  ['ANN', 'DJF', 'JJA']
+
+# comparison grid(s) ('latlon', 'antarctic') on which to plot analysis
+comparisonGrids = ['latlon']
+
+# reference lat/lon for sea ice plots in the northern hemisphere
+minimumLatitude = 50
+referenceLongitude = 0
+
+# arrange subplots vertically?
+vertical = False
+
+
+[climatologyMapSeaIceProductionSH]
+# options related to plotting horizontally remapped climatologies of
+# sea ice production against control model results and observations
+# in the southern hemisphere (SH)
+
+# colormap for model/observations
+colormapNameResult = dense
+# whether the colormap is indexed or continuous
+colormapTypeResult = indexed
+# color indices into colormapName for filled contours
+colormapIndicesResult = [20, 80, 110, 140, 170, 200, 230, 255]
+# colormap levels/values for contour boundaries
+colorbarLevelsResult = [0, 0.1, 0.25, 0.5, 1., 2., 3., 5., 10.]
+
+# colormap for differences
+colormapNameDifference = balance
+# whether the colormap is indexed or continuous
+colormapTypeDifference = indexed
+# color indices into colormapName for filled contours
+colormapIndicesDifference = [0, 32, 64, 96, 128, 128, 160, 192, 224, 255]
+# colormap levels/values for contour boundaries
+colorbarLevelsDifference = [-3., -2.5, -2, -0.5, -0.1, 0, 0.1, 0.5, 2, 2.5, 3.]
+
+# Times for comparison times
+seasons =  ['ANN', 'DJF', 'JJA']
+
+# comparison grid(s) ('latlon', 'antarctic') on which to plot analysis
+comparisonGrids = ['latlon']
+
+# reference lat/lon for sea ice plots in the northern hemisphere
+minimumLatitude = -50
+referenceLongitude = 180
+
+# arrange subplots vertically?
+vertical = False
+
+# observations files
+productionSH = Haumann/aniceflux_v01_clim_mean_1982-2008.nc
+
+
+[climatologyMapSeaIceMeltingNH]
+# options related to plotting horizontally remapped climatologies of
+# sea ice melting against control model results and observations
+# in the northern hemisphere (NH)
+
+# colormap for model/observations
+colormapNameResult = amp
+# whether the colormap is indexed or continuous
+colormapTypeResult = indexed
+# color indices into colormapName for filled contours
+colormapIndicesResult = [20, 80, 110, 140, 170, 200, 230, 255]
+# colormap levels/values for contour boundaries
+colorbarLevelsResult = [0, 0.1, 0.25, 0.5, 1., 2., 3., 5., 10.]
+
+# colormap for differences
+colormapNameDifference = balance
+# whether the colormap is indexed or continuous
+colormapTypeDifference = indexed
+# color indices into colormapName for filled contours
+colormapIndicesDifference = [0, 32, 64, 96, 128, 128, 160, 192, 224, 255]
+# colormap levels/values for contour boundaries
+colorbarLevelsDifference = [-3., -2.5, -2, -0.5, -0.1, 0, 0.1, 0.5, 2, 2.5, 3.]
+
+# Times for comparison times
+seasons =  ['ANN', 'DJF', 'JJA']
+
+# comparison grid(s) ('latlon', 'antarctic') on which to plot analysis
+comparisonGrids = ['latlon']
+
+# reference lat/lon for sea ice plots in the northern hemisphere
+minimumLatitude = 50
+referenceLongitude = 0
+
+# arrange subplots vertically?
+vertical = False
+
+
+[climatologyMapSeaIceMeltingSH]
+# options related to plotting horizontally remapped climatologies of
+# sea ice melting against control model results and observations
+# in the southern hemisphere (SH)
+
+# colormap for model/observations
+colormapNameResult = amp
+# whether the colormap is indexed or continuous
+colormapTypeResult = indexed
+# color indices into colormapName for filled contours
+colormapIndicesResult = [20, 80, 110, 140, 170, 200, 230, 255]
+# colormap levels/values for contour boundaries
+colorbarLevelsResult = [0, 0.1, 0.25, 0.5, 1., 2., 3., 5., 10.]
+
+# colormap for differences
+colormapNameDifference = balance
+# whether the colormap is indexed or continuous
+colormapTypeDifference = indexed
+# color indices into colormapName for filled contours
+colormapIndicesDifference = [0, 32, 64, 96, 128, 128, 160, 192, 224, 255]
+# colormap levels/values for contour boundaries
+colorbarLevelsDifference = [-3., -2.5, -2, -0.5, -0.1, 0, 0.1, 0.5, 2, 2.5, 3.]
+
+# Times for comparison times
+seasons =  ['ANN', 'DJF', 'JJA']
+
+# comparison grid(s) ('latlon', 'antarctic') on which to plot analysis
+comparisonGrids = ['latlon']
+
+# reference lat/lon for sea ice plots in the northern hemisphere
+minimumLatitude = -50
+referenceLongitude = 180
+
+# arrange subplots vertically?
+vertical = False
+
+# observations files
+meltingSH = Haumann/aniceflux_v01_clim_mean_1982-2008.nc

mpas_analysis/sea_ice/__init__.py

@@ -5,3 +5,7 @@
 from mpas_analysis.sea_ice.time_series import TimeSeriesSeaIce
 from mpas_analysis.sea_ice.climatology_map_berg_conc import \
     ClimatologyMapIcebergConc
+from mpas_analysis.sea_ice.climatology_map_production import \
+    ClimatologyMapSeaIceProduction
+from mpas_analysis.sea_ice.climatology_map_melting import \
+    ClimatologyMapSeaIceMelting