diff --git a/new_masks/basins.yml b/new_masks/basins.yml
index 201d6fe02166b44f491f270f472db1486961e59d..0dfd98d8b086b39e6975f4e29ca8e14ada98baff 100644
--- a/new_masks/basins.yml
+++ b/new_masks/basins.yml
@@ -1,3 +1,5 @@
+version: 1
+
 basins:
   Global_Ocean:
     - {lat: [-90, 90], lon: [-180, 180]}
@@ -429,10 +431,47 @@ basins:
   Nino4:
     - {lat: [-5, 5], lon: [-180, -150]}
     - {lat: [-5, 5], lon: [160, 180]}
-    
+
   AMV_North_Atlantic:
+    - {basin: Atlantic_Ocean}
     - {lat: [0, 60], lon:[-80, 0]}
-  
+
+  AMV_North_Atlantic2:
+    - {basin: Atlantic_Ocean}
+    - {lat: [30, 60], lon:[-80, 0]}
+
   AMV_trend:
     - {lat: [-60, 60], lon: [-180, 180]}
-    
\ No newline at end of file
+
+  Subpolar_Gyre_Box:
+    - {lat: [50, 65], lon:[-60, 0]}
+
+  West_Subpolar_Gyre_Box:
+    - {lat: [50, 65], lon:[-60, -30]}
+
+  East_Subpolar_Gyre_Box:
+    - {lat: [50, 65], lon:[-30, 0]}
+
+  Pacific_TPI1:
+    - {lat: [25, 45], lon:[140, 180]}
+    - {lat: [25, 45], lon:[-180, -145]}
+
+  Pacific_TPI2:
+    - {lat: [-10, 10], lon:[170, 180]}
+    - {lat: [-10, 10], lon:[-180, -90]}
+
+  Pacific_TPI3:
+    - {lat: [-50, -15], lon:[150, 180]}
+    - {lat: [-50, -15], lon:[-180, -160]}
+
+  Indian_dipole_west:
+    - {lat: [-10, 10], lon:[50,70]}
+
+  Indian_dipole_east:
+    - {lat: [-10, 0], lon:[90, 110]}
+
+  Labrador_Sea_conv:
+    - {lat: [55, 65], lon: [-56, -46]}
+
+  Labrador_Sea_conv2:
+    - {lat: [52, 65], lon: [-58, -43]}
diff --git a/new_masks/create_basins_mask.py b/new_masks/create_basins_mask.py
index 0f48da3dbe7f1d42fc930b7317810543c1e8c168..f2b4c6ca201840746048e950115b812a913a9249 100644
--- a/new_masks/create_basins_mask.py
+++ b/new_masks/create_basins_mask.py
@@ -72,6 +72,8 @@ def create_mask(path, version, config_file=None, plot=False):
         DimCoord(np.arange(lat.shape[1]), var_name='x'), (2, ))
     basins.add_dim_coord(
         DimCoord(np.arange(lat.shape[0]), var_name='y'), (1, ))
+    basins.attributes['version'] = int(config['version'])
+    basins.attributes['grid'] = version
     print('Basins cube:')
     print(basins)
     iris.save(
@@ -178,6 +180,6 @@ def intersection(mask1, mask2):
 
 
 if __name__ == "__main__":
-    mesh_path = "/esarchive/autosubmit/con_files"
-    # create_mask(mesh_path, 'Ec3.2_O1L75')
-    create_mask(mesh_path, 'Ec3.2_O25L75')
+    MESH_PATH = "/esarchive/autosubmit/con_files"
+    create_mask(MESH_PATH, 'Ec3.2_O1L75')
+    create_mask(MESH_PATH, 'Ec3.2_O25L75')
diff --git a/new_masks/load_modules.txt b/new_masks/load_modules.txt
new file mode 100644
index 0000000000000000000000000000000000000000..b4161944be7847cd135e8bfa1476457db82a93e5
--- /dev/null
+++ b/new_masks/load_modules.txt
@@ -0,0 +1,4 @@
+module load earthdiagnostics/3.1.1-foss-2015a-Python-3.7.3
+module load CDO/1.5.3-foss-2015a
+module load Iris/2.2.1-foss-2015a-Python-3.7.3
+module load numba/0.44.1-foss-2015a-Python-3.7.3