Thermal plume model Generic PCM - Revision history

Alesaux: /* In a 1D rcm1d run */

2025-04-30T09:18:36Z

‎In a 1D rcm1d run

Alesaux: /* In a 1D rcm1d run */

2025-04-30T09:09:52Z

‎In a 1D rcm1d run

Alesaux at 09:05, 30 April 2025

2025-04-30T09:05:17Z

Alesaux at 10:03, 5 March 2025

2025-03-05T10:03:59Z

Alesaux at 10:03, 5 March 2025

2025-03-05T10:03:23Z

Alesaux at 09:54, 5 March 2025

2025-03-05T09:54:05Z

Alesaux at 09:52, 5 March 2025

2025-03-05T09:52:52Z

Alesaux: /* How to use it? */

2025-03-05T09:33:05Z

‎How to use it?

Alesaux: /* How to use it? */

2025-03-05T09:12:26Z

‎How to use it?

Alesaux: Created page with " == Description of the Thermal plume model == == Implementation in the Generic PCM == == How to use it? == == References =="

2025-03-05T09:10:48Z

Created page with " == Description of the Thermal plume model == == Implementation in the Generic PCM == == How to use it? == == References =="

New page

== Description of the Thermal plume model ==

== Implementation in the Generic PCM ==

== How to use it? ==

== References ==

@@ Line 17: / Line 17: @@
   dqimpl               = .True.   !  .true.   Flag for thermcell_dq version (True : implicit scheme || False : explicit scheme)
   dvimpl               = .False.  !  .false.  Flag for specific u, v mixing (True : thermcell_dv2 || False : thermcell_dq)
-  r_aspect_thermals    = 2.0
+  r_aspect_thermals    = 2.0      !  1.0      Aspect ratio of the thermals (width / height)
-  tau_thermals         = 0.0
+  tau_thermals         = 0.0      !  0.       Relaxation time
-  betalpha             = 0.9
+  betalpha             = 0.9      !  0.9      - between 0 (e=d=0) and 1 (rho*fraca=cst)
-  afact                = 0.9
+  afact                = 0.9      !  2./3.    - buoyancy acceleration efficiency, between 0 and 1
-  fact_epsilon         = 1.e-4
+  fact_epsilon         = 1.e-4    !  2.e-3    - turbulence and friction deceleration
-  alpha_max            = 0.7
+  alpha_max            = 0.7      !  0.7      Maximal permitted updraft fraction
-  fomass_max           = 0.5
+  fomass_max           = 0.5      !  0.5      Maximal permitted outgoing layer mass fraction
-  pres_limit           = 1.e3
+  pres_limit           = 1.e3     !  1.e5     Minimum pressure level above which thermals cannot be triggered
-  nu                   = 0.0
+  nu                   = 0.0      !  0.       Geometrical contributions to entrainment and detrainment
 === In a 3D DYNAMICO run ===

@@ Line 15: / Line 15: @@
 Then, several options may be used to tune the model
-  dvimpl               = .False.   ! Implicit (.True.) or explicit (.False) computation of vertical advection in the environment of thermal plume
+  dqimpl               = .True.   !  .true.   Flag for thermcell_dq version (True : implicit scheme || False : explicit scheme)
-  dqimpl               = .True.
+  dvimpl               = .False.  !  .false.  Flag for specific u, v mixing (True : thermcell_dv2 || False : thermcell_dq)
   r_aspect_thermals    = 2.0
   tau_thermals         = 0.0
@@ Line 26: / Line 26: @@
   pres_limit           = 1.e3
   nu                   = 0.0
 === In a 3D DYNAMICO run ===

@@ Line 15: / Line 15: @@
 Then, several options may be used to tune the model
-  dvimpl               = .False.
+  dvimpl               = .False.   ! Implicit (.True.) or explicit (.False) computation of vertical advection in the environment of thermal plume
   dqimpl               = .True.
   r_aspect_thermals    = 2.0

@@ Line 11: / Line 11: @@
   calltherm = .True.
-Remember to also turn off the convective adjustment scheme to use only one parametrisation of convection
+Remember to also turn off the convective adjustment scheme (as it default value is .True.), to use only one parametrisation of convection
-  calladj  = .True.
+  calladj  = .False.
 Then, several options may be used to tune the model

@@ Line 1: / Line 1: @@
-== Description of the Thermal plume model ==
+== Overview of the Thermal plume model ==
 == Implementation in the Generic PCM ==
@@ Line 41: / Line 70: @@
      Implementation to take into account generic tracers
 == References ==

@@ Line 9: / Line 9: @@
 === Edition history ===
-   Authors: F. Hourdin, C. Rio, A. Mathieu, A. Boissinot, A. Le Saux
+Authors: F. Hourdin, C. Rio, A. Mathieu, A. Boissinot, A. Le Saux
     Calculation of vertical transport in the boundary layer in the presence of explicitly represented “thermals” with cloud processes
     Rewriting from a paper listing in Habas, 14/02/00
     Thermal energy is assumed to be homogeneous and dissipated by mixing with its surroundings.
     The length l_mix controls the mixing efficiency.
     The transport of the different species is calculated by taking into account:
 . upward mass flux
@@ Line 25: / Line 23: @@
 . a detrainment
- Modif 2013/01/04 (FH hourdin@lmd.jussieu.fr)
+Modif 2013/01/04 (FH hourdin@lmd.jussieu.fr)
      Introduction of an implicit computation of vertical advection in the environment of thermal plumes in thermcell_dq
      impl = 0 : explicit ; impl = 1 : implicit ; impl =-1 : old version
@@ Line 33: / Line 32: @@
 and 17 correspond to the activation of the stratocumulus "bidouille"
- Major changes 2018-19 (AB alexandre.boissinot@lmd.jussieu.fr)
+Major changes 2018-19 (AB alexandre.boissinot@lmd.jussieu.fr)
      New detr and entre formulae (no longer alimentation)
      lmin can be greater than 1
@@ Line 39: / Line 38: @@
      Can stack verticaly multiple plumes (it makes thermcell_dv2 unusable for the moment)
- Modif 2024 (ALS, arthur.le-saux@lmd.ipsl.fr)
+Modif 2024 (ALS, arthur.le-saux@lmd.ipsl.fr)
      Implementation to take into account generic tracers

@@ Line 7: / Line 7: @@
 == How to use it? ==
-The TPM is activated in callphys.def with:
+=== In a 1D rcm1d run ===
   calltherm = .true.
 == References ==