Add moc script (postprocessing) and plot_vertical_section

config.default

@@ -61,11 +61,12 @@ logsSubdirectory = logs
 mpasClimatologySubdirectory = clim/mpas
 mpasRegriddedClimSubdirectory = clim/mpas/regridded
 mappingSubdirectory = mapping
+timeSeriesSubdirectory = timeseries
 
 # a list of analyses to generate.  Valid names are:
 #   'timeSeriesOHC', 'timeSeriesSST', 'regriddedSST',
-#   'regriddedSSS', 'regriddedMLD', 'timeSeriesSeaIceAreaVol',
-#   'regriddedSeaIceConcThick'
+#   'regriddedSSS', 'regriddedMLD', 'streamfunctionMOC',
+#   'timeSeriesSeaIceAreaVol', 'regriddedSeaIceConcThick'
 # the following shortcuts exist:
 #   'all' -- all analyses will be run
 #   'all_timeSeries' -- all time-series analyses will be run
@@ -85,7 +86,8 @@ generate = ['all']
 
 # alternative examples that would perform all analysis except
 #   'timeSeriesOHC'
-#generate = ['timeSeriesSST', 'all_regriddedHorizontal', 'all_seaIce']
+#generate = ['timeSeriesSST', 'streamfunctionMOC',
+#            'all_regriddedHorizontal', 'all_seaIce']
 #generate = ['all', 'no_timeSeriesOHC']
 # Each subsequent list entry can be used to alter previous list entries. For
 # example, the following would produce all analyses except regriddedSST,
@@ -304,19 +306,47 @@ regionIndicesToPlot = [6]
 # window)
 movingAveragePoints = 12
 
-[timeSeriesMOC]
-## options related to plotting time series of meridional overturning
+[streamfunctionMOC]
+## options related to plotting the streamfunction of the meridional overturning 
 ## circulation (MOC)
 
+# Region names for basin MOC calculation.
+# Supported options are Atlantic and IndoPacific
+regionNames = ['Atlantic']
+
+# Mask file for post-processing regional MOC computation
+regionMaskFiles = /path/to/MOCregional/mapping/file
+
+# Size of latitude bins over which MOC streamfunction is integrated
+latBinSizeGlobal = 1.
+latBinSizeAtlantic = 0.5
+latBinSizeIndoPacific = 0.5
+
+# colormap for model results
+colormapNameGlobal = RdYlBu_r
+colormapNameAtlantic = RdYlBu_r
+colormapNameIndoPacific = RdYlBu_r
+# colormap indices for contour color
+colormapIndicesGlobal = [0, 40, 80, 110, 140, 170, 200, 230, 255]
+colormapIndicesAtlantic = [0, 40, 80, 110, 140, 170, 200, 230, 255]
+colormapIndicesIndoPacific = [0, 40, 80, 110, 140, 170, 200, 230, 255]
+# colorbar levels/values for contour boundaries
+colorbarLevelsGlobal = [-20, -10, -5, -2, 2, 5, 10, 20, 30, 40]
+colorbarLevelsAtlantic = [-10, -5, -2, 0, 5, 8, 10, 14, 18, 22]
+colorbarLevelsIndoPacific = [-10, -5, -2, 0, 5, 8, 10, 14, 18, 22]
+# contour line levels
+contourLevelsGlobal = np.arange(-25.1,35.1,10)
+contourLevelsAtlantic = np.arange(-8,20.1,2)
+contourLevelsIndoPacific = np.arange(-8,20.1,2)
+
 ## compare to output from another model run?
 #compareWithModel = True
 # compare to observations?
 compareWithObservations = True
-# list of region indices to plot from the region list in [regions] below
-regionIndicesToPlot = [6]
-# Number of points over which to compute moving average (e.g., for monthly
-# output, movingAveragePoints=12 corresponds to a 12-month moving average
-# window)
+
+# Number of points over which to compute moving average for
+# MOC timeseries (e.g., for monthly output, movingAveragePoints=12
+# corresponds to a 12-month moving average window)
 movingAveragePoints = 12
 
 [timeSeriesSeaIceAreaVol]
@@ -342,18 +372,18 @@ polarPlot = False
 ## (SST) against reference model results and observations
 
 # colormap for model/observations
-resultColormap = RdYlBu_r
-# indices into resultColormap for contour color
-resultColormapIndices = [0, 40, 80, 110, 140, 170, 200, 230, 255]
+colormapNameResult = RdYlBu_r
+# color indices into colormapName for filled contours
+colormapIndicesResult = [0, 40, 80, 110, 140, 170, 200, 230, 255]
 # colormap levels/values for contour boundaries
