Added SubpolarNorthAtlantic region

valeriupredoi · Mar 25, 2024 · a2ec969 · a2ec969
1 parent a880296
commit a2ec969
Show file tree

Hide file tree

Showing 3 changed files with 9 additions and 1 deletion.
diff --git a/bgcval2/bgcvaltools/makeMask.py b/bgcval2/bgcvaltools/makeMask.py
@@ -245,6 +245,12 @@ def makeMask(name, newSlice, xt, xz, xy, xx, xd, debug=False):
             xy, 60., 80.).mask
         return mx
 
+    if newSlice in ['SubpolarNorthAtlantic', 'SPNA',]:
+        #  Based on  SPNA region here: https://www.nature.com/articles/s43247-021-00120-y#citeas
+        mx = np.ma.masked_outside(xx, -35., -10.).mask + np.ma.masked_outside(
+            xy, 40., 65.).mask
+        return mx
+
     if newSlice == 'AtlanticSOcean':
         mx = np.ma.masked_outside(xx, -40., 20.).mask + np.ma.masked_outside(
             xy, -50., -75.).mask

diff --git a/bgcval2/bgcvaltools/pftnames.py b/bgcval2/bgcvaltools/pftnames.py
@@ -527,6 +527,8 @@ def makeLongNameDict():
     lnd['Remainder'] = "Oligotrophic Gyres"
     lnd['ArcticOcean'] = "Arctic Ocean"
     lnd['NorthernSubpolarAtlantic'] = "Northern Subpolar Atlantic"
+    lnd['SPNA'] = "Subpolar North Atlantic"
+    lnd['SubpolarNorthAtlantic'] =  lnd['SPNA']
     lnd['NorthernSubpolarPacific'] = "Northern Subpolar Pacific"
 
     lnd['SouthernOcean'] = "Southern Ocean"

diff --git a/key_files/temperature.yml b/key_files/temperature.yml
@@ -21,7 +21,7 @@ data_convert    : NoChange
 data_tdict      : ZeroToZero
 
 layers          : Surface
-regions         : Global ignoreInlandSeas SouthernOcean ArcticOcean Equator10 NorthPacificOcean SouthPacificOcean NorthAtlanticOcean  SouthAtlanticOcean 
+regions         : Global ignoreInlandSeas SouthernOcean ArcticOcean Equator10 NorthPacificOcean SouthPacificOcean NorthAtlanticOcean  SouthAtlanticOcean SPNA