Changeset 2144 for trunk/LMDZ.GENERIC/libf/phystd/thermcell_mod.F90

Timestamp:

Jun 20, 2019, 7:07:28 PM (5 years ago)

Author:

aboissinot

Message:

Thermal plume model parameters are now set in callphys.def and (re)add flag to choose between thermcell_dv2 or thermcell_dq to mix horizontal momentum.

File:

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

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

@@ -6 +6 @@
 ! Flags for computations
                                                  !  default
-      INTEGER,SAVE :: dqimpl                     !  1        flag for thermcell_dq version (1 : implicit scheme || 0 : explicit scheme)
+      LOGICAL,SAVE :: dqimpl                     !  .true.   Flag for thermcell_dq version (True : implicit scheme || False : explicit scheme)
+      LOGICAL,SAVE :: dvimpl                     !  .false.  Flag for specific u, v mixing (True : thermcell_dv2 || False : thermcell_dq)
       
 ! Physical parameters
@@ -69 +70 @@
          RD = r
          
+         print *, 'Implicit mixing scheme ?'
+         dqimpl = .true.
+         call getin_p('dqimpl', dqimpl)
+         print *, 'dqimpl = ', dqimpl
+         
+         print *, 'Use specific horizontal momentum scheme ?'
+         dvimpl = .false.
+         call getin_p('dvimpl', dvimpl)
+         print *, 'dvimpl = ', dvimpl
+         
          print *, 'Plume aspect ratio ?'
          r_aspect_thermals = 1.