Skip to content
Snippets Groups Projects
Commit 0fff4552 authored by Evgeny Mortikov's avatar Evgeny Mortikov
Browse files

removing redundant closure option

parent a8809157
No related branches found
No related tags found
No related merge requests found
......@@ -295,7 +295,6 @@ program obl_main
do while (time_current < time_end )
! ----------------------------------------------------------------------------
if (closure_mode.ne.5) then
!< define fluxes & dynamic scales [surface]
! ----------------------------------------------------------------------------
......@@ -396,14 +395,6 @@ program obl_main
call solve_scalar_eq (Salin_dev, Kh, nz, dz, dt, flux_sal_surf_res, flux_sal_bot_res, f_sal_right)
call solve_vector_eq (U, Km, nz, dz, dt, flux_u_surf_res, flux_u_bot_res, f_u_right)
call solve_vector_eq (V, Km, nz, dz, dt, flux_v_surf_res, flux_v_bot_res, f_v_right)
!else if (closure_mode.eq.5) then
! call solve_state_eq (Rho, Theta_dev, Salin_dev, nz)
! call vermix_init(rho, u, v, Tau_x_surf(i), Tau_y_surf(i), Ri_grad, Kh, Km)
! call solve_scalar_eq (Theta_dev, Kh, nz, dz, dt, flux_heat_surf(i), flux_heat_bot(i), f_heat_right)
! 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)
endif
!> advance time
i = i + 1
......@@ -435,10 +426,12 @@ program obl_main
endif
enddo
if (output_mode > 0) then
output_Theta%data(:,1:output_Theta%num) = output_Theta%data(:,1:output_Theta%num) + Theta_ref
output_Salin%data(:,1:output_Salin%num) = output_Salin%data(:,1:output_Salin%num) + Salin_ref
output_TinC%data(:,1:output_TinC%num) = output_TinC%data(:,1:output_TinC%num) + Theta_ref - 273.15
endif
if (output_mode.eq.1) then
write(*, *) ' >> writing netcdf output ...'
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment