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

Timestamp:

Nov 3, 2023, 2:36:01 PM (20 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

Location:

trunk/LMDZ.MARS/libf/phymars/dyn1d

Files:

• init_testphys1d_mod.F90 (modified) (2 diffs)
• testphys1d.F90 (modified) (2 diffs)

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

@@ -20 +20 @@
 use iostart,                  only: open_startphy, get_var, close_startphy
 use physics_distribution_mod, only: init_physics_distribution
-use comsoil_h,                only: volcapa, nsoilmx, inertiesoil, inertiedat, layer, mlayer, flux_geo, tsoil
+use comsoil_h,                only: volcapa, nsoilmx, inertiesoil, inertiedat, layer, mlayer, flux_geo, tsoil, qsoil
 use comvert_mod,              only: ap, bp, aps, bps, pa, preff, presnivs, pseudoalt, scaleheight
 use dimradmars_mod,           only: tauvis, totcloudfrac, albedo
@@ -661 +661 @@
 flux_geo(:,:) = flux_geo_tmp
 
+! Initialize soil content
+! -----------------
+if (.not. startfiles_1D) then
+    qsoil(:,:,:,:) = 0.
+endif
+
 ! Initialize depths
 ! -----------------

trunk/LMDZ.MARS/libf/phymars/dyn1d/testphys1d.F90

@@ -1 +1 @@
 PROGRAM testphys1d
 
-use comsoil_h,           only: inertiedat, inertiesoil, nsoilmx, tsoil
+use comsoil_h,           only: inertiedat, inertiesoil, nsoilmx, tsoil, nqsoil, qsoil
 use surfdat_h,           only: albedodat, perenial_co2ice, watercap, tsurf, emis, qsurf
 use comslope_mod,        only: def_slope, subslope_dist
@@ -127 +127 @@
                   llm,nq,dttestphys,float(day0),0.,cell_area,    &
                   albedodat,inertiedat,def_slope,subslope_dist)
-    call physdem1("startfi.nc",nsoilmx,ngrid,llm,nq,dttestphys,time,       &
-                  tsurf,tsoil,inertiesoil,albedo,emis,q2,qsurf,tauscaling, &
+    call physdem1("startfi.nc",nsoilmx,ngrid,llm,nq,nqsoil,dttestphys,time,       &
+                  tsurf,tsoil,inertiesoil,albedo,emis,q2,qsurf,qsoil,tauscaling, &
                   totcloudfrac,wstar,watercap,perenial_co2ice)
 endif !(.not. therestartfi)