Fix bad regridding in dynamics extension and fix bad description

dynamics/dynamics/cisma_racmo23_scatter_standalone.py

@@ -22,7 +22,7 @@ def make_plot(config=None, out_path='.',
     # which were saved in my local directory NCLplot
 
     # ---------------- Read Data -----------------------------------
-    input_file1 = os.path.join(racmo_path, 'climos/racmo23_GRN_monthly.smb.1980-1999.ANN.nc')
+    input_file1 = os.path.join(racmo_path, 'racmo23.ANN_with_latlon_smb.nc')
     input_file2 = os.path.join(cism_path, 'CISM-Albany/postproc/remapped_cism/CISM-Albany.acab.remap2racmo.ANN.nc')
     input_file3 = os.path.join(racmo_path, 'RACMO23_masks_ZGRN11.nc')
     # --------------------------------------------------------------

dynamics/dynamics/cisma_racmo23_standalone.py

@@ -19,7 +19,7 @@ def make_plot(config=None, out_path='.',
     # which were saved in my local directory NCLplot
 
     # ---------------- Read Data -----------------------------------
-    input_file1 = os.path.join(racmo_path, 'climos/racmo23_GRN_monthly.smb.1980-1999.ANN.nc')
+    input_file1 = os.path.join(racmo_path, 'racmo23.ANN_with_latlon_smb.nc')
     input_file2 = os.path.join(cism_path, 'CISM-Albany/postproc/remapped_cism/CISM-Albany.acab.remap2racmo.ANN.nc')
     input_file3 = os.path.join(racmo_path, 'RACMO23_masks_ZGRN11.nc')
     input_file4 = os.path.join(cism_path, 'CISM-Albany/postproc/albany.out.tstep.ANN_acab.nc')
@@ -68,7 +68,7 @@ def make_plot(config=None, out_path='.',
     thk_mask = ma.masked_less(thk, 0.01)
     cisma_ori = ma.masked_array(cisma_ori, mask=thk_mask.mask)
 
-    # diff = cisma - racmo
+    diff = cisma - racmo
     # print(np.max(diff))
     # print(np.min(diff))
 
@@ -166,22 +166,21 @@ def make_plot(config=None, out_path='.',
     res22.lbLabelFontHeightF = 0.02  # Make fonts smaller.
     res22.pmLabelBarOrthogonalPosF = -0.04  # move label bar closer
 
-    res22.sfXArray = lon1
-    res22.sfYArray = lat1
+    res22.sfXArray = lon2
+    res22.sfYArray = lat2
 
     racmo_plot = Ngl.contour(wks, racmo, res1)
     cisma_plot = Ngl.contour(wks, cisma_ori, res2)
-    # diff_plot = Ngl.contour(wks,diff,res22)
-    diff_plot = Ngl.contour(wks, cisma, res22)
+    diff_plot = Ngl.contour(wks, diff, res22)
 
     # Creat multiple figures and draw, which now contains the elevation and temperature
     # "[1,3]" indicates 1 row, 3 columns.
-    # map_title = ["CISM-A","RACMO23", "CISM-A-RACMO23"]
+    map_title = ["CISM-A","RACMO23", "CISM-A--RACMO23"]
 
     nmap = 3
     plot = []
     for i in range(nmap):
-        # mres.tiMainString      = map_title[i]
+        mres.tiMainString = map_title[i]
         plot.append(Ngl.map(wks, mres))
 
     # Overlay everything on the map plot.

dynamics/dynamics_cisma.py

@@ -26,7 +26,7 @@
 # OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
 # OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
 
-"""An analysis of CESM's cloud cover over Greenland."""
+"""An analysis of CISM-Albany's Greenland dynamics."""
 
 from __future__ import absolute_import, division, print_function, unicode_literals