diff --git a/esmvaltool/recipes/recipe_ocean_quadmap_NEMO.yml b/esmvaltool/recipes/recipe_ocean_quadmap_NEMO.yml
index c65d7e40cb..4d0f03c1f1 100644
--- a/esmvaltool/recipes/recipe_ocean_quadmap_NEMO.yml
+++ b/esmvaltool/recipes/recipe_ocean_quadmap_NEMO.yml
@@ -64,3 +64,24 @@ diagnostics:
         control_model: {project: NEMO, dataset: GC5}
         exper_model: {project: NEMO, dataset: GC3p1}
         observational_dataset: {project: OBS2, dataset: observation}
+
+  diag_map_2:
+    description: Global surface quad plots sos
+    variables:
+      sos:
+        preprocessor: prep_quad_map
+        mip: Omon
+        additional_datasets:
+          - {project: NEMO, dataset: GC5, suite: cp656, model: nemo, freq: 1y, grid: T, timerange: '199501/200501'}
+          - {project: NEMO, dataset: GC3p1, suite: co957, model: nemo, freq: 1y, grid: T, timerange: '199501/200501'}
+      sos_obs:
+        mip: Omon
+        short_name: sos
+        additional_datasets:
+          - {project: OBS3, dataset: observation, suite: EN4_v1.1, model: nemo, freq: annualmean, grid: T}
+    scripts:
+      Global_Ocean_map:  # &Global_Ocean_map SOS
+        script: ocean/diagnostic_ocean_marineassess_sss.py
+        control_model: {project: NEMO, dataset: GC5}
+        exper_model: {project: NEMO, dataset: GC3p1}
+        observational_dataset: {project: OBS3, dataset: observation}