Changeset 3103

Timestamp:

Oct 25, 2023, 1:46:15 PM (13 months ago)

Author:

emillour

Message:

changed tau calculation + comments

File:

• trunk/LMDZ.GENERIC/libf/phystd/condensation_generic_mod.F90 (modified) (2 diffs)

trunk/LMDZ.GENERIC/libf/phystd/condensation_generic_mod.F90

@@ -197 +197 @@
                         if ((perfect_vap_profile) .and. (ngrid.eq.1)) then
                                 ! perfect_vap_profile is a mode that should a priori only be used in 1D:
-                                ! it aims to force the vap profile:
+                                ! as it is written below, it aims to force the vap profile:
                                 ! - below condensation, it is forced to qvap_deep
                                 ! - at condensation levels, it is forced to 99% of sat
                                 ! - above the cold trap, it is forced to the value allowed by the cold trap
-                                tau=3*24*3600
+
+                                ! perfect_vap_profile can be customed as you want
+                                
+                                tau = 10. * ptimestep ! tau must not be lower than the physical step time. 
 
                                 k_cold_trap = 2
@@ -240 +243 @@
 
                         if (qvap_deep >= 0.) then
-                                if (ngrid.eq.1) then ! if 1D
-                                        tau=3*24*3600  ! tau must not be lower than the physical step time. In 1D time step is very long
-                                else
-                                        tau=3600 ! for 3D
-                                endif
-                                !brings lower generic vapor ratio to a fixed value.
+                                tau = 10. * ptimestep ! tau must not be lower than the physical step time. 
+                                ! brings lower generic vapor ratio to a fixed value.
                                 ! tau is in seconds and must not be lower than the physical step time.
                                 pdqvaplsc(1:ngrid,1,igcm_generic_vap) = (qvap_deep - pq(1:ngrid,1,igcm_generic_vap))/tau - pdq(1:ngrid,1,igcm_generic_vap)
                         endif
                         if (qvap_top >= 0.) then
-                                if (ngrid.eq.1) then ! if 1D
-                                        tau=3*24*3600  ! tau must not be lower than the physical step time. In 1D time step is very long
-                                else
-                                        tau=3600 ! for 3D
-                                endif
-                                !brings lower generic vapor ratio to a fixed value.
+                                tau = 10. * ptimestep ! tau must not be lower than the physical step time. 
+                                ! brings lower generic vapor ratio to a fixed value.
                                 ! tau is in seconds and must not be lower than the physical step time.
                                 pdqvaplsc(1:ngrid,nlayer,igcm_generic_vap) = (qvap_top - pq(1:ngrid,nlayer,igcm_generic_vap))/tau - pdq(1:ngrid,nlayer,igcm_generic_vap)