Merge pull request #1083 from xylar/fix-sic-diff-over-under

Set under/over indices for sea ice conc and thickness

Author: xylar

docs/users_guide/tasks/climatologyMapSeaIceConcNH.rst

@@ -18,42 +18,43 @@ The following configuration options are available for this task::
 
   [climatologyMapSeaIceConcNH]
   ## options related to plotting horizontally remapped climatologies of
-  ## sea ice concentration against reference model results and observations
+  ## sea ice concentration against control model results and observations
   ## in the northern hemisphere (NH)
 
   # colormap for model/observations
   colormapNameResult = ice
+  # 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.15, 0.3, 0.5, 0.7, 0.8, 0.85, 0.9, 0.95, 1]
 
   # 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, 112, 128, 128, 144, 160, 192,
-                               224, 255]
+  colormapIndicesDifference = [0, 0,  26,  51,  77, 102, 128, 128, 153, 179, 204, 230, 255, 255]
   # colormap levels/values for contour boundaries
-  colorbarLevelsDifference = [-1., -0.8, -0.6, -0.4, -0.2, -0.1, 0, 0.1, 0.2,
-                              0.4, 0.6, 0.8, 1.]
+  colorbarLevelsDifference = [-0.5, -0.4, -0.3, -0.2, -0.1, -0.05, 0, 0.05, 0.1, 0.2, 0.3, 0.4, 0.5]
 
   # Months or seasons to plot (These should be left unchanged, since
   # observations are only available for these seasons)
   seasons =  ['JFM', 'JAS']
 
-  # 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
+  # comparison grid(s) (typically 'arctic_extended') on which to plot analysis
+  comparisonGrids = ['arctic_extended']
 
   # a list of prefixes describing the sources of the observations to be used
   observationPrefixes = ['NASATeam', 'Bootstrap']
 
   # arrange subplots vertically?
   vertical = False
 
+  # the minimum threshold below which concentration is masked out
+  minConcentration = 0.15
+
   # observations files
   concentrationNASATeamNH_JFM = SSMI/NASATeam_NSIDC0051/SSMI_NASATeam_gridded_concentration_NH_jfm.interp0.5x0.5_20180710.nc
   concentrationNASATeamNH_JAS = SSMI/NASATeam_NSIDC0051/SSMI_NASATeam_gridded_concentration_NH_jas.interp0.5x0.5_20180710.nc

docs/users_guide/tasks/climatologyMapSeaIceConcSH.rst

@@ -18,42 +18,43 @@ The following configuration options are available for this task::
 
   [climatologyMapSeaIceConcSH]
   ## options related to plotting horizontally remapped climatologies of
-  ## sea ice concentration against reference model results and observations
+  ## sea ice concentration against control model results and observations
   ## in the southern hemisphere (SH)
 
   # colormap for model/observations
   colormapNameResult = ice
+  # 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.15, 0.3, 0.5, 0.7, 0.8, 0.85, 0.9, 0.95, 1]
 
   # 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, 112, 128, 128, 144, 160, 192,
-                               224, 255]
+  colormapIndicesDifference = [0, 0,  26,  51,  77, 102, 128, 128, 153, 179, 204, 230, 255, 255]
   # colormap levels/values for contour boundaries
