Changeset 1033 for trunk/LMDZ.MARS/libf/phymars/comtherm_h.F90

Timestamp:

Sep 8, 2013, 11:43:08 PM (11 years ago)

Author:

aslmd

Message:

LMDZ.MARS additional comments of thermal plume model. moved tuning variables in comtherm_h

File:

• trunk/LMDZ.MARS/libf/phymars/comtherm_h.F90 (modified) (2 diffs)

trunk/LMDZ.MARS/libf/phymars/comtherm_h.F90

@@ +1 @@
+! *****************************************************
+! control parameters of the Martian thermal plume model
+! *****************************************************
+! Reference paper:
+! A. Colaïtis, A. Spiga, F. Hourdin, C. Rio, F. Forget, and E. Millour. 
+! A thermal plume model for the Martian convective boundary layer. 
+! Journal of Geophysical Research (Planets), 118:1468-1487, July 2013.
+! http://dx.doi.org/10.1002/jgre.20104
+! http://arxiv.org/abs/1306.6215
+! -----------------------------------------------------------------------
+! Author : A. Colaitis 2011-01-05 (with updates 2011-2013)
+! Institution : Laboratoire de Meteorologie Dynamique (LMD) Paris, France
+! -----------------------------------------------------------------------
+! Corresponding author : A. Spiga aymeric.spiga_AT_upmc.fr
+! -----------------------------------------------------------------------
+
        module comtherm_h
        implicit none
-
-       ! *******************************************
-       ! control parameters of the thermal model
-       ! SEE COMMENTS IN CALLTHERM_INTERFACE
-       ! *******************************************
-
-       !------------------------------------------------------------
-       !------------------------------------------------------------ 
-       ! thermverbose  ! make thermal plume model more verbose
-       LOGICAL,PARAMETER :: thermverbose = .false.
-
-       !------------------------------------------------------------
-       !------------------------------------------------------------ 
-       ! qtransport_thermals  ! logical to activate tracer transport in thermals
-       !
-       LOGICAL,PARAMETER :: qtransport_thermals = .true. 
-
-       !------------------------------------------------------------
-       !------------------------------------------------------------ 
-       ! dtke_thermals  ! logical to activate TKE in thermals
-                        ! -- still experimental, for testing purposes only.
-                        ! -- not used in current thermal plume models both on Earth and Mars.
-       LOGICAL,PARAMETER :: dtke_thermals = .false.
-
-       !------------------------------------------------------------
-       !------------------------------------------------------------ 
-       ! r_aspect_thermals     ! Mainly control the shape of the temperature profile
-                               ! in the surface layer. Decreasing it goes toward
-                               ! a convective-adjustment like profile.
-                               ! (see paragraph 45 of paper and appendix S4)
-       REAL,PARAMETER :: r_aspect_thermals = 1.
 
        !------------------------------------------------------------
@@ -55 +40 @@
 #endif
 
+       !------------------------------------------------------------
+       !------------------------------------------------------------ 
+       ! r_aspect_thermals     ! Mainly control the shape of the temperature profile
+                               ! at the bottom of the mixed layer. Decreasing it goes toward
+                               ! a convective-adjustment like profile.
+                               ! (see paragraph 45 of paper and appendix S4)
+       REAL,PARAMETER :: r_aspect_thermals = 1.
+
+       !------------------------------------------------------------
+       !------------------------------------------------------------ 
+       ! qtransport_thermals  ! logical to activate tracer transport in thermals
+       !
+       LOGICAL,PARAMETER :: qtransport_thermals = .true. 
+
+       !------------------------------------------------------------
+       !------------------------------------------------------------ 
+       ! dtke_thermals  ! logical to activate TKE transport in thermals
+                        ! -- still experimental, for testing purposes only.
+                        ! -- not used in current thermal plume models both on Earth and Mars.
+       LOGICAL,PARAMETER :: dtke_thermals = .false.
+
+       !------------------------------------------------------------
+       !------------------------------------------------------------ 
+       ! thermverbose  ! make thermal plume model more verbose
+       LOGICAL,PARAMETER :: thermverbose = .false.
+
+
+       ! ------------------------------------------------------------------------------------
+       ! -------------- TUNING PARAMETERS FOR MARTIAN THERMALS MODEL ------------------------
+       ! ------------------------------------------------------------------------------------
+       ! Detrainment
+       REAL,PARAMETER :: ad = 0.0004    ! D_2 in paper, see paragraph 44
+       REAL,PARAMETER :: bd = -0.6697   ! D_1 in paper, see paragraph 44
+       ! Entrainment
+       REAL,PARAMETER :: ae = 0.03683   ! E_1 in paper, see paragraph 43
+       REAL,PARAMETER :: be = 0.631631  ! E_2 in paper, see paragraph 43
+       ! Downdraft
+       REAL,PARAMETER :: fdfu=-0.8      ! downdraft to updraft mass flux ratio
+                                        ! see paper paragraph 48
+       REAL,PARAMETER :: omega=-0.03    ! omega. see paper paragraph 48
+       ! Vertical velocity equation
+       REAL,PARAMETER :: a1=1.          ! a in paper, see paragraph 41
+       REAL,PARAMETER :: b1=0.0001      ! b in paper, see paragraph 41
+       ! Inversion layer
+       REAL :: a1inv ! a1 coeff in inversion layer
+       REAL :: b1inv ! b1 coeff in inversion layer
+       a1inv = a1  
+       b1inv = b1  
+       ! ------------------------------------------------------------------------------------
+       ! ------------------------------------------------------------------------------------
+       ! ------------------------------------------------------------------------------------
+
        end module comtherm_h