minor update

obl_main.f90

@@ -352,9 +352,6 @@ program obl_main
         !    call solve_scalar_eq (Salin_dev, Kh, nz, dz, dt, flux_sal_surf(i), flux_sal_bot(i), f_sal_right)
         !    call solve_vector_eq(U, Km, nz, dz, dt, flux_u_surf(i), flux_u_bot(i), f_u_right)
         !    call solve_vector_eq(V, Km, nz, dz, dt, flux_v_surf(i), flux_v_bot(i), f_v_right)
-
-            call get_mld(mld, N2, dz, nz, z)
-            call get_lab_mld(lab_mld, u_dynH, N2_0, time_current, z)
         endif
 
         !> advance time
@@ -363,6 +360,9 @@ program obl_main
 
         !> advance screen output
         if (mod(i, nscreen) == 0) then
+            call get_mld(mld, N2, dz, nz, z)
+            call get_lab_mld(lab_mld, u_dynH, N2_0, time_current, z)
+
             write(*, '(a,g0)') ' mld = ', mld
             write(*, '(a,g0)') ' Theta(surface) = ', Theta_dev(grid%cz) + Theta_ref
             write(*, '(a,g0,a,g0,a)') ' current time = ', time_current / 3600.0, ' HRS [ ', & 
@@ -374,6 +374,9 @@ program obl_main
         if (mod(i, noutput) == 0) then
             call push_state_vec(nz)
 
+            call get_mld(mld, N2, dz, nz, z)
+            call get_lab_mld(lab_mld, u_dynH, N2_0, time_current, z)
+            
             call push_value_to_tseries(output_mld, mld)
             call push_value_to_tseries(output_mld_ref, lab_mld)