Changeset 1033 for trunk/LMDZ.MARS/libf/phymars/thermcell_main_mars.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/thermcell_main_mars.F90 (modified) (5 diffs)

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

@@ -1 +1 @@
 !=======================================================================
 ! THERMCELL_MAIN_MARS
-! Author : A. Colaitis after C. Rio and F. Hourdin
-! 
+!=======================================================================
 ! This routine is called by calltherm_interface and is inside a sub-timestep
 ! loop. It computes thermals properties from parametrized entrainment and
@@ -8 +7 @@
 ! Mass flux are then computed and temperature and CO2 MMR are transported.
 !=======================================================================
-! SEE COMMENTS IN CALLTHERM_INTERFACE
+! Author : A. Colaitis 2011-01-05 (with updates 2011-2013)
+!          after C. Rio and F. Hourdin
+! Institution : Laboratoire de Meteorologie Dynamique (LMD) Paris, France
+! -----------------------------------------------------------------------
+! Corresponding author : A. Spiga aymeric.spiga_AT_upmc.fr
+! -----------------------------------------------------------------------
+! ASSOCIATED FILES
+! --> calltherm_interface.F90
+! --> thermcell_dqup.F90
+! --> comtherm_h.F90
 !=======================================================================
-!
-!
+! 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
+! -----------------------------------------------------------------------
+! Reference paper for terrestrial plume model:
+! C. Rio and F. Hourdin.
+! A thermal plume model for the convective boundary layer : Representation of cumulus clouds.
+! Journal of the Atmospheric Sciences, 65:407-425, 2008.
+! -----------------------------------------------------------------------
+
       SUBROUTINE thermcell_main_mars(ngrid,nlayer,nq &
      &                  ,tracer,igcm_co2 &
@@ -124 +143 @@
       REAL zdz,zbuoy(ngrid,nlayer),zw2m
       LOGICAL activecell(ngrid),activetmp(ngrid)
-      REAL a1,b1,ae,be,ad,bd,fdfu,b1inv,a1inv,omega
       INTEGER tic
 
@@ -191 +209 @@
 ! This is especially important for modelling polar convection.
 ! **********************************************************************
-      if (igcm_co2.ne.0) then
+       if (igcm_co2.ne.0) then
 
          m_co2 = 44.01E-3  ! CO2 molecular mass (kg/mol)
@@ -332 +350 @@
       wa_moy(:,:)=0. !vertical velocity
 
-! --------------------------------------------------------------------------
-! -------------- MAIN PARAMETERS FOR THERMALS MODEL ------------------------
-
-! Detrainment
-      ad = 0.0004 !D_2 in paper, see paragraph 44
-      bd = -0.6697 !D_1 in paper, see paragraph 44
-
-! Entrainment
-      ae = 0.03683 !E_1 in paper, see paragraph 43
-      be = 0.631631 !E_2 in paper, see paragraph 43
-
-! Downdraft
-      fdfu=-0.8 !downdraft to updraft mass flux ratio, see paper paragraph 48
-      omega=-0.03 !see paper paragraph 48
-
-! Vertical velocity equation
-      a1=1. !a in paper, see paragraph 41
-      b1=0.0001 !b in paper, see paragraph 41
-
-! Inversion layer
-      a1inv=a1 !a1 in inversion layer
-      b1inv=b1 !b1 in inversion layer
-! --------------------------------------------------------------------------
-! --------------------------------------------------------------------------
-! --------------------------------------------------------------------------
-
 ! Some more initializations
       wmaxa(:)=0.