Fix uninitialized variables in find_surface when NDIFF_INTERIOR_ONLY …

…= False
NOAA-GFDL · Sep 26, 2019 · af41095 · af41095
1 parent e806fba
commit af41095
Showing 1 changed file with 3 additions and 0 deletions.
diff --git a/src/tracer/MOM_neutral_diffusion.F90 b/src/tracer/MOM_neutral_diffusion.F90
@@ -1177,6 +1177,9 @@ subroutine find_neutral_surface_positions_discontinuous(CS, nk, Pres_l, hcol_l,
     p_init_L = zeta_bot_L; p_init_R = zeta_bot_R
     lastP_left = zeta_bot_L; lastP_right = zeta_bot_R
     kl_left = k_bot_L; kl_right = k_bot_R
+  else
+    k_init_L = 1  ; k_init_R = 1
+    p_init_L = 0. ; p_init_R = 0.
   endif
   ! Loop over each neutral surface, working from top to bottom
   neutral_surfaces: do k_surface = 1, ns