Changeset 1033 for trunk/LMDZ.MARS/libf
- Timestamp:
- Sep 8, 2013, 11:43:08 PM (12 years ago)
- Location:
- trunk/LMDZ.MARS/libf/phymars
- Files:
-
- 4 edited
-
calltherm_interface.F90 (modified) (4 diffs)
-
comtherm_h.F90 (modified) (2 diffs)
-
thermcell_dqup.F90 (modified) (2 diffs)
-
thermcell_main_mars.F90 (modified) (5 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.MARS/libf/phymars/calltherm_interface.F90
r1032 r1033 1 ! -----------------------------------------------------------------------1 !======================================================================= 2 2 ! CALLTHERM_INTERFACE 3 ! -----------------------------------------------------------------------3 !======================================================================= 4 4 ! Main interface to the Martian thermal plume model 5 5 ! This interface handles sub-timesteps for this model … … 12 12 ! radiative transfer and (orographic) gravity wave drag. 13 13 ! ----------------------------------------------------------------------- 14 ! Author : A. Colaitis 2011-01-05 (with updates 2011-2013) 15 ! after C. Rio and F. Hourdin 16 ! Institution : Laboratoire de Meteorologie Dynamique (LMD) Paris, France 17 ! ----------------------------------------------------------------------- 18 ! Corresponding author : A. Spiga aymeric.spiga_AT_upmc.fr 19 ! ----------------------------------------------------------------------- 14 20 ! ASSOCIATED FILES 21 ! --> thermcell_main_mars.F90 15 22 ! --> thermcell_dqup.F90 16 ! --> thermcell_main_mars.F9017 23 ! --> comtherm_h.F90 18 24 ! ----------------------------------------------------------------------- … … 21 27 ! A thermal plume model for the Martian convective boundary layer. 22 28 ! Journal of Geophysical Research (Planets), 118:1468-1487, July 2013. 29 ! http://dx.doi.org/10.1002/jgre.20104 30 ! http://arxiv.org/abs/1306.6215 23 31 ! ----------------------------------------------------------------------- 24 32 ! Reference paper for terrestrial plume model: … … 26 34 ! A thermal plume model for the convective boundary layer : Representation of cumulus clouds. 27 35 ! Journal of the Atmospheric Sciences, 65:407-425, 2008. 28 ! -----------------------------------------------------------------------29 ! Author : A. Colaitis 2011-01-05 (with updates 2011-2013)30 ! Institution : Laboratoire de Meteorologie Dynamique (LMD) Paris, France31 ! Corresponding author : A. Spiga aymeric.spiga_AT_upmc.fr32 36 ! ----------------------------------------------------------------------- 33 37 -
trunk/LMDZ.MARS/libf/phymars/comtherm_h.F90
r1032 r1033 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 1 17 module comtherm_h 2 18 implicit none 3 4 ! *******************************************5 ! control parameters of the thermal model6 ! SEE COMMENTS IN CALLTHERM_INTERFACE7 ! *******************************************8 9 !------------------------------------------------------------10 !------------------------------------------------------------11 ! thermverbose ! make thermal plume model more verbose12 LOGICAL,PARAMETER :: thermverbose = .false.13 14 !------------------------------------------------------------15 !------------------------------------------------------------16 ! qtransport_thermals ! logical to activate tracer transport in thermals17 !18 LOGICAL,PARAMETER :: qtransport_thermals = .true.19 20 !------------------------------------------------------------21 !------------------------------------------------------------22 ! dtke_thermals ! logical to activate TKE in thermals23 ! -- 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 profile30 ! in the surface layer. Decreasing it goes toward31 ! a convective-adjustment like profile.32 ! (see paragraph 45 of paper and appendix S4)33 REAL,PARAMETER :: r_aspect_thermals = 1.34 19 35 20 !------------------------------------------------------------ … … 55 40 #endif 56 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 86 REAL :: a1inv ! a1 coeff in inversion layer 87 REAL :: b1inv ! b1 coeff in inversion layer 88 a1inv = a1 89 b1inv = b1 90 ! ------------------------------------------------------------------------------------ 91 ! ------------------------------------------------------------------------------------ 92 ! ------------------------------------------------------------------------------------ 93 57 94 end module comtherm_h 58 95 -
trunk/LMDZ.MARS/libf/phymars/thermcell_dqup.F90
r1032 r1033 1 subroutine thermcell_dqup(ngrid,nlayer,ptimestep,fm,entr,detr, & 2 & masse0,q_therm,dq_therm,ndt,zlmax) 3 implicit none 4 1 !======================================================================= 2 ! THERMCELL_DQUP 5 3 !======================================================================= 6 4 ! … … 13 11 ! 14 12 !======================================================================= 15 ! Arnaud Colaitis 2011-01-05 (modified 2011-2013) 16 ! SEE COMMENTS IN CALLTHERM_INTERFACE 17 !======================================================================= 13 ! Author : A. Colaitis 2011-01-05 (with updates 2011-2013) 14 ! Institution : Laboratoire de Meteorologie Dynamique (LMD) Paris, France 15 ! ----------------------------------------------------------------------- 16 ! Corresponding author : A. Spiga aymeric.spiga_AT_upmc.fr 17 ! ----------------------------------------------------------------------- 18 ! Reference paper: 19 ! A. Colaïtis, A. Spiga, F. Hourdin, C. Rio, F. Forget, and E. Millour. 20 ! A thermal plume model for the Martian convective boundary layer. 21 ! Journal of Geophysical Research (Planets), 118:1468-1487, July 2013. 22 ! http://dx.doi.org/10.1002/jgre.20104 23 ! http://arxiv.org/abs/1306.6215 24 ! ----------------------------------------------------------------------- 25 26 subroutine thermcell_dqup(ngrid,nlayer,ptimestep,fm,entr,detr, & 27 & masse0,q_therm,dq_therm,ndt,zlmax) 28 implicit none 18 29 19 30 ! ============================ INPUTS ============================ -
trunk/LMDZ.MARS/libf/phymars/thermcell_main_mars.F90
r1032 r1033 1 1 !======================================================================= 2 2 ! THERMCELL_MAIN_MARS 3 ! Author : A. Colaitis after C. Rio and F. Hourdin 4 ! 3 !======================================================================= 5 4 ! This routine is called by calltherm_interface and is inside a sub-timestep 6 5 ! loop. It computes thermals properties from parametrized entrainment and … … 8 7 ! Mass flux are then computed and temperature and CO2 MMR are transported. 9 8 !======================================================================= 10 ! SEE COMMENTS IN CALLTHERM_INTERFACE 9 ! Author : A. Colaitis 2011-01-05 (with updates 2011-2013) 10 ! after C. Rio and F. Hourdin 11 ! Institution : Laboratoire de Meteorologie Dynamique (LMD) Paris, France 12 ! ----------------------------------------------------------------------- 13 ! Corresponding author : A. Spiga aymeric.spiga_AT_upmc.fr 14 ! ----------------------------------------------------------------------- 15 ! ASSOCIATED FILES 16 ! --> calltherm_interface.F90 17 ! --> thermcell_dqup.F90 18 ! --> comtherm_h.F90 11 19 !======================================================================= 12 ! 13 ! 20 ! Reference paper: 21 ! A. Colaïtis, A. Spiga, F. Hourdin, C. Rio, F. Forget, and E. Millour. 22 ! A thermal plume model for the Martian convective boundary layer. 23 ! Journal of Geophysical Research (Planets), 118:1468-1487, July 2013. 24 ! http://dx.doi.org/10.1002/jgre.20104 25 ! http://arxiv.org/abs/1306.6215 26 ! ----------------------------------------------------------------------- 27 ! Reference paper for terrestrial plume model: 28 ! C. Rio and F. Hourdin. 29 ! A thermal plume model for the convective boundary layer : Representation of cumulus clouds. 30 ! Journal of the Atmospheric Sciences, 65:407-425, 2008. 31 ! ----------------------------------------------------------------------- 32 14 33 SUBROUTINE thermcell_main_mars(ngrid,nlayer,nq & 15 34 & ,tracer,igcm_co2 & … … 124 143 REAL zdz,zbuoy(ngrid,nlayer),zw2m 125 144 LOGICAL activecell(ngrid),activetmp(ngrid) 126 REAL a1,b1,ae,be,ad,bd,fdfu,b1inv,a1inv,omega127 145 INTEGER tic 128 146 … … 191 209 ! This is especially important for modelling polar convection. 192 210 ! ********************************************************************** 193 if (igcm_co2.ne.0) then211 if (igcm_co2.ne.0) then 194 212 195 213 m_co2 = 44.01E-3 ! CO2 molecular mass (kg/mol) … … 332 350 wa_moy(:,:)=0. !vertical velocity 333 351 334 ! --------------------------------------------------------------------------335 ! -------------- MAIN PARAMETERS FOR THERMALS MODEL ------------------------336 337 ! Detrainment338 ad = 0.0004 !D_2 in paper, see paragraph 44339 bd = -0.6697 !D_1 in paper, see paragraph 44340 341 ! Entrainment342 ae = 0.03683 !E_1 in paper, see paragraph 43343 be = 0.631631 !E_2 in paper, see paragraph 43344 345 ! Downdraft346 fdfu=-0.8 !downdraft to updraft mass flux ratio, see paper paragraph 48347 omega=-0.03 !see paper paragraph 48348 349 ! Vertical velocity equation350 a1=1. !a in paper, see paragraph 41351 b1=0.0001 !b in paper, see paragraph 41352 353 ! Inversion layer354 a1inv=a1 !a1 in inversion layer355 b1inv=b1 !b1 in inversion layer356 ! --------------------------------------------------------------------------357 ! --------------------------------------------------------------------------358 ! --------------------------------------------------------------------------359 360 352 ! Some more initializations 361 353 wmaxa(:)=0.
Note: See TracChangeset
for help on using the changeset viewer.