-  colorbarLevelsDifference = [-1., -0.8, -0.6, -0.4, -0.2, -0.1, 0, 0.1, 0.2,
-                              0.4, 0.6, 0.8, 1.]
+  colorbarLevelsDifference = [-0.5, -0.4, -0.3, -0.2, -0.1, -0.05, 0, 0.05, 0.1, 0.2, 0.3, 0.4, 0.5]
 
   # Months or seasons to plot (These should be left unchanged, since
   # observations are only available for these seasons)
   seasons =  ['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
+  # comparison grid(s) (typically 'antarctic_extended') on which to plot analysis
+  comparisonGrids = ['antarctic_extended']
 
   # a list of prefixes describing the sources of the observations to be used
   observationPrefixes = ['NASATeam', 'Bootstrap']
 
   # arrange subplots vertically?
   vertical = False
 
+  # the minimum threshold below which concentration is masked out
+  minConcentration = 0.15
+
   # observations files
   concentrationNASATeamSH_DJF = SSMI/NASATeam_NSIDC0051/SSMI_NASATeam_gridded_concentration_SH_djf.interp0.5x0.5_20180710.nc
   concentrationNASATeamSH_JJA = SSMI/NASATeam_NSIDC0051/SSMI_NASATeam_gridded_concentration_SH_jja.interp0.5x0.5_20180710.nc

docs/users_guide/tasks/climatologyMapSeaIceThickNH.rst

@@ -19,33 +19,33 @@ The following configuration options are available for this task::
 
   [climatologyMapSeaIceThickNH]
   ## options related to plotting horizontally remapped climatologies of
-  ## sea ice thickness against reference model results and observations
+  ## sea ice thickness against control model results and observations
   ## in the northern hemisphere (NH)
 
   # colormap for model/observations
-  colormapNameResult = ice
+  colormapNameResult = davos
+  # 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.25, 0.5, 1, 1.5, 2, 2.5, 3, 3.5]
 
   # 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]
+  colormapIndicesDifference = [0, 0, 32, 64, 96, 128, 128, 160, 192, 223, 255, 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.]
 
   # Months or seasons to plot (These should be left unchanged, since
   # observations are only available for these seasons)
   seasons =  ['FM', 'ON']
 
-  # 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
+  # comparison grid(s) (typically 'arctic_extended') on which to plot analysis
+  comparisonGrids = ['arctic_extended']
 
   # a list of prefixes describing the sources of the observations to be used
   observationPrefixes = ['']

docs/users_guide/tasks/climatologyMapSeaIceThickSH.rst

@@ -19,33 +19,33 @@ The following configuration options are available for this task::
 
   [climatologyMapSeaIceThickSH]
   ## options related to plotting horizontally remapped climatologies of
-  ## sea ice thickness against reference model results and observations
+  ## sea ice thickness against control model results and observations
   ## in the southern hemisphere (SH)
 
   # colormap for model/observations
-  colormapNameResult = ice
+  colormapNameResult = davos
+  # 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.2, 0.4, 0.6, 0.8, 1, 1.5, 2, 2.5]
 
   # 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]
+  colormapIndicesDifference = [0, 0, 32, 64, 96, 128, 128, 160, 192, 223, 255, 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.]
 
   # Months or seasons to plot (These should be left unchanged, since
   # observations are only available for these seasons)
   seasons =  ['FM', 'ON']
 
-  # 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
+  # comparison grid(s) (typically 'antarctic_extended') on which to plot analysis
+  comparisonGrids = ['antarctic_extended']
 
   # a list of prefixes describing the sources of the observations to be used
   observationPrefixes = ['']

mpas_analysis/default.cfg

@@ -4711,7 +4711,7 @@ colormapNameDifference = balance
 # whether the colormap is indexed or continuous
 colormapTypeDifference = indexed
 # color indices into colormapName for filled contours
-colormapIndicesDifference = [0, 32, 64, 96, 112, 128, 128, 144, 160, 192, 224, 255]
+colormapIndicesDifference = [0, 0,  26,  51,  77, 102, 128, 128, 153, 179, 204, 230, 255, 255]
 # colormap levels/values for contour boundaries
 colorbarLevelsDifference = [-0.5, -0.4, -0.3, -0.2, -0.1, -0.05, 0, 0.05, 0.1, 0.2, 0.3, 0.4, 0.5]
 
@@ -4757,7 +4757,7 @@ colormapNameDifference = balance
 # whether the colormap is indexed or continuous
 colormapTypeDifference = indexed
 # color indices into colormapName for filled contours
-colormapIndicesDifference = [0, 32, 64, 96, 112, 128, 128, 144, 160, 192, 224, 255]
+colormapIndicesDifference = [0, 0,  26,  51,  77, 102, 128, 128, 153, 179, 204, 230, 255, 255]
 # colormap levels/values for contour boundaries
 colorbarLevelsDifference = [-0.5, -0.4, -0.3, -0.2, -0.1, -0.05, 0, 0.05, 0.1, 0.2, 0.3, 0.4, 0.5]
 
@@ -4803,7 +4803,7 @@ 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]
+colormapIndicesDifference = [0, 0, 32, 64, 96, 128, 128, 160, 192, 223, 255, 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.]
 
@@ -4844,7 +4844,7 @@ 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]
+colormapIndicesDifference = [0, 0, 32, 64, 96, 128, 128, 160, 192, 223, 255, 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.]
 

mpas_analysis/sea_ice/climatology_map_sea_ice_conc.py

@@ -184,7 +184,7 @@ def _add_obs_tasks(self, seasons, comparisonGridNames, hemisphere,
                         galleryName='Observations: SSM/I {}'.format(
                             prefix),
                         maskMinThreshold=minConcentration,
-                        extend='neither',
+                        extend='both',
                         prependComparisonGrid=False)
 
                     self.add_subtask(subtask)
@@ -231,7 +231,7 @@ def _add_ref_tasks(self, seasons, comparisonGridNames, hemisphere,
                     groupLink='{}_conc'.format(hemisphere.lower()),
                     galleryName=galleryName,
                     maskMinThreshold=minConcentration,
-                    extend='neither',
+                    extend='both',
                     prependComparisonGrid=False)
 
                 self.add_subtask(subtask)

mpas_analysis/sea_ice/climatology_map_sea_ice_thick.py

@@ -166,7 +166,7 @@ def __init__(self, config, mpasClimatologyTask, hemisphere,
                     groupLink=f'{hemisphere.lower()}_thick',
                     galleryName=galleryName,
                     maskMinThreshold=0,
-                    extend='neither',
+                    extend='both',
                     prependComparisonGrid=False)
 
                 self.add_subtask(subtask)