Changeset 2942 for trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars

Timestamp:

Apr 14, 2023, 6:25:51 PM (21 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/dynphy_lonlat/phymars/newstart.F (modified) (2 diffs)

trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/newstart.F

@@ -31 +31 @@
      &                     emis, hmons, summit, base, watercap,
      &               ini_surfdat_h_slope_var,end_surfdat_h_slope_var
-      use comsoil_h, only: inertiedat, layer, mlayer, nsoilmx, tsoil,
-     & ini_comsoil_h_slope_var, end_comsoil_h_slope_var
+      use comsoil_h, only: inertiedat, inertiesoil,layer, mlayer,
+     & nsoilmx, tsoil,ini_comsoil_h_slope_var, end_comsoil_h_slope_var
       use control_mod, only: day_step, iphysiq, anneeref, planet_type
       use geometry_mod, only: longitude,latitude,cell_area
@@ -1387 +1387 @@
             inertiedat(1:ngridmx,isoil)=surfithfi(1:ngridmx)
           enddo
+      do islope = 1,nslope
+        inertiesoil(:,:,islope) =  inertiedat(:,:)
+      enddo
           write(*,*)'OK: Soil thermal inertia has been reset to referenc
      &e surface values'