Changeset 3249 for trunk/LMDZ.MARS/libf

Timestamp:

Mar 4, 2024, 8:23:10 AM (10 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

File:

• trunk/LMDZ.MARS/libf/phymars/dyn1d/init_testphys1d_mod.F90 (modified) (2 diffs)

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
 ! -------------------