Skip to content
Snippets Groups Projects

Changes in layer dynamics and enthalpy calculations

Merged Fuchs, Niels requested to merge debug_summertime_expulsion into master
3 files
+ 937
934
Compare changes
  • Side-by-side
  • Inline
Files
3
+ 1
1
@@ -287,7 +287,7 @@ CONTAINS
@@ -287,7 +287,7 @@ CONTAINS
!##########################################################################################
!##########################################################################################
!Brine flux due to Expulsion
!Brine flux due to Expulsion
!##########################################################################################
!##########################################################################################
CALL expulsion_flux (thick, V_ex, Nlayer, N_active, psi_g, fl_m, m)
CALL expulsion_flux (thick, V_ex, Nlayer, N_active, psi_g, fl_m, m, H_abs)
IF (i .NE. 1) THEN !is disabled on first timestep to avoid expulsion due to unbalanced initial conditions
IF (i .NE. 1) THEN !is disabled on first timestep to avoid expulsion due to unbalanced initial conditions
!This ensures that the salinity is not effected if the initial mass of the layers is to high
!This ensures that the salinity is not effected if the initial mass of the layers is to high
CALL mass_transfer (Nlayer, N_active, T, H_abs, S_abs, S_bu, T_bottom, S_bu_bottom, fl_m)
CALL mass_transfer (Nlayer, N_active, T, H_abs, S_abs, S_bu, T_bottom, S_bu_bottom, fl_m)
Loading