Changeset 3733 for trunk/LMDZ.MARS/libf/phymars/dyn1d/init_testphys1d_mod.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/init_testphys1d_mod.F90 (modified) (4 diffs)

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

@@ -779 +779 @@
 ! ---------------------------------------------------
 ! Adding an option to force atmospheric water values JN
-atm_wat_profile = -1. ! Default: free atm wat profile
+atm_wat_profile = -1. ! Default: free atmospheric water profile
 if (water) then
-    write(*,*)'Force atmospheric water vapor profile?'
+    write(*,*)'Force atmospheric water vapor profile (uniform column value in kg.m-2)?'
     call getin('atm_wat_profile',atm_wat_profile)
     write(*,*) 'atm_wat_profile = ', atm_wat_profile
@@ -789 +789 @@
     else if (abs(atm_wat_profile) < 1.e-15) then ! if == 0.
         write(*,*) 'Dry atmospheric water vapor profile'
-    else if (0. < atm_wat_profile .and. atm_wat_profile <= 1.) then
+    else if (0. < atm_wat_profile) then
         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
@@ -806 +802 @@
 atm_wat_tau = -1. ! Default: no time relaxation
 if (water) then
-    write(*,*) 'Relax atmospheric water vapor profile?'
+    write(*,*) 'Relax atmospheric water vapor profile with time constant (s)?'
     call getin('atm_wat_tau',atm_wat_tau)
     write(*,*) 'atm_wat_tau = ', atm_wat_tau
     if (atm_wat_tau < 0.) then
-        write(*,*) 'Atmospheric water vapor profile is not relaxed.'
+        write(*,*) 'Atmospheric water vapor profile is not relaxed (forcing).'
     else
-        if (0. <= atm_wat_profile .and. atm_wat_profile <= 1.) then
+        if (0. <= atm_wat_profile) then
             write(*,*) 'Relaxed atmospheric water vapor profile towards ', atm_wat_profile
             write(*,*) 'Unless it reaches saturation (maximal value)'
@@ -821 +817 @@
     endif
 endif
+print*, 'looooooooooool', q(1,1,igcm_h2o_vap)
 
 END SUBROUTINE init_testphys1d