Changeset 2942 for trunk/LMDZ.MARS/libf/phymars/comsoil_h.F90

Timestamp:

Apr 14, 2023, 6:25:51 PM (20 months ago)

Author:

llange

Message:

Mars PCM
Add the possibility to use a different thermal inertia (field
'inertiesoil') than inertiedat in the PCM (for paleoclimate studies). By defaut, if there is not
inertiesoil, inertiedat is used. Soil_tifeedback still work with
inertiedat
Newstart adapted, start2archive will be modified by Romain
Vandemeulebrouck.
LL

File:

• trunk/LMDZ.MARS/libf/phymars/comsoil_h.F90 (modified) (5 diffs)

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

@@ -8 +8 @@
   real,save,allocatable,dimension(:) :: layer      ! soil layer depths
   real,save,allocatable,dimension(:) :: mlayer     ! soil mid-layer depths
-  real,save,allocatable,dimension(:,:) :: inertiedat ! soil thermal inertia
+  real,save,allocatable,dimension(:,:) :: inertiedat ! soil thermal inertia for present climate
+  real,save,allocatable,dimension(:,:,:) :: inertiesoil ! soil thermal inertia
   real,save :: volcapa    ! soil volumetric heat capacity
        ! NB: volcapa is read fromn control(35) from physicq start file
@@ -38 +39 @@
     allocate(layer(nsoilmx)) !soil layer depths
     allocate(mlayer(0:nsoilmx-1)) ! soil mid-layer depths
-    allocate(inertiedat(ngrid,nsoilmx)) ! soil thermal inertia
- 
+    allocate(inertiedat(ngrid,nsoilmx)) ! soil thermal inertia for present climate
+    allocate(inertiesoil(ngrid,nsoilmx,nslope)) ! soil thermal inertia
     allocate(tsoil(ngrid,nsoilmx,nslope)) ! soil temperatures
-
     allocate(mthermdiff(ngrid,0:nsoilmx-1,nslope))
     allocate(thermdiff(ngrid,nsoilmx-1,nslope))
@@ -60 +60 @@
     if (allocated(mlayer)) deallocate(mlayer)
     if (allocated(inertiedat)) deallocate(inertiedat)
+    if (allocated(inertiesoil)) deallocate(inertiesoil)
     if (allocated(tsoil)) deallocate(tsoil)
     if (allocated(mthermdiff)) deallocate(mthermdiff)
@@ -77 +78 @@
  
     allocate(tsoil(ngrid,nsoilmx,nslope)) ! soil temperatures
+    allocate(inertiesoil(ngrid,nsoilmx,nslope)) ! soil thermal inertia
     allocate(mthermdiff(ngrid,0:nsoilmx-1,nslope))
     allocate(thermdiff(ngrid,nsoilmx-1,nslope))
@@ -91 +93 @@
 
     if (allocated(tsoil)) deallocate(tsoil)
+    if (allocated(inertiesoil)) deallocate(inertiesoil)
     if (allocated(mthermdiff)) deallocate(mthermdiff)
     if (allocated(thermdiff)) deallocate(thermdiff)