-resultContourValues = [-2, 0, 2, 6, 10, 16, 22, 26, 28, 32]
+colorbarLevelsResult = [-2, 0, 2, 6, 10, 16, 22, 26, 28, 32]
 
 # colormap for differences
-differenceColormap = RdBu_r
-# indices into differenceColormap for contour color
-differenceColormapIndices = [0, 28, 57, 85, 113, 142, 170, 198, 227, 255]
+colormapNameDifference = RdBu_r
+# color indices into colormapName for filled contours
+colormapIndicesDifference = [0, 28, 57, 85, 113, 142, 170, 198, 227, 255]
 # colormap levels/values for contour boundaries
-differenceContourValues = [-5, -3, -2, -1, 0, 1, 2, 3, 5]
+colorbarLevelsDifference = [-5, -3, -2, -1, 0, 1, 2, 3, 5]
 
 # Times for comparison times (Jan, Feb, Mar, Apr, May, Jun, Jul, Aug, Sep, Oct, Nov, Dec, JFM, AMJ, JAS, OND, ANN)
 comparisonTimes =  ['JFM', 'JAS', 'ANN']
@@ -363,18 +393,18 @@ comparisonTimes =  ['JFM', 'JAS', 'ANN']
 ## (SSS) against reference model results and observations
 
 # colormap for model/observations
-resultColormap = RdYlBu_r
-# indices into resultColormap for contour color
-resultColormapIndices = [0, 40, 80, 110, 140, 170, 200, 230, 255]
+colormapNameResult = RdYlBu_r
+# color indices into colormapName for filled contours
+colormapIndicesResult = [0, 40, 80, 110, 140, 170, 200, 230, 255]
 # colormap levels/values for contour boundaries
-resultContourValues = [28, 29, 30, 31, 32, 33, 34, 35, 36, 38]
+colorbarLevelsResult = [28, 29, 30, 31, 32, 33, 34, 35, 36, 38]
 
 # colormap for differences
-differenceColormap = RdBu_r
-# indices into differenceColormap for contour color
-differenceColormapIndices = [0, 28, 57, 85, 113, 142, 170, 198, 227, 255]
+colormapNameDifference = RdBu_r
+# color indices into colormapName for filled contours
+colormapIndicesDifference = [0, 28, 57, 85, 113, 142, 170, 198, 227, 255]
 # colormap levels/values for contour boundaries
-differenceContourValues = [-3, -2, -1, -0.5, 0, 0.5, 1, 2, 3]
+colorbarLevelsDifference = [-3, -2, -1, -0.5, 0, 0.5, 1, 2, 3]
 
 # Times for comparison times (Jan, Feb, Mar, Apr, May, Jun, Jul, Aug, Sep, Oct, Nov, Dec, JFM, AMJ, JAS, OND, ANN)
 comparisonTimes =  ['JFM', 'JAS', 'ANN']
@@ -384,18 +414,18 @@ comparisonTimes =  ['JFM', 'JAS', 'ANN']
 ## (MLD) against reference model results and observations
 
 # colormap for model/observations
-resultColormap = viridis
-# indices into resultColormap for contour color
-resultColormapIndices = [0, 40, 80, 110, 140, 170, 200, 230, 255]
+colormapNameResult = viridis
+# color indices into colormapName for filled contours
+colormapIndicesResult = [0, 40, 80, 110, 140, 170, 200, 230, 255]
 # colormap levels/values for contour boundaries
-resultContourValues = [0, 20, 40, 60, 80, 100, 150, 200, 400, 800]
+colorbarLevelsResult = [0, 20, 40, 60, 80, 100, 150, 200, 400, 800]
 
 # colormap for differences
-differenceColormap = RdBu_r
-# indices into differenceColormap for contour color
-differenceColormapIndices = [0, 28, 57, 85, 113, 142, 170, 198, 227, 255]
+colormapNameDifference = RdBu_r
+# color indices into colormapName for filled contours
+colormapIndicesDifference = [0, 28, 57, 85, 113, 142, 170, 198, 227, 255]
 # colormap levels/values for contour boundaries
-differenceContourValues = [-150, -80, -30, -10, 0, 10, 30, 80, 150]
+colorbarLevelsDifference = [-150, -80, -30, -10, 0, 10, 30, 80, 150]
 
 # Times for comparison times (Jan, Feb, Mar, Apr, May, Jun, Jul, Aug, Sep, Oct, Nov, Dec, JFM, AMJ, JAS, OND, ANN)
 comparisonTimes =  ['JFM', 'JAS', 'ANN']
