diff --git a/mo_grotz.f90 b/mo_grotz.f90 index bf151820a78bf2cda4501d1fc0b4aa030acd915f..b0f69ca22ff5b2b3bd26425d26b22d02012edfdc 100644 --- a/mo_grotz.f90 +++ b/mo_grotz.f90 @@ -79,6 +79,7 @@ CONTAINS !! Basic thermodynamics and layer_dynamics for fixed boundaries seem stable, backup made. by griewank (2010-08-10) \n !! Add some more outputs, changed routine names and arguments with respect to newly introduces flags by Niels Fuchs, MPIMET (2017-03-01) \n !! Added a bit of description with the run down of what happends by Philipp Griewank, Uni K (2018-08-08) + !! Edited by Niels Fuchs, UHH (2024-05-03), correct layer enthalpy balance in brine expulsion routine SUBROUTINE grotz () USE mo_parameters diff --git a/mo_layer_dynamics.f90 b/mo_layer_dynamics.f90 index a4243133708de8ca39dd7746e3c022095246e07f..660237010555a73deeca3a5505916916c6e75282 100644 --- a/mo_layer_dynamics.f90 +++ b/mo_layer_dynamics.f90 @@ -185,7 +185,7 @@ CONTAINS !! !! @par Revision History !! Pasted by Philipp Griewank, IMPRS (2011-05-10) - !! Edited by Niels Fuchs, UHH (2024-05-03), changed thick_0 to thick(i) to work better with ice core initialization + !! Edited by Niels Fuchs, UHH (2024-05-03), changed thick_0 to thick(k) in middle layers to work better with ice core initialization SUBROUTINE top_melt (Nlayer,N_active,N_middle,N_top,thick_0,m,S_abs,H_abs,thick,N_bgc,bgc_abs) diff --git a/mo_mass.f90 b/mo_mass.f90 index 895c13c1c767de8b52170efe41ab92f9fe2d16ac..7cc3f158e66ef461014c247fa751912d9b23d8b9 100644 --- a/mo_mass.f90 +++ b/mo_mass.f90 @@ -122,7 +122,7 @@ CONTAINS INTEGER::k - H = H_abs(:)/m(:) + H = H_abs(:) / m(:) fl_m(1:Nlayer+1) = 0._wp fl_m(2) = -V_ex(1)*rho_l DO k = 2,N_active