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

Timestamp:

Mar 5, 2024, 11:36:32 AM (12 months ago)

Author:

llange

Message:

Mars PCM
Following -r 3098; cleaning of vdfic for the management of subsurface water ice.
Fixing some errors (wrong interpolation to compute the water ice temperature, wrong boundary conditions to compute qvap(1))
LL

File:

• trunk/LMDZ.MARS/libf/phymars/dyn1d/init_testphys1d_mod.F90 (modified) (4 diffs)

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

@@ -42 +42 @@
 use nonoro_gwd_ran_mod,       only: du_nonoro_gwd, dv_nonoro_gwd
 use conf_phys_mod,            only: conf_phys
-use paleoclimate_mod,         only: paleoclimate, ini_paleoclimate_h, end_paleoclimate_h
+use paleoclimate_mod,         only: paleoclimate, ini_paleoclimate_h, end_paleoclimate_h, h2o_ice_depth
 use comslope_mod,             only: nslope, subslope_dist, ini_comslope_h, end_comslope_h
 use co2condens_mod,           only: CO2cond_ps
@@ -541 +541 @@
     call getin("subsurface_ice_depth",ice_depth)
     write(*,*) " subsurface_ice_depth = ",ice_depth
-
+    h2o_ice_depth(:,:) = ice_depth
     z0(1) = z0_default ! default value for roughness
     write(*,*) 'Surface roughness length z0 (m)?'
@@ -668 +668 @@
 
 if (.not. therestartfi) qsoil = 0.
-
 
 if (.not. therestartfi) then
@@ -716 +715 @@
         inertiedat(1,:) = inertiedat(1,1) ! soil thermal inertia
     endif ! ice_depth > 0
-
+    
     do isoil = 1,nsoil
         inertiesoil(1,isoil,:) = inertiedat(1,isoil)