@@ -405,28 +435,40 @@ comparisonTimes =  ['JFM', 'JAS', 'ANN']
 ## and thickness against reference model results and observations
 
 # colormap for model/observations
-resultColormap = inferno
-# indices into resultColormap for contour color
-resultColormapIndices = [20, 80, 110, 140, 170, 200, 230, 255]
+colormapNameConcResultWinter = inferno
+colormapNameConcResultSummer = inferno
+colormapNameThickResultNH = inferno
+colormapNameThickResultSH = inferno
+# color indices into colormapName for filled contours
+colormapIndicesConcResultWinter = [20, 80, 110, 140, 170, 200, 230, 255]
+colormapIndicesConcResultSummer = [20, 80, 110, 140, 170, 200, 230, 255]
+colormapIndicesThickResultNH = [20, 80, 110, 140, 170, 200, 230, 255]
+colormapIndicesThickResultSH = [20, 80, 110, 140, 170, 200, 230, 255]
 # colormap levels/values for contour boundaries for:
 # concentration in winter and summer
-resultConcWinterContourValues = [0.15, 0.4, 0.7, 0.9, 0.94, 0.96, 0.98, 0.99, 1]
-resultConcSummerContourValues = [0.15, 0.3, 0.5, 0.7, 0.8, 0.85, 0.9, 0.95, 1]
+colorbarLevelsConcResultWinter = [0.15, 0.4, 0.7, 0.9, 0.94, 0.96, 0.98, 0.99, 1]
+colorbarLevelsConcResultSummer = [0.15, 0.3, 0.5, 0.7, 0.8, 0.85, 0.9, 0.95, 1]
 # thickness in the northern and southern hemispheres
-resultThickNHContourValues = [0, 0.25, 0.5, 1, 1.5, 2, 2.5, 3, 3.5]
-resultThickSHContourValues = [0, 0.2, 0.4, 0.6, 0.8, 1, 1.5, 2, 2.5]
+colorbarLevelsThickResultNH = [0, 0.25, 0.5, 1, 1.5, 2, 2.5, 3, 3.5]
+colorbarLevelsThickResultSH = [0, 0.2, 0.4, 0.6, 0.8, 1, 1.5, 2, 2.5]
 
 # colormap for differences
-differenceColormap = RdBu_r
-# indices into differenceColormap for contour color
-differenceColormapIndices = [0, 40, 80, 127, 127, 170, 210, 255]
+colormapNameConcDifferenceWinter = RdBu_r
+colormapNameConcDifferenceSummer = RdBu_r
+colormapNameThickDifferenceNH = RdBu_r
+colormapNameThickDifferenceSH = RdBu_r
+# color indices into colormapName for filled contours
+colormapIndicesConcDifferenceWinter = [0, 40, 80, 127, 127, 170, 210, 255]
+colormapIndicesConcDifferenceSummer = [0, 40, 80, 127, 127, 170, 210, 255]
+colormapIndicesThickDifferenceNH = [0, 40, 80, 127, 127, 170, 210, 255]
+colormapIndicesThickDifferenceSH = [0, 40, 80, 127, 127, 170, 210, 255]
 # colormap levels/values for contour boundaries for:
 # concentration in winter and summer
-differenceConcWinterContourValues = [-0.8, -0.6, -0.4, -0.2, 0, 0.2, 0.4, 0.6, 0.8]
-differenceConcSummerContourValues = [-0.8, -0.6, -0.4, -0.2, 0, 0.2, 0.4, 0.6, 0.8]
+colorbarLevelsConcDifferenceWinter = [-0.8, -0.6, -0.4, -0.2, 0, 0.2, 0.4, 0.6, 0.8]
+colorbarLevelsConcDifferenceSummer = [-0.8, -0.6, -0.4, -0.2, 0, 0.2, 0.4, 0.6, 0.8]
 # thickness in the northern and southern hemispheres
-differenceThickNHContourValues = [-2.5, -2, -0.5, -0.1, 0, 0.1, 0.5, 2, 2.5]
-differenceThickSHContourValues = [-2.5, -2, -0.5, -0.1, 0, 0.1, 0.5, 2, 2.5]
+colorbarLevelsThickDifferenceNH = [-2.5, -2, -0.5, -0.1, 0, 0.1, 0.5, 2, 2.5]
+colorbarLevelsThickDifferenceSH = [-2.5, -2, -0.5, -0.1, 0, 0.1, 0.5, 2, 2.5]
 
 # reference lat/lon for sea ice plots in the northern hemisphere
 minimumLatitudeNH = 50

