Merge pull request #434 from wwieder/clm-diags

LDF corrections for CUPiD

Author: wwieder

lib/ldf_variable_defaults.yaml

@@ -369,7 +369,7 @@ GPP: # Gross Primary Production
   mpl:
     colorbar:  
       label : "gC m$^{-2}$ d$^{-1}$"
-  pct_diff_contour_levels: [-100,-75,-50,-40,-30,-20,-10,-8,-6,-4,-2,0,2,4,6,8,10,20,30,40,50,75,100]
+  pct_diff_contour_levels: [-200,-100,-75,-50,-25,-10,-5,0,5,10,25,50,75,100,200]
   pct_diff_colormap: "PiYG"
   scale_factor_table: 0.000000365  #days to years, g/m2 to Pg globally
   avg_method: 'sum'

lib/regions_lnd.yaml

@@ -91,6 +91,14 @@ Europe:
   lat_bounds: [45, 60]
   lon_bounds: [-10, 30]
   region_category: Temperate 
+East Asia:
+  lat_bounds: [25, 45]
+  lon_bounds: [85, 120]
+  region_category: Temperate
+Eurasia:
+  lat_bounds: [15, 75]
+  lon_bounds: [30, 120]
+  region_category: Temperate
 Mediterranean:
   lat_bounds: [34, 45]
   lon_bounds: [-10, 30]

scripts/plotting/polar_map.py

@@ -339,6 +339,9 @@ def polar_map(adfobj):
                                     elif (var == 'FAREA_BURNED') and (s == 'SON'):
                                         print("\t\t Skipping FAREA_BURNED in SON plot")
                                         continue
+                                    elif (var == 'GPP') and (s == 'SON'):
+                                        print("\t\t Skipping GPP in SON plot")
+                                        continue
                                     # not working for regular grids?
                                     elif (var == 'NPP'):
                                         print("\t\t Skipping NPP polar plot")

scripts/plotting/regional_climatology.py

@@ -336,6 +336,10 @@ def regional_climatology(adfobj):
                     map_ax.set_extent([-140, -55, 10, 70],crs=ccrs.PlateCarree())
                 elif region_list[iReg]=='S America':
                     map_ax.set_extent([-100, -20, -45, 20],crs=ccrs.PlateCarree())
+                elif region_list[iReg]=='East Asia':
+                    map_ax.set_extent([70, 145, 5, 55],crs=ccrs.PlateCarree())
+                elif region_list[iReg]=='Eurasia':
+                    map_ax.set_extent([20, 130, 10, 80],crs=ccrs.PlateCarree())
                 else: 
                     if ((box_south >= 30) & (box_east<=-5) ):
                         map_ax.set_extent([-180, 0, 30, 90],crs=ccrs.PlateCarree())

scripts/plotting/regional_timeseries.py

@@ -337,6 +337,10 @@ def regional_timeseries(adfobj):
                     map_ax.set_extent([-140, -55, 10, 70],crs=ccrs.PlateCarree())
                 elif region_list[iReg]=='S America':
                     map_ax.set_extent([-100, -20, -45, 20],crs=ccrs.PlateCarree())
+                elif region_list[iReg]=='East Asia':
+                    map_ax.set_extent([70, 145, 5, 55],crs=ccrs.PlateCarree())
+                elif region_list[iReg]=='Eurasia':
+                    map_ax.set_extent([20, 130, 10, 80],crs=ccrs.PlateCarree())
                 else: 
                     if ((box_south >= 30) & (box_east<=-5) ):
                         map_ax.set_extent([-180, 0, 30, 90],crs=ccrs.PlateCarree())