Changeset 3103 for trunk/LMDZ.GENERIC


Ignore:
Timestamp:
Oct 25, 2023, 1:46:15 PM (14 months ago)
Author:
emillour
Message:

changed tau calculation + comments

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/LMDZ.GENERIC/libf/phystd/condensation_generic_mod.F90

    r3044 r3103  
    197197                        if ((perfect_vap_profile) .and. (ngrid.eq.1)) then
    198198                                ! perfect_vap_profile is a mode that should a priori only be used in 1D:
    199                                 ! it aims to force the vap profile:
     199                                ! as it is written below, it aims to force the vap profile:
    200200                                ! - below condensation, it is forced to qvap_deep
    201201                                ! - at condensation levels, it is forced to 99% of sat
    202202                                ! - above the cold trap, it is forced to the value allowed by the cold trap
    203                                 tau=3*24*3600
     203
     204                                ! perfect_vap_profile can be customed as you want
     205                               
     206                                tau = 10. * ptimestep ! tau must not be lower than the physical step time.
    204207
    205208                                k_cold_trap = 2
     
    240243
    241244                        if (qvap_deep >= 0.) then
    242                                 if (ngrid.eq.1) then ! if 1D
    243                                         tau=3*24*3600  ! tau must not be lower than the physical step time. In 1D time step is very long
    244                                 else
    245                                         tau=3600 ! for 3D
    246                                 endif
    247                                 !brings lower generic vapor ratio to a fixed value.
     245                                tau = 10. * ptimestep ! tau must not be lower than the physical step time.
     246                                ! brings lower generic vapor ratio to a fixed value.
    248247                                ! tau is in seconds and must not be lower than the physical step time.
    249248                                pdqvaplsc(1:ngrid,1,igcm_generic_vap) = (qvap_deep - pq(1:ngrid,1,igcm_generic_vap))/tau - pdq(1:ngrid,1,igcm_generic_vap)
    250249                        endif
    251250                        if (qvap_top >= 0.) then
    252                                 if (ngrid.eq.1) then ! if 1D
    253                                         tau=3*24*3600  ! tau must not be lower than the physical step time. In 1D time step is very long
    254                                 else
    255                                         tau=3600 ! for 3D
    256                                 endif
    257                                 !brings lower generic vapor ratio to a fixed value.
     251                                tau = 10. * ptimestep ! tau must not be lower than the physical step time.
     252                                ! brings lower generic vapor ratio to a fixed value.
    258253                                ! tau is in seconds and must not be lower than the physical step time.
    259254                                pdqvaplsc(1:ngrid,nlayer,igcm_generic_vap) = (qvap_top - pq(1:ngrid,nlayer,igcm_generic_vap))/tau - pdq(1:ngrid,nlayer,igcm_generic_vap)
Note: See TracChangeset for help on using the changeset viewer.