Changeset 2231 for trunk/LMDZ.GENERIC/libf/phystd/thermcell_dv2.F90

Timestamp:

Jan 30, 2020, 10:15:32 AM (5 years ago)

Author:

aboissinot

Message:

Fix a bug in thermcell_dv2: now the plume height (zmin-zmax) is used for computations
instead of maximal altitude (zmax).

File:

• trunk/LMDZ.GENERIC/libf/phystd/thermcell_dv2.F90 (modified) (5 diffs)

trunk/LMDZ.GENERIC/libf/phystd/thermcell_dv2.F90

@@ -8 +8 @@
 !===============================================================================
 !
-!   Calcul du transport verticale dans la couche limite en presence
+!   Calcul du transport vertical dans la couche limite en presence
 !   de "thermiques" explicitement representes
-!   calcul du dq/dt une fois qu'on connait les ascendances
+!   Calcul du dq/dt une fois qu'on connait les ascendances
 !
 ! Vectorisation, FH : 2010/03/08
 !
 !===============================================================================
-      
-      USE print_control_mod, ONLY: prt_level, lunout
       
       IMPLICIT none
@@ -45 +43 @@
 !     --------
       
-      REAL, INTENT(out) :: ua(ngrid,nlay)
-      REAL, INTENT(out) :: va(ngrid,nlay)
-      REAL, INTENT(out) :: du(ngrid,nlay)
-      REAL, INTENT(out) :: dv(ngrid,nlay)
+      REAL, INTENT(out) :: ua(ngrid,nlay)       ! u in the plume
+      REAL, INTENT(out) :: va(ngrid,nlay)       ! v in the plume
+      REAL, INTENT(out) :: du(ngrid,nlay)       ! large scale u variation
+      REAL, INTENT(out) :: dv(ngrid,nlay)       ! large scale v variation
       
 !     Local:
@@ -91 +89 @@
          DO ig=1,ngrid
             ltherm(ig,l) = (fm(ig,l+1) + detr(ig,l)) * ptimestep > 1.e-5 * masse(ig,l)
-            IF (ltherm(ig,l).and.(zmax(ig) > 0.)) THEN
-               gamma0(ig,l) = masse(ig,l) * 0.5 / zmax(ig)                    &
+!            IF (ltherm(ig,l).and.(zmax(ig) > 0.)) THEN
+            IF (ltherm(ig,l).and.(plume_height(ig) > 0.)) THEN
+!               gamma0(ig,l) = masse(ig,l) * 0.5 / zmax(ig)                    &
+               gamma0(ig,l) = masse(ig,l) * 0.5 / plume_height(ig)            &
                &            * sqrt(0.5 * (fraca(ig,l+1) + fraca(ig,l)))
             ELSE
@@ -110 +110 @@
          DO ig=1,ngrid
             IF (ltherm(ig,l)) THEN
-               dua(ig,l)=ua(ig,l-1)-u(ig,l-1)
-               dva(ig,l)=va(ig,l-1)-v(ig,l-1)
+               dua(ig,l) = ua(ig,l-1) - u(ig,l-1)
+               dva(ig,l) = va(ig,l-1) - v(ig,l-1)
             ELSE
-               ua(ig,l)=u(ig,l)
-               va(ig,l)=v(ig,l)
-               ue(ig,l)=u(ig,l)
-               ve(ig,l)=v(ig,l)
+               ua(ig,l) = u(ig,l)
+               va(ig,l) = v(ig,l)
+               ue(ig,l) = u(ig,l)
+               ve(ig,l) = v(ig,l)
             ENDIF
          ENDDO
@@ -169 +169 @@
       DO ig=1,ngrid
          wud(ig,1) = 0.
+         wvd(ig,1) = 0.
          wud(ig,nlay+1) = 0.
-         wvd(ig,1) = 0.
          wvd(ig,nlay+1) = 0.
       ENDDO