| 1 | ! ***************************************************** |
|---|
| 2 | ! control parameters of the Martian thermal plume model |
|---|
| 3 | ! ***************************************************** |
|---|
| 4 | ! Reference paper: |
|---|
| 5 | ! A. Colaïtis, A. Spiga, F. Hourdin, C. Rio, F. Forget, and E. Millour. |
|---|
| 6 | ! A thermal plume model for the Martian convective boundary layer. |
|---|
| 7 | ! Journal of Geophysical Research (Planets), 118:1468-1487, July 2013. |
|---|
| 8 | ! http://dx.doi.org/10.1002/jgre.20104 |
|---|
| 9 | ! http://arxiv.org/abs/1306.6215 |
|---|
| 10 | ! ----------------------------------------------------------------------- |
|---|
| 11 | ! Author : A. Colaitis 2011-01-05 (with updates 2011-2013) |
|---|
| 12 | ! Institution : Laboratoire de Meteorologie Dynamique (LMD) Paris, France |
|---|
| 13 | ! ----------------------------------------------------------------------- |
|---|
| 14 | ! Corresponding author : A. Spiga aymeric.spiga_AT_upmc.fr |
|---|
| 15 | ! ----------------------------------------------------------------------- |
|---|
| 16 | |
|---|
| 17 | module comtherm_h |
|---|
| 18 | implicit none |
|---|
| 19 | |
|---|
| 20 | !------------------------------------------------------------ |
|---|
| 21 | !------------------------------------------------------------ |
|---|
| 22 | ! nsplit_thermals ! Sub-timestep for the thermal plume model. |
|---|
| 23 | ! Dependent on the timestep chosen for radiative transfer. |
|---|
| 24 | ! It is recommended to run with 96 timestep per day and |
|---|
| 25 | ! call radiative transfer at each timestep, |
|---|
| 26 | ! configuration in which thermals can run |
|---|
| 27 | ! very well with a sub-timestep of 10. |
|---|
| 28 | #ifdef MESOSCALE |
|---|
| 29 | ! ---- MESOSCALE (timesteps < 200s) |
|---|
| 30 | ! -- only a few subtimestep are needed because physical timestep is low |
|---|
| 31 | INTEGER,PARAMETER :: nsplit_thermals = 4 |
|---|
| 32 | #else |
|---|
| 33 | ! ---- GCM |
|---|
| 34 | ! -- we recommend here a value for 96 physical timesteps per day in the GCM |
|---|
| 35 | ! (a larger value will not change results... |
|---|
| 36 | ! but would cost more for negligible benefit in stability and accuracy) |
|---|
| 37 | INTEGER,PARAMETER :: nsplit_thermals = 10 |
|---|
| 38 | ! -- if there is less than 96 physical timesteps per day we recommend |
|---|
| 39 | !INTEGER,PARAMETER :: nsplit_thermals=35 |
|---|
| 40 | #endif |
|---|
| 41 | |
|---|
| 42 | !------------------------------------------------------------ |
|---|
| 43 | !------------------------------------------------------------ |
|---|
| 44 | ! r_aspect_thermals ! Mainly control the shape of the temperature profile |
|---|
| 45 | ! at the bottom of the mixed layer. Decreasing it goes toward |
|---|
| 46 | ! a convective-adjustment like profile. |
|---|
| 47 | ! (see paragraph 45 of paper and appendix S4) |
|---|
| 48 | REAL,PARAMETER :: r_aspect_thermals = 1. |
|---|
| 49 | |
|---|
| 50 | !------------------------------------------------------------ |
|---|
| 51 | !------------------------------------------------------------ |
|---|
| 52 | ! qtransport_thermals ! logical to activate tracer transport in thermals |
|---|
| 53 | ! |
|---|
| 54 | LOGICAL,PARAMETER :: qtransport_thermals = .true. |
|---|
| 55 | |
|---|
| 56 | !------------------------------------------------------------ |
|---|
| 57 | !------------------------------------------------------------ |
|---|
| 58 | ! dtke_thermals ! logical to activate TKE transport in thermals |
|---|
| 59 | ! -- still experimental, for testing purposes only. |
|---|
| 60 | ! -- not used in current thermal plume models both on Earth and Mars. |
|---|
| 61 | LOGICAL,PARAMETER :: dtke_thermals = .false. |
|---|
| 62 | |
|---|
| 63 | !------------------------------------------------------------ |
|---|
| 64 | !------------------------------------------------------------ |
|---|
| 65 | ! thermverbose ! make thermal plume model more verbose |
|---|
| 66 | LOGICAL,PARAMETER :: thermverbose = .false. |
|---|
| 67 | |
|---|
| 68 | |
|---|
| 69 | ! ------------------------------------------------------------------------------------ |
|---|
| 70 | ! -------------- TUNING PARAMETERS FOR MARTIAN THERMALS MODEL ------------------------ |
|---|
| 71 | ! ------------------------------------------------------------------------------------ |
|---|
| 72 | ! Detrainment |
|---|
| 73 | REAL,PARAMETER :: ad = 0.0004 ! D_2 in paper, see paragraph 44 |
|---|
| 74 | REAL,PARAMETER :: bd = -0.6697 ! D_1 in paper, see paragraph 44 |
|---|
| 75 | ! Entrainment |
|---|
| 76 | REAL,PARAMETER :: ae = 0.03683 ! E_1 in paper, see paragraph 43 |
|---|
| 77 | REAL,PARAMETER :: be = 0.631631 ! E_2 in paper, see paragraph 43 |
|---|
| 78 | ! Downdraft |
|---|
| 79 | REAL,PARAMETER :: fdfu=-0.8 ! downdraft to updraft mass flux ratio |
|---|
| 80 | ! see paper paragraph 48 |
|---|
| 81 | REAL,PARAMETER :: omega=-0.03 ! omega. see paper paragraph 48 |
|---|
| 82 | ! Vertical velocity equation |
|---|
| 83 | REAL,PARAMETER :: a1=1. ! a in paper, see paragraph 41 |
|---|
| 84 | REAL,PARAMETER :: b1=0.0001 ! b in paper, see paragraph 41 |
|---|
| 85 | ! Inversion layer (same as a1 b1 actually) |
|---|
| 86 | REAL :: a1inv=1. ! a1 coeff in inversion layer |
|---|
| 87 | REAL :: b1inv=0.0001 ! b1 coeff in inversion layer |
|---|
| 88 | ! ------------------------------------------------------------------------------------ |
|---|
| 89 | ! ------------------------------------------------------------------------------------ |
|---|
| 90 | ! ------------------------------------------------------------------------------------ |
|---|
| 91 | |
|---|
| 92 | end module comtherm_h |
|---|
| 93 | |
|---|
