Changeset 3115 for trunk/LMDZ.MARS/libf/phymars/dyn1d

Timestamp:

Nov 3, 2023, 6:02:56 PM (20 months ago)

Author:

llange

Message:

Mars PCM

• Introducing the possibily to compute water adsorption / desorption -routine soilwater.F90) -. FOR NOW IT WORKS ONLY IN 1D, DON'T TEST IN 3D (by default, adsorption is not called, if not specified in the callphys.def). By default, when using the adsorption, the number of subtimestep used in the water subimation scheme is 1 (otherwise, it crashes)

• New grid for the soil layers (better resolution) as it is needed to solve the diffusion equations. It does not increase the CPU time.

PYM & LL

File:

• trunk/LMDZ.MARS/libf/phymars/dyn1d/init_testphys1d_mod.F90 (modified) (1 diff)

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

@@ -670 +670 @@
 ! -----------------
 do isoil = 0,nsoil - 1
-    mlayer(isoil) = 2.e-4*(2.**(isoil - 0.5)) ! mid-layer depth
-    layer(isoil + 1) = 2.e-4*(2.**isoil)      ! layer depth
+    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