Plot pressure adjustment in grounded ice region

Author: cbegeman

compass/ocean/tests/isomip_plus/viz/__init__.py

@@ -99,6 +99,13 @@ def run(self):
                                   True, vmin=min_column_thickness + 1e-10,
                                   vmax=700, cmap_set_under='r',
                                   cmap_scale='log')
+        plotter.plot_horiz_series(dsOut.surfacePressure,
+                                  'surfacePressure', 'surfacePressure',
+                                  True, vmin=1e5, vmax=1e7, cmap_scale='log')
+        delSurfacePressure = dsOut.landIcePressure - dsOut.surfacePressure
+        plotter.plot_horiz_series(delSurfacePressure,
+                                  'delSurfacePressure', 'delSurfacePressure',
+                                  True, vmin=1e5, vmax=1e7, cmap_scale='log')
 
         if 'tidal' in expt:
             delssh = dsOut.ssh - dsOut.ssh[0, :]
@@ -111,10 +118,10 @@ def run(self):
         wct_thin = wct[:, idx_thin]
         wct_mean = wct_thin.mean(dim='nCells').values
         time = dsOut.daysSinceStartOfSim.values
-        fig = plt.figure()
+        plt.figure()
         plt.plot(time, wct_mean, '.')
-        fig.set_xlabel('Time (days)')
-        fig.set_ylabel('Mean thickness of thin film (m)')
+        plt.xlabel('Time (days)')
+        plt.ylabel('Mean thickness of thin film (m)')
         plt.savefig('wct_thin_t.png')
         plt.close()