Changeset 3098 for trunk/LMDZ.MARS/libf/phymars/dyn1d/init_testphys1d_mod.F90

Timestamp:

Oct 23, 2023, 7:02:06 PM (14 months ago)

Author:

evos

Message:

Mars PCM

We added the to module vdifc the possibily of subsurface intercation, mostly to have the option of buried glaicer that can lose ice and create polar layers.
if no ice is on the surface the ice the atmosphere can directly interact with the subsurface ice.
if ice is on the surface it can interact with the subsurface ice.
For now it is not working with more then one subslope. I will add this feature soon
In order to use it there are two flags that needs to be set to true.
paleoclimate=.true.
lag_layer=.true.
in addition one need to give the following parametrs in startfi or run.def:
d_coef=X for the diffusion coeficent (4e-4 is the defualt)
<h2o_ice_depth=X> if its zero or lower there is no subsurface ice (default -1).
outputs: zdq_ssi, momentarily zdqsdif after the interaction with the SSIzdq_subtimestep,  momentarily flux after the interaction with the SSIzdq_ssi_frost momentarily zdqsdif after the interaction of the frost with the SSI
q1, atmospheric first layer water vapor quatity
qeq, SSI water vapor quantity

Eran Vos

File:

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

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

@@ -477 +477 @@
     call getin("subsurface_ice_depth",ice_depth)
 
+
     z0(1) = z0_default ! default value for roughness
     write(*,*) 'Surface roughness length z0 (m)?'
@@ -716 +717 @@
         write(*,*) 'Prescribed atmospheric water vapor profile'
         write(*,*) 'Unless it reaches saturation (maximal value)'
+    else if (atm_wat_profile .eq. 2) then
+        write(*,*) 'Prescribed atmospheric water vapor profile'
+        write(*,*) 'like present day northern midlatitude'
+        write(*,*) 'Unless it reaches saturation (maximal value)'
     else
         error stop 'Water vapor profile value not correct!'