Changeset 2887 for trunk/LMDZ.MARS/libf/phymars/phyetat0_mod.F90

Timestamp:

Feb 1, 2023, 9:34:37 AM (22 months ago)

Author:

llange

Message:

Mars PCM
Include geothermal flux as boundary conditions when solving the
diffusion equation for the soil temperature
By default, the geothermal flux is set to 0. To have geothermal flux >
0., it must be specified in the startfi. (I suggest to use a value of
30e-3 in this case)
Will be use with the PEM for studies of "paleosubsurfaces"
LL

File:

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

trunk/LMDZ.MARS/libf/phymars/phyetat0_mod.F90

@@ -24 +24 @@
   use dust_param_mod, only: dustscaling_mode
   USE ioipsl_getin_p_mod, ONLY : getin_p
-
+  use comsoil_h, only: flux_geo
   implicit none
   
@@ -644 +644 @@
   ! as well as thermal inertia and volumetric heat capacity
   call soil_settings(nid_start,ngrid,nsoil,tsurf,tsoil,indextime)
+else
+    flux_geo(:) = 0.
 endif ! of if (startphy_file)