Thermal plume model 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