hotfix to make python plotting work for the mpas2024052700 case (#191)

NOAA-EMC · Oct 11, 2024 · 261fab3 · 261fab3
1 parent 689818a
commit 261fab3
Show file tree

Hide file tree

Showing 3 changed files with 5 additions and 4 deletions.
diff --git a/rrfs-test/scripts/link_mpasjedi_expr.sh b/rrfs-test/scripts/link_mpasjedi_expr.sh
@@ -43,6 +43,7 @@ sed -e "s#@RDASApp@#${RDASApp}#" ${RDASApp}/rrfs-test/scripts/templates/mpasjedi
 sed -e "s#@RDASApp@#${RDASApp}#" ${RDASApp}/rrfs-test/scripts/templates/mpasjedi_expr/run_jedi.sh > run_jedi.sh
 cp ${RDASApp}/rrfs-test/ush/colormap.py .
 cp ${RDASApp}/rrfs-test/ush/mpasjedi_increment_singleob.py .
+cp ${RDASApp}/rrfs-test/ush/mpasjedi_increment_fulldom.py .
 cp ${RDASApp}/rrfs-test/ush/mpasjedi_spread.py .
 
 mkdir -p data

diff --git a/rrfs-test/ush/mpasjedi_increment_fulldom.py b/rrfs-test/ush/mpasjedi_increment_fulldom.py
@@ -19,7 +19,7 @@
 
 ############ USER INPUT ##########################################################
 plot_var = "Increment"
-lev = 1                 # 1=sfc, 55=toa
+lev = 23                 # 1=sfc, 55=toa
 clevmax_incr = 2.0      # max contour level for colorbar increment plots
 decimals = 2            # number of decimals to round for text boxes
 plot_box_width = 100.   # define size of plot domain (units: lat/lon degrees)
@@ -119,7 +119,7 @@ def plot_T_inc(var_n, clevmax):
 
 # Plot the data using triangulation
 clevs, cm, units, longname = plot_T_inc(jedi_inc, clevmax_incr)
-c1 = m1.tricontourf(triang, jedi_inc, clevs, cmap=cm, extend='both', transform=ccrs.PlateCarree())
+c1 = m1.tricontourf(triang, jedi_inc, clevs, cmap=cm, extend='both')
 
 # Add colorbar
 cbar1 = fig.colorbar(c1, orientation="horizontal", fraction=0.046, pad=0.07)

diff --git a/rrfs-test/ush/mpasjedi_increment_singleob.py b/rrfs-test/ush/mpasjedi_increment_singleob.py
@@ -19,7 +19,7 @@
 
 ############ USER INPUT ##########################################################
 plot_var = "Increment"
-lev = 1                 # 1=sfc, 55=toa
+lev = 23                # 1=sfc, 55=toa
 clevmax_incr = 2.0      # max contour level for colorbar increment plots
 decimals = 2            # number of decimals to round for text boxes
 plot_box_width = 8.     # define size of plot domain (units: lat/lon degrees)
@@ -169,7 +169,7 @@ def plot_T_inc(var_n, clevmax):
 
 # Plot the data using triangulation
 clevs, cm, units, longname = plot_T_inc(jedi_inc, clevmax_incr)
-c1 = m1.tricontourf(triang, jedi_inc, clevs, cmap=cm, extend='both', transform=ccrs.PlateCarree())
+c1 = m1.tricontourf(triang, jedi_inc, clevs, cmap=cm, extend='both')
 
 # Scatter the single ob location
 m1.scatter(singleob_lon, singleob_lat, color='g', marker='o', s=2)