Changeset 2143 for trunk/LMDZ.GENERIC/libf/phystd/thermcell_env.F90

Timestamp:

Jun 20, 2019, 4:11:54 PM (5 years ago)

Author:

aboissinot

Message:

Update the thermal plume model:

• check formulae consistency between thermcell_plume and thercell_closure
• compute correctly thermal plume height
• fix alimentation computation in the first unstable layer

File:

• trunk/LMDZ.GENERIC/libf/phystd/thermcell_env.F90 (modified) (4 diffs)

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

@@ -51 +51 @@
 !     --------
       
-      REAL zqt(ngrid,nlay)                      ! q    environment
-      REAL zql(ngrid,nlay)                      ! q    environment
+      REAL zqt(ngrid,nlay)                      ! qt   environment
+      REAL zql(ngrid,nlay)                      ! ql   environment
       REAL zt(ngrid,nlay)                       ! T    environment
       REAL ztv(ngrid,nlay)                      ! TRPV environment
@@ -76 +76 @@
       zhl(:,:) = pt(:,:) / zpopsk(:,:)
       
+      zqt(:,:) = pq(:,:,igcm_h2o_vap)
+      zql(:,:) = 0.
+      
 !===============================================================================
 ! Condensation and latent heat release
@@ -81 +84 @@
       
       IF (water) THEN
-         
-         zqt(:,:) = pq(:,:,igcm_h2o_vap)
          
          DO l=1,nlay
@@ -101 +102 @@
       ELSE
          
-         zqt(:,:) = 0.
-         zql(:,:) = 0.
          zt(:,:) = pt(:,:)
-         ztv(:,:) = zhl(:,:)
+         ztv(:,:) = pt(:,:) / zpopsk(:,:)
          
       ENDIF