Difference between revisions of "Thermal plume model Generic PCM"
From Planets
Line 1: | Line 1: | ||
− | == | + | == Overview of the Thermal plume model == |
+ | The Thermal Plume Model (TPM) is a parametrisation of the vertical transport related to the small-scale convective processes. It is based on a mass flux approach, meaning it is a nonlocal scheme which allows mixing between layers that are not directly adjacent. | ||
+ | |||
+ | == How to use it? == | ||
+ | |||
+ | === In a 1D rcm1d run === | ||
+ | |||
+ | The TPM is activated in callphys.def with | ||
+ | calltherm = .True. | ||
+ | |||
+ | Remember to also turn off the convective adjustment scheme to use only one parametrisation of convection | ||
+ | calladj = .True. | ||
+ | |||
+ | Then, several options may be used to tune the model | ||
+ | dvimpl = .False. | ||
+ | dqimpl = .True. | ||
+ | r_aspect_thermals = 2.0 | ||
+ | tau_thermals = 0.0 | ||
+ | betalpha = 0.9 | ||
+ | afact = 0.9 | ||
+ | fact_epsilon = 1.e-4 | ||
+ | alpha_max = 0.7 | ||
+ | fomass_max = 0.5 | ||
+ | pres_limit = 1.e3 | ||
+ | nu = 0.0 | ||
+ | |||
+ | |||
+ | === In a 3D DYNAMICO run === | ||
+ | |||
+ | == Physical description of the TPM == | ||
== Implementation in the Generic PCM == | == Implementation in the Generic PCM == | ||
Line 41: | Line 70: | ||
Implementation to take into account generic tracers | Implementation to take into account generic tracers | ||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
== References == | == References == |
Revision as of 12:03, 5 March 2025
Contents
Overview of the Thermal plume model
The Thermal Plume Model (TPM) is a parametrisation of the vertical transport related to the small-scale convective processes. It is based on a mass flux approach, meaning it is a nonlocal scheme which allows mixing between layers that are not directly adjacent.
How to use it?
In a 1D rcm1d run
The TPM is activated in callphys.def with
calltherm = .True.
Remember to also turn off the convective adjustment scheme to use only one parametrisation of convection
calladj = .True.
Then, several options may be used to tune the model
dvimpl = .False. dqimpl = .True. r_aspect_thermals = 2.0 tau_thermals = 0.0 betalpha = 0.9 afact = 0.9 fact_epsilon = 1.e-4 alpha_max = 0.7 fomass_max = 0.5 pres_limit = 1.e3 nu = 0.0
In a 3D DYNAMICO run
Physical description of the TPM
Implementation in the Generic PCM
Code
Edition history
Authors: F. Hourdin, C. Rio, A. Mathieu, A. Boissinot, A. Le Saux
Version du 09.02.07
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: 1. upward mass flux 2. downward mass flux 3. entrainment 4. a detrainment
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 controled by iflag_thermals = 15, 16 run with impl=-1 : numerical convergence with NPv3 17, 18 run with impl=1 : more stable 15 and 17 correspond to the activation of the stratocumulus "bidouille"
Major changes 2018-19 (AB alexandre.boissinot@lmd.jussieu.fr)
New detr and entre formulae (no longer alimentation) lmin can be greater than 1 Mix every tracer Can stack verticaly multiple plumes (it makes thermcell_dv2 unusable for the moment)
Modif 2024 (ALS, arthur.le-saux@lmd.ipsl.fr)
Implementation to take into account generic tracers