Changeset 2093

Timestamp:

Feb 7, 2019, 12:01:36 PM (6 years ago)

Author:

aboissinot

Message:

d_temp is now taken into acocunt in thermcell_plume to trigger the plume and compute first unstable layer speed.
Modified temperature is computed in thermcell_plume and thermcell_alim get it via its arguments.

Location:

trunk/LMDZ.GENERIC

Files:

• README (modified) (1 diff)
• libf/phystd/thermcell_alim.F90 (modified) (4 diffs)
• libf/phystd/thermcell_plume.F90 (modified) (6 diffs)

trunk/LMDZ.GENERIC/README

@@ -1448 +1448 @@
 - uncomment two "corrections" in thermcell_flux. They are used only if iflag_thermals_optflux is set to 0
 - remove useless parameter fact_shell in thermcell_mod
+- now take d_temp into acocunt in thermcell_plume to trigger the plume and compute first unstable layer speed.

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

@@ -2 +2 @@
 !
 !
-      SUBROUTINE thermcell_alim(ngrid,klev,ztv,zlev,alim_star,lalim,lmin)
+      SUBROUTINE thermcell_alim(ngrid,klev,ztv2,zlev,alim_star,lalim,lmin)
       
       
@@ -27 +27 @@
       INTEGER, INTENT(IN) :: lmin(ngrid)        ! plume initial level
       
-      REAL, INTENT(IN) :: ztv(ngrid,klev)       ! Large scale virtual potential temperature
-! Virtual potential temperature offset in linf layer
+      REAL, INTENT(IN) :: ztv2(ngrid,klev)      ! Virtual potential temperature
       REAL, INTENT(IN) :: zlev(ngrid,klev+1)    ! levels altitude
       
@@ -44 +43 @@
       
       REAL :: alim_star_tot(ngrid)              ! integrated alimentation
-      REAL :: ztv2(ngrid,klev)                  ! ztv + d_temp
       
 !==============================================================================
@@ -53 +51 @@
       
       alim_star_tot(:) = 0.
-      
-      ztv2(:,:) = ztv(:,:)
-      ztv2(:,linf) = ztv2(:,linf) + d_temp
       
 !==============================================================================

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

@@ -96 +96 @@
       REAL zqsat(ngrid)                         ! saturation vapor pressure (before mixing)
       REAL zdz                                  ! layers height
+      REAL ztv2(ngrid,klev)                     ! ztv + d_temp * Dirac(l=lmin)
       
       REAL zalpha                               !
@@ -159 +160 @@
       lmin(:)          = linf
       
+      ztv2(:,:)        = ztv(:,:)
+      ztv2(:,linf)     = ztv2(:,linf) + d_temp
+      
 !==============================================================================
 ! 0. Calcul de l'alimentation
@@ -175 +179 @@
          l = linf
          DO WHILE ((.not.active(ig)) .and. pplev(ig,l+1).gt.pres_limit .and. l.lt.klev)
-            IF (ztv(ig,l).gt.ztv(ig,l+1)) THEN
+            IF (ztv2(ig,l).gt.ztv2(ig,l+1)) THEN
                active(ig) = .true.
                lmin(ig) = l
@@ -186 +190 @@
 ! AB : On pourrait n'appeler thermcell_alim que si la plume est active
 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
-      CALL thermcell_alim(ngrid,klev,ztv,zlev,alim_star,lalim,lmin)
+      CALL thermcell_alim(ngrid,klev,ztv2,zlev,alim_star,lalim,lmin)
       
 !==============================================================================
@@ -208 +212 @@
             
             zw2(ig,l+1) = 2. * RG * (zlev(ig,l+1) - zlev(ig,l))               &
-            &                * (ztv(ig,l) - ztv(ig,l+1)) / ztv(ig,l+1)
+            &                * (ztv2(ig,l) - ztv2(ig,l+1)) / ztv2(ig,l+1)
             
             w_est(ig,l+1) = zw2(ig,l+1)
-            
          ENDIF
       ENDDO
@@ -488 +491 @@
       ENDDO
       
-#undef wrgrads_thermcell
-#ifdef wrgrads_thermcell
-      call wrgradsfi(1,klev,entr_star(igout,1:klev),'esta    ','esta    ')
-      call wrgradsfi(1,klev,detr_star(igout,1:klev),'dsta    ','dsta    ')
-      call wrgradsfi(1,klev,zbuoy(igout,1:klev)    ,'buoy    ','buoy    ')
-      call wrgradsfi(1,klev,zdqt(igout,1:klev)     ,'dqt     ','dqt     ')
-      call wrgradsfi(1,klev,w_est(igout,1:klev)    ,'w_est   ','w_est   ')
-      call wrgradsfi(1,klev,w_est(igout,2:klev+1)  ,'w_es2   ','w_es2   ')
-      call wrgradsfi(1,klev,zw2(igout,1:klev)      ,'zw2A    ','zw2A    ')
-#endif
-      
       
 RETURN