Changeset 3249 for trunk/LMDZ.MARS

Timestamp:

Mar 4, 2024, 8:23:10 AM (9 months ago)

Author:

llange

Message:

Mars PCM
Fixing bug in initestphys: the soil depth were initialized using the old grid
(but with 57 layers). It was not detected in soil_settings -to be investigated

• and therefore soil thermal inertia, when adding subsurface ice, were

completely wrong
LL

Location:

trunk/LMDZ.MARS

Files:

• changelog.txt (modified) (1 diff)
• libf/phymars/dyn1d/init_testphys1d_mod.F90 (modified) (2 diffs)

trunk/LMDZ.MARS/changelog.txt

@@ -4535 +4535 @@
 evolution without  adsorption
 
+
+== 04/03/2024 == LL
+Fixing bug in initestphys: the soil depth were initialized using the old grid
+(but with 57 layers). It was not detected in soil_settings -to be investigated
+- and therefore soil thermal inertia, when adding subsurface ice, were
+completely wrong

trunk/LMDZ.MARS/libf/phymars/dyn1d/init_testphys1d_mod.F90

@@ -661 +661 @@
     ! Initialize depths
     ! -----------------
-    do isoil = 1,nsoil
-        layer(isoil) = 2.e-4*(2.**(isoil - 1)) ! layer depth
+    do isoil = 0,nsoil - 1
+       mlayer(isoil) = 2.e-4*(1 + isoil**2.9*(1 - exp(-real(isoil)/20.))) ! mid layer depth
     enddo
+    do isoil = 1,nsoil - 1
+       layer(isoil) = (mlayer(isoil) + mlayer(isoil - 1))/2
+    enddo
+    layer(nsoil) = 2*mlayer(nsoil - 1) - mlayer(nsoil - 2)
 
     ! Creating the new soil inertia table if there is subsurface ice:
@@ -706 +710 @@
 if (.not. therestartfi) qsoil = 0.
 
-! Initialize depths
-! -----------------
-do isoil = 0,nsoil - 1
-    mlayer(isoil) = 2.e-4*(1 + isoil**2.9*(1 - exp(-real(isoil)/20.))) ! mid layer depth
-enddo
-do isoil = 1,nsoil - 1
-    layer(isoil) = (mlayer(isoil) + mlayer(isoil - 1))/2
-enddo
-layer(nsoil) = 2*mlayer(nsoil - 1) - mlayer(nsoil - 2)
-
 ! Initialize traceurs
 ! -------------------