| 1 | module comtherm_h |
|---|
| 2 | implicit none |
|---|
| 3 | |
|---|
| 4 | ! ******************************************* |
|---|
| 5 | ! control parameters of the thermal model |
|---|
| 6 | ! SEE COMMENTS IN CALLTHERM_INTERFACE |
|---|
| 7 | ! ******************************************* |
|---|
| 8 | |
|---|
| 9 | !------------------------------------------------------------ |
|---|
| 10 | !------------------------------------------------------------ |
|---|
| 11 | ! thermverbose ! make thermal plume model more verbose |
|---|
| 12 | LOGICAL,PARAMETER :: thermverbose = .false. |
|---|
| 13 | |
|---|
| 14 | !------------------------------------------------------------ |
|---|
| 15 | !------------------------------------------------------------ |
|---|
| 16 | ! qtransport_thermals ! logical to activate tracer transport in thermals |
|---|
| 17 | ! |
|---|
| 18 | LOGICAL,PARAMETER :: qtransport_thermals = .true. |
|---|
| 19 | |
|---|
| 20 | !------------------------------------------------------------ |
|---|
| 21 | !------------------------------------------------------------ |
|---|
| 22 | ! dtke_thermals ! logical to activate TKE in thermals |
|---|
| 23 | ! -- still experimental, for testing purposes only. |
|---|
| 24 | ! -- not used in current thermal plume models both on Earth and Mars. |
|---|
| 25 | LOGICAL,PARAMETER :: dtke_thermals = .false. |
|---|
| 26 | |
|---|
| 27 | !------------------------------------------------------------ |
|---|
| 28 | !------------------------------------------------------------ |
|---|
| 29 | ! r_aspect_thermals ! Mainly control the shape of the temperature profile |
|---|
| 30 | ! in the surface layer. Decreasing it goes toward |
|---|
| 31 | ! a convective-adjustment like profile. |
|---|
| 32 | ! (see paragraph 45 of paper and appendix S4) |
|---|
| 33 | REAL,PARAMETER :: r_aspect_thermals = 1. |
|---|
| 34 | |
|---|
| 35 | !------------------------------------------------------------ |
|---|
| 36 | !------------------------------------------------------------ |
|---|
| 37 | ! nsplit_thermals ! Sub-timestep for the thermal plume model. |
|---|
| 38 | ! Dependent on the timestep chosen for radiative transfer. |
|---|
| 39 | ! It is recommended to run with 96 timestep per day and |
|---|
| 40 | ! call radiative transfer at each timestep, |
|---|
| 41 | ! configuration in which thermals can run |
|---|
| 42 | ! very well with a sub-timestep of 10. |
|---|
| 43 | #ifdef MESOSCALE |
|---|
| 44 | ! ---- MESOSCALE (timesteps < 200s) |
|---|
| 45 | ! -- only a few subtimestep are needed because physical timestep is low |
|---|
| 46 | INTEGER,PARAMETER :: nsplit_thermals = 4 |
|---|
| 47 | #else |
|---|
| 48 | ! ---- GCM |
|---|
| 49 | ! -- we recommend here a value for 96 physical timesteps per day in the GCM |
|---|
| 50 | ! (a larger value will not change results... |
|---|
| 51 | ! but would cost more for negligible benefit in stability and accuracy) |
|---|
| 52 | INTEGER,PARAMETER :: nsplit_thermals = 10 |
|---|
| 53 | ! -- if there is less than 96 physical timesteps per day we recommend |
|---|
| 54 | !INTEGER,PARAMETER :: nsplit_thermals=35 |
|---|
| 55 | #endif |
|---|
| 56 | |
|---|
| 57 | end module comtherm_h |
|---|
| 58 | |
|---|
