Planets - User contributions [en]

Thermal plume model Generic PCM

2025-04-30T09:18:36Z

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

== 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 (as it default value is .True.), to use only one parametrisation of convection
calladj = .False.

Then, several options may be used to tune the model
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 ! 1.0 Aspect ratio of the thermals (width / height)
tau_thermals = 0.0 ! 0. Relaxation time
betalpha = 0.9 ! 0.9 - between 0 (e=d=0) and 1 (rho*fraca=cst)
afact = 0.9 ! 2./3. - buoyancy acceleration efficiency, between 0 and 1
fact_epsilon = 1.e-4 ! 2.e-3 - turbulence and friction deceleration
alpha_max = 0.7 ! 0.7 Maximal permitted updraft fraction
fomass_max = 0.5 ! 0.5 Maximal permitted outgoing layer mass fraction
pres_limit = 1.e3 ! 1.e5 Minimum pressure level above which thermals cannot be triggered
nu = 0.0 ! 0. Geometrical contributions to entrainment and detrainment
linf = 1.0 ! 1.0 linf is used to set the lowest possible first level because we allow it
! to begin higher than the surface. It is set to 2 in order to remove the
! first layer for gas giant.
! If there is a surface, it has to be set to 1.
! If someone want to call more than once the thermal plume model in some
! grid points, this variable may become a saved array of INTEGER with size ngrid.

=== 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

== References ==

[[Category:Generic-Model]]
[[Category:Generic-DYNAMICO]]

Thermal plume model Generic PCM

2025-04-30T09:09:52Z

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

== 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 (as it default value is .True.), to use only one parametrisation of convection
calladj = .False.

Then, several options may be used to tune the model
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
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

== References ==

[[Category:Generic-Model]]
[[Category:Generic-DYNAMICO]]

Thermal plume model Generic PCM

2025-04-30T09:05:17Z

Alesaux:

== 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 (as it default value is .True.), to use only one parametrisation of convection
calladj = .False.

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.
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

== References ==

[[Category:Generic-Model]]
[[Category:Generic-DYNAMICO]]

Physics of the Generic PCM

2025-03-05T10:04:31Z

Alesaux:

This page describes the various physical parametrizations of the Generic PCM and the chronology of their call through the physical iteration. This chronology is important because some variables need to be updated by certain processes before others (examples).
During one physical iteration, the code passes through multiple sub-routines each encapsulating a parametrization. The sub-routines usually take as arguments:
*the dynamical values of the state variables
*the dynamical tendencies of the state variables
*any additional relevant variable
*any additional relevant tendency
and it returns in general tendencies (of the state variables as well as of any other relevant variable).

Work in progress. Need to add links.

=Initialization=
==First call==
Some initializations only need to be done at the very first iteration (e.g. examples).
==Each call==
Some other initializations need to be done at each iteration (e.g. examples).

=Radiative transfer=
==Correlated-k==
The main radiative transfer solver of the Generic PCM implements the correlated-k method, which provides a flexible and quick way to solve radiative transfer equations, particularly suited for GCMs. More informations here.
==Newtonian relaxation==
Newton was a very relaxed physicist, still inspiring us today.
==No atmosphere==
If you have no atmosphere, why do you need a GCM?

=Vertical diffusion=
Compute the vertical diffusion due to turbulence in the planetary boundary layer
==vdifc==
The "old" vertical diffusion routine.
==turbdiff==
The new (and improved!) vertical diffusion routine

=Convection=
Convective mixing in an atmosphere column involves non-resolved, sub-grid processes. These processes are parameterized using the following routines:

==Thermal plume==
This module implements a more accurate and sophisticated parametrization of convection. More info [[Thermal_plume_model_Generic_PCM|here]].

==Dry convection==
More info [[Convective_adjustment_scheme_in_the_generic_PCM|here]].
==Non-orographic gravity waves==
More info here.

=CO2 condensation=
Inheriting from the Mars PCM where it is the background species, CO2 condensation is treated as a dedicated step in the Generic PCM's physics. Does that only concern CO2 as a background gas? Or does it work also if CO2 is a non-background variable gas? In any case more info here.

=Tracers=
Many things can be advected by tracers in the Generic PCM, like chemical species or aerosols. Physical processes involving tracers are parameterized using the following routines:

==Volcano==
This routine parameterizes a source of tracers corresponding to volcanic eruption. More information here.

==Water/ice==
Water aerosols (liquid or solid) are created (resp. consumed) by condensation (resp. vaporization or sublimation), consuming (resp. releasing) latent heat in the atmosphere. In the Generic PCM, the (atmospheric part of the) water cycle is handled by various routines, as explained here.

==Photochemistry==
Chemistry can turn molecules into other molecules, by the action of temperature (thermochemistry) or UV light (photochemistry). This is handled by the photochemistry routine, described [[Photochemistry|here]].

==Generic condensation==
On Earth, only water condenses in the atmosphere, but on other planets (which the Generic PCM aims at simulating), many other chemicals can undergo state change. To take that into account, the Generic PCM has a flexible scheme to deal with any arbitrary species, as desbribed [[Radiative_Generic_Condensable_Specie|here]].

==Sedimentation==
What goes up must come down, as explained [[Sedimentation_of_tracers_in_the_generic_PCM|here]].

==Updates==
This section essentially takes care that the condensation of a major species (e.g. water vapor in steam-rich atmospheres) affects other species, as described here.

==Slab ocean==
This routine solves for big fish eating small fish, as desbried [[Slab_ocean_model|here]].

==Surface==
The surface part of the water cycle is handled here, as the hydrology page describes.

=Subsurface thermal=
Thermal diffusion in the subsurface is solved here, as described there.

=Diagnostics/write outputs=
Nothing very physical here, just writting the outputs! By the way, if you want to know how to output new variables in the diagfi.nc file, check out [[Outputs|this]] page.

Thermal plume model Generic PCM

2025-03-05T10:03:59Z

Alesaux:

Thermal plume model Generic PCM

2025-03-05T10:03:23Z

Alesaux:

Thermal plume model Generic PCM

2025-03-05T09:54:05Z

Alesaux:

== Description of the Thermal plume model ==

== 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

== 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 ===

== References ==

[[Category:Generic-Model]]
[[Category:Generic-DYNAMICO]]

Thermal plume model Generic PCM

2025-03-05T09:52:52Z

Alesaux:

== Description of the Thermal plume model ==

== 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

== 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 ===

== References ==

[[Category:Generic-Model]]
[[Category:Generic-DYNAMICO]]

Thermal plume model Generic PCM

2025-03-05T09:33:05Z

Alesaux: /* How to use it? */

== Description of the Thermal plume model ==

== Implementation in the Generic PCM ==

== How to use it? ==

=== In a 1D rcm1d run ===

The TPM is activated in callphys.def with
calltherm = .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 ===

== References ==

Thermal plume model Generic PCM

2025-03-05T09:12:26Z

Alesaux: /* How to use it? */

== Description of the Thermal plume model ==

== Implementation in the Generic PCM ==

== How to use it? ==

The TPM is activated in callphys.def with:
calltherm = .true.

== References ==

Thermal plume model Generic PCM

2025-03-05T09:10:48Z

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

== Description of the Thermal plume model ==

== Implementation in the Generic PCM ==

== How to use it? ==

== References ==