From 23be5ffd601f7363b027f7eff5905f1ec0447680 Mon Sep 17 00:00:00 2001
From: "Fuchs, Niels" <niels.fuchs@uni-hamburg.de>
Date: Fri, 3 May 2024 09:07:21 +0000
Subject: [PATCH] correct enthalpy balance in expulsion
---
mo_grotz.f90 | 1 +
mo_layer_dynamics.f90 | 2 +-
mo_mass.f90 | 2 +-
3 files changed, 3 insertions(+), 2 deletions(-)
diff --git a/mo_grotz.f90 b/mo_grotz.f90
index bf15182..b0f69ca 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 a424313..6602370 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 895c13c..7cc3f15 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
--
GitLab