Changeset 2113 for trunk/LMDZ.GENERIC/libf

Timestamp:

Mar 6, 2019, 7:10:00 PM (6 years ago)

Author:

aboissinot

Message:

Add new formulae to compute vertical speed in thermcell_plume (only consistent with mix0=0)

Add pdt and pdq to pt and pq for computing zdttherm and zdqtherm in physiq_mod when water=false.
That in order to be consistent with the case water=true

Location:

trunk/LMDZ.GENERIC/libf/phystd

Files:

• physiq_mod.F90 (modified) (1 diff)
• thermcell_plume.F90 (modified) (5 diffs)

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

@@ -1193 +1193 @@
             CALL evap(ngrid,nlayer,nq,ptimestep,pt,pq,pdq,pdt,dqevap,dtevap,zqtherm,zttherm)
          ELSE
-            zttherm(:,:) = pt(:,:)
-            zqtherm(:,:,:) = pq(:,:,:)
+            zttherm(:,:)   = pt(:,:)   + pdt(:,:)
+            zqtherm(:,:,:) = pq(:,:,:) + pdq(:,:,:)
          ENDIF
          

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

@@ -31 +31 @@
 !==============================================================================
       
-!      inputs:
-!      -------
+!     Inputs:
+!     -------
       
       INTEGER itap
@@ -52 +52 @@
       REAL zpspsk(ngrid,nlay)                   ! Exner function
       
-!      outputs:
-!      --------
+!     Outputs:
+!     --------
       
       INTEGER lmin(ngrid)                       ! plume base level (first unstable level)
@@ -81 +81 @@
       REAL zqsatth(ngrid,nlay)                  ! saturation vapor pressure (after mixing)
       
-!      local:
-!      ------
+!     Local:
+!     ------
       
       INTEGER ig, l, k
@@ -305 +305 @@
 ! AB : own derivation for w_est (Rio 2010 formula with e=d=0)
 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
-               zw2fact = 2. * fact_epsilon * zdz
-               zdw2 = 2. * afact * zbuoy(ig,l) * zdz
+!               zw2fact = 2. * fact_epsilon * zdz
+!               zdw2 = 2. * afact * zbuoy(ig,l) * zdz
+               zw2fact = 2. * fact_epsilon * zdz / (1. + betalpha)
+               zdw2 = 2. * afact * zbuoy(ig,l) * zdz / (1. + betalpha)
                w_est(ig,l+1) = Max(0., exp(-zw2fact) * w_est(ig,l) + zdw2)
                
@@ -459 +461 @@
 !               zw2fact = fact_epsilon * 2. * zdz / (1. + betalpha)
 !               zdw2 = afact * zbuoy(ig,l) / fact_epsilon
-!               zdw2bis = afact * zbuoy(ig,l-1) / fact_epsilon
 !               zw2(ig,l+1) = Max(0.0001,exp(-zw2fact)*(zw2(ig,l)-zdw2)+zdw2)
 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 ! AB : own derivation for zw2 (Rio 2010 formula)
 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
-               zw2fact = 2. * (fact_epsilon * zdz  + entr_star(ig,l) / f_star(ig,l))
-               zdw2 = 2. * afact * zbuoy(ig,l) * zdz
+!               zw2fact = 2. * (fact_epsilon * zdz  + entr_star(ig,l) / f_star(ig,l))
+!               zdw2 = 2. * afact * zbuoy(ig,l) * zdz
+               zw2fact = 2. * fact_epsilon * zdz / (1. + betalpha)
+               zdw2 = 2. * afact * zbuoy(ig,l) * zdz / (1. + betalpha)
                zw2(ig,l+1) = Max(0., exp(-zw2fact) * zw2(ig,l) + zdw2)