Changeset 3113 for trunk/LMDZ.MARS/libf/phymars/physiq_mod.F

Timestamp:

Nov 3, 2023, 2:36:01 PM (13 months ago)

Author:

llange

Message:

Mars PCM
Introducing qsoil to model H2O adsorption/desorption in the subsurface. For now, I've fixed the number of tracers in the subsurface to three (H2O vap, H2O ice, H2O ads).
The model of adsorption/desorption will follow later.
LL

File:

• trunk/LMDZ.MARS/libf/phymars/physiq_mod.F (modified) (3 diffs)

trunk/LMDZ.MARS/libf/phymars/physiq_mod.F

@@ -48 +48 @@
       use comsoil_h, only: inertiedat, inertiesoil,! dat: soil thermal inertia for present climate, inertiesoil is the TI read in the start
      &                     tsoil, nsoilmx,!number of subsurface layers
-     &                     mlayer,layer ! soil mid layer depths
+     &                     mlayer,layer, ! soil mid layer depths
+     &                     nqsoil,qsoil  ! adsorption
       use geometry_mod, only: longitude, latitude, cell_area,
      &                        cell_area_for_lonlat_outputs,longitude_deg
@@ -599 +600 @@
 ! GCM. Read netcdf initial physical parameters.
          CALL phyetat0 ("startfi.nc",0,0,
-     &         nsoilmx,ngrid,nlayer,nq,
+     &         nsoilmx,ngrid,nlayer,nq,nqsoil,
      &         day_ini,time_phys,
      &         tsurf,tsoil,albedo,emis,
-     &         q2,qsurf,tauscaling,totcloudfrac,wstar,
+     &         q2,qsurf,qsoil,tauscaling,totcloudfrac,wstar,
      &         watercap,perenial_co2ice,
      &         def_slope,def_slope_mean,subslope_dist)
@@ -2621 +2622 @@
      .          icount,' date=',ztime_fin
             
-          call physdem1("restartfi.nc",nsoilmx,ngrid,nlayer,nq,
+          call physdem1("restartfi.nc",nsoilmx,ngrid,nlayer,nq,nqsoil,
      .                ptimestep,ztime_fin,
      .                tsurf,tsoil,inertiesoil,albedo,
-     .                emis,q2,qsurf,tauscaling,totcloudfrac,wstar,
-     .                watercap,perenial_co2ice)
+     .                emis,q2,qsurf,qsoil,tauscaling,totcloudfrac,
+     .                wstar,watercap,perenial_co2ice)
           ENDIF ! of IF (write_restart)