configs/edison/config.20161006bugfix.alpha8.A_WCYCL1850S.ne30_oEC_ICG.edison

@@ -4,10 +4,6 @@
 
 # mainRunName is a name that identifies the simulation being analyzed.
 mainRunName = 20161006bugfix.alpha8.A_WCYCL1850S.ne30_oEC_ICG.edison
-# referenceRunName is the name of a reference run to compare against (or None
-# to turn off comparison with a reference, e.g. if no reference case is
-# available)
-referenceRunName = None
 # preprocessedReferenceRunName is the name of a reference run that has been
 # preprocessed to compare against (or None to turn off comparison).  Reference
 # runs of this type would have preprocessed results because they were not
@@ -20,11 +16,6 @@ preprocessedReferenceRunName = B1850C5_ne30_v0.4
 
 # directory containing model results
 baseDirectory = /global/cscratch1/sd/jonbob/ACME_simulations/20161006bugfix.alpha8.A_WCYCL1850S.ne30_oEC_ICG.edison/run
-# names of namelist and streams files.  If not in baseDirectory, give full path
-oceanNamelistFileName = mpas-o_in
-oceanStreamsFileName = streams.ocean
-seaIceNamelistFileName = mpas-cice_in
-seaIceStreamsFileName = streams.cice
 
 # names of ocean and sea ice meshes (e.g. EC60to30, QU240, RRS30to10, etc.)
 mpasMeshName = EC60to30
@@ -121,3 +112,15 @@ baseDirectory = /global/project/projectdirs/acme/observations/SeaIce
 # directory where ocean reference simulation results are stored
 baseDirectory =  /global/project/projectdirs/acme/ACMEv0_lowres/B1850C5_ne30_v0.4/ice/postprocessing
 
+[timeSeriesSeaIceAreaVol]
+## options related to plotting time series of sea ice area and volume
+
+# plot on polar plot
+polarPlot = False
+
+[streamfunctionMOC]
+## options related to plotting the streamfunction of the meridional overturning 
+## circulation (MOC)
+
+# Mask file for ocean basin regional computation
+regionMaskFiles = /global/project/projectdirs/acme/mapping/grids/EC60to30v1_SingleRegionAtlanticWTransportTransects_masks.nc

configs/edison/config.20161117.beta0.A_WCYCL1850.ne30_oEC.edison

@@ -4,10 +4,6 @@
 
 # mainRunName is a name that identifies the simulation being analyzed.
 mainRunName = 20161117.beta0.A_WCYCL1850.ne30_oEC.edison
-# referenceRunName is the name of a reference run to compare against (or None
-# to turn off comparison with a reference, e.g. if no reference case is
-# available)
-referenceRunName = None
 # preprocessedReferenceRunName is the name of a reference run that has been
 # preprocessed to compare against (or None to turn off comparison).  Reference
 # runs of this type would have preprocessed results because they were not
@@ -20,11 +16,6 @@ preprocessedReferenceRunName = B1850C5_ne30_v0.4
 
 # directory containing model results
 baseDirectory = /scratch2/scratchdirs/golaz/ACME_simulations/20161117.beta0.A_WCYCL1850.ne30_oEC.edison/run
-# names of namelist and streams files.  If not in baseDirectory, give full path
-oceanNamelistFileName = mpas-o_in
-oceanStreamsFileName = streams.ocean
-seaIceNamelistFileName = mpas-cice_in
-seaIceStreamsFileName = streams.cice
 
 # names of ocean and sea ice meshes (e.g. EC60to30, QU240, RRS30to10, etc.)
 mpasMeshName = EC60to30
@@ -122,3 +113,15 @@ baseDirectory = /global/project/projectdirs/acme/observations/SeaIce
 # directory where ocean reference simulation results are stored
 baseDirectory =  /global/project/projectdirs/acme/ACMEv0_lowres/B1850C5_ne30_v0.4/ice/postprocessing
 
+[timeSeriesSeaIceAreaVol]
+## options related to plotting time series of sea ice area and volume
+
+# plot on polar plot
+polarPlot = False
+
+[streamfunctionMOC]
+## options related to plotting the streamfunction of the meridional overturning 
+## circulation (MOC)
+
+# Mask file for ocean basin regional computation
+regionMaskFiles = /global/project/projectdirs/acme/mapping/grids/EC60to30v1_SingleRegionAtlanticWTransportTransects_masks.nc