Changes in layer dynamics and enthalpy calculations

• /mo_grotz.f90
• /mo_layer_dynamics.f90: made middle layer thickness more dynamic to improve handling of ice core initialization runs
• /mo_mass.f90: remove brine mass from total enthalpy in expulsion function
• /mo_data.f90: corrected confused enthalpy descriptions

mo_data.f90

@@ -21,6 +21,7 @@
 !! @par Revision History
 !! Initialized by Philipp Griewank, IMPRS (2010-07-14) \n
 !! Add several variables by Niels Fuchs, MPIMET (2017-03-01)
+!! edited by Niels Fuchs, UHH (2024-05-03), corrected Enthalpy declaration (before specific and total enthalpy were confused)
 
 !!
 MODULE mo_data
@@ -31,8 +32,8 @@ MODULE mo_data
 
   
   !----Arrays
-  REAL(wp), DIMENSION(:), ALLOCATABLE   :: H                 !<  Enthalpy [J]
-  REAL(wp), DIMENSION(:), ALLOCATABLE   :: H_abs             !<  specific Enthalpy [J/kg]
+  REAL(wp), DIMENSION(:), ALLOCATABLE   :: H                 !<  specific Enthalpy [J/kg]
+  REAL(wp), DIMENSION(:), ALLOCATABLE   :: H_abs             !<  total Enthalpy [J]
   REAL(wp), DIMENSION(:), ALLOCATABLE   :: Q                 !<  Heat in layer [J]
   REAL(wp), DIMENSION(:), ALLOCATABLE   :: fl_Q              !<  Heat flux between layers [J/s]
   REAL(wp), DIMENSION(:), ALLOCATABLE   :: T                 !<  Temperature [C]