Changeset 3527 for trunk/LMDZ.COMMON/libf/evolution/NS_fast_subs_univ.F90

Timestamp:

Nov 20, 2024, 3:53:19 PM (3 days ago)

Author:

jbclement

Message:

PEM:

• Removing redundant Norbert Schorghofer's subroutines/parameters (follow-up of r3526);
• Making all Norbert Schorghofer's subroutines with modern explicit interface via modules;
• Cleaning of "glaciers_mod.F90";
• Optimization for the computation of ice density according to temperature by using a function.

JBC

File:

• trunk/LMDZ.COMMON/libf/evolution/NS_fast_subs_univ.F90 (modified) (8 diffs)

trunk/LMDZ.COMMON/libf/evolution/NS_fast_subs_univ.F90

@@ +1 @@
+MODULE fast_subs_univ
+
+implicit none
+
+CONTAINS
+
 !*****************************************************
 ! Commonly used subroutines for fast method
@@ -19 +25 @@
 !  this is not the true (final) equilibrium depth
 !***********************************************************************
-  use allinterfaces, except_this_one => equildepth
   implicit none
   integer, intent(IN) :: nz
@@ -54 +59 @@
 !  B = Diff*bigstep/(porosity*icedensity)  [SI units]
 !***********************************************************************
-  use allinterfaces, except_this_one => depths_avmeth
   use math_mod, only: deriv2_simple, deriv1_onesided, deriv1, colint
+  use ice_table_mod, only: constriction
   implicit none
   integer, intent(IN) :: nz, typeT
@@ -178 +183 @@
 
 
-pure function constriction(porefill)
-! specify constriction function here, 0<=eta<=1
-  implicit none
-  real(8), intent(IN) :: porefill
-  real(8) eta, constriction
-  if (porefill<=0.) eta = 1.
-  if (porefill>0. .and. porefill<1.) then
-     ! eta = 1.
-     ! eta = 1-porefill
-     eta = (1-porefill)**2  ! Hudson et al., JGR, 2009
-  endif
-  if (porefill>=1.) eta = 0.
-  constriction = eta
-end function constriction
-
-
-
 pure subroutine icechanges(nz,z,typeF,avdrho,avdrhoP,ypp, &
      & Diff,porosity,icefrac,bigstep,zdepthT,porefill,typeG)
@@ -202 +190 @@
 ! a crucial subroutine
 !***********************************************************
-  use miscparameters, only : icedensity
-  use allinterfaces, except_this_one => icechanges
   use math_mod, only: colint
+  use ice_table_mod, only: rho_ice
   implicit none
   integer, intent(IN) :: nz, typeF, typeG
@@ -213 +200 @@
   real(8) B, beta, integ
 
-  B = Diff*bigstep*86400.*365.24/(porosity*icedensity)
+  B = Diff*bigstep*86400.*365.24/(porosity*927.)
+  !B = Diff*bigstep*86400.*365.24/(porosity*rho_ice(T,'h2o'))
 
   ! advance ice table, avdrho>0 is retreat
@@ -232 +220 @@
      if (typeP==typeT) then   ! new 2011-09-01
         beta = (1-icefrac)/(1-porosity)/icefrac
-        beta = Diff*bigstep*beta*86400*365.24/icedensity
+        beta = Diff*bigstep*beta*86400*365.24/927.
+        !beta = Diff*bigstep*beta*86400*365.24/rho_ice(T,'h2o')
         zdepthT = sqrt(2*beta*avdrho*18./8314. + zdepthT**2)
      endif
@@ -283 +272 @@
   end if
 end subroutine icechanges
+
+END MODULE fast_subs_univ