Changeset 3733 for trunk/LMDZ.MARS/libf/phymars/dyn1d/testphys1d.F90

Timestamp:

Apr 22, 2025, 6:08:46 PM (8 weeks ago)

Author:

jbclement

Message:

Mars PCM:
Fixing a sanity check for the prescribtion of the atmospheric water profile in 1D and making comments more explicit.
JBC

File:

• trunk/LMDZ.MARS/libf/phymars/dyn1d/testphys1d.F90 (modified) (3 diffs)

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

@@ -150 +150 @@
                      time,psurf,u,v,temp,ndt,ptif,pks,dttestphys,zqsat,dq,dqdyn,day0,day,gru,grv,w, &
                      play,plev,latitude,longitude,cell_area,atm_wat_profile,atm_wat_tau)
-
+print*, 'yooooooooooooooooooo', q(1,1,igcm_h2o_vap)
 ! Write a "startfi" file
 ! ----------------------
@@ -163 +163 @@
                   totcloudfrac,wstar,watercap,perennial_co2ice)
 endif !(.not. therestartfi)
+print*, 'yooooooooooooooooooo', q(1,1,igcm_h2o_vap)
 
 !=======================================================================
@@ -200 +201 @@
     if (water) then
         call watersat(nlayer,temp,play,zqsat)
-        if (0. <= atm_wat_profile .and. atm_wat_profile <= 1.) then
-        ! If atmospheric water is monitored
-            if (atm_wat_tau < 0.) then ! Prescribed atm_wat_profile: wet if >0, dry if =0
+        if (0. <= atm_wat_profile) then ! If atmospheric water is monitored
+            if (atm_wat_tau < 0.) then ! Case for prescribed atm_wat_profile: wet if >0, dry if =0
                 q(1,:,igcm_h2o_vap) = min(zqsat,atm_wat_profile*g/psurf)
                 q(1,:,igcm_h2o_ice) = 0. ! reset h2o ice
-            else ! Relaxation towards the value atm_wat_profile with relaxation time atm_wat_tau
+            else !  Case for 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_ice) = 0. ! reset h2o ice