Only register N2, S2 diagnostics when the configuration supports them. #147

Author: Nicholas Hannah

src/parameterizations/lateral/MOM_lateral_mixing_coeffs.F90

@@ -824,14 +824,17 @@ subroutine VarMix_init(Time, G, param_file, diag, CS)
        'Length scale squared for mixing coefficient, at u-points', 'm^2')
     CS%id_L2v = register_diag_field('ocean_model', 'L2v', diag%axesCv1, Time, &
        'Length scale squared for mixing coefficient, at v-points', 'm^2')
-    CS%id_N2_u = register_diag_field('ocean_model', 'N2_u', diag%axesCui, Time, &
-       'Square of Brunt-Vaisala frequency, N^2, at u-points, as used in Visbeck et al.', 's^-2')
-    CS%id_N2_v = register_diag_field('ocean_model', 'N2_v', diag%axesCvi, Time, &
-       'Square of Brunt-Vaisala frequency, N^2, at v-points, as used in Visbeck et al.', 's^-2')
-    CS%id_S2_u = register_diag_field('ocean_model', 'S2_u', diag%axesCu1, Time, &
-       'Depth average square of slope magnitude, S^2, at u-points, as used in Visbeck et al.', 's^-2')
-    CS%id_S2_v = register_diag_field('ocean_model', 'S2_v', diag%axesCv1, Time, &
-       'Depth average square of slope magnitude, S^2, at v-points, as used in Visbeck et al.', 's^-2')
+
+    if (CS%use_stored_slopes) then
+      CS%id_N2_u = register_diag_field('ocean_model', 'N2_u', diag%axesCui, Time, &
+         'Square of Brunt-Vaisala frequency, N^2, at u-points, as used in Visbeck et al.', 's^-2')
+      CS%id_N2_v = register_diag_field('ocean_model', 'N2_v', diag%axesCvi, Time, &
+         'Square of Brunt-Vaisala frequency, N^2, at v-points, as used in Visbeck et al.', 's^-2')
+      CS%id_S2_u = register_diag_field('ocean_model', 'S2_u', diag%axesCu1, Time, &
+         'Depth average square of slope magnitude, S^2, at u-points, as used in Visbeck et al.', 's^-2')
+      CS%id_S2_v = register_diag_field('ocean_model', 'S2_v', diag%axesCv1, Time, &
+         'Depth average square of slope magnitude, S^2, at v-points, as used in Visbeck et al.', 's^-2')
+    endif
   endif
 
   call wave_speed_init(Time, G, param_file, diag, CS%wave_speed_CSp)