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

Timestamp:

Feb 20, 2024, 10:17:54 PM (11 months ago)

Author:

llange

Message:

Mars PCM
Fixing a bug in 1D: atm_wat_profile should evolve freely when considering
relaxation, and should not be set to qsat in cold conditions.
LL

File:

• trunk/LMDZ.MARS/libf/phymars/dyn1d/testphys1d.F90 (modified) (1 diff)

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

@@ -204 +204 @@
             else ! Relaxation towards the value atm_wat_profile with relaxation time atm_wat_tau
                 q(1,:,igcm_h2o_vap) = atm_wat_profile*g/psurf + (q(1,:,igcm_h2o_vap) - atm_wat_profile*g/psurf)*dexp(-dttestphys/atm_wat_tau)
-                q(1,:,igcm_h2o_vap) = min(zqsat,q(1,:,igcm_h2o_vap))
                 q(1,:,igcm_h2o_ice) = 0. ! reset h2o ice
             endif