Changeset 2205 for LMDZ5/trunk/libf/phylmd/cv3_routines.F90

Timestamp:

Feb 13, 2015, 8:22:32 AM (10 years ago)

Author:

fhourdin

Message:

Introduction of a bi-gausian descrption of the subgrid scale
distribution of water associated with deep convection.

Introduction d'un schéma bi-gausien pour la distribution sous maille
de l'eau associée à la convection profonde.

Contrôlé par
iflag_cld_cv=2 (1 pour le schéma de Bony et Emanuel 2001)
autres paramètres de contrôle : tau_cld_cv et coefw_cld_cv

Arnaud Jam

File:

• LMDZ5/trunk/libf/phylmd/cv3_routines.F90 (modified) (11 diffs)

LMDZ5/trunk/libf/phylmd/cv3_routines.F90

@@ -1 +1 @@
 
 ! $Id$
+
 
 
@@ -2763 +2764 @@
                      cbmf, upwd, dnwd, dnwd0, ma, mip, &
                      tls, tps, qcondc, wd, &
-                     ftd, fqd)
+                     ftd, fqd, qnk, qtc, sigt, tau_cld_cv, coefw_cld_cv)
 
   IMPLICIT NONE
@@ -2800 +2801 @@
 !input/output:
       INTEGER iflag(nloc)
+      REAL,INTENT(in) :: tau_cld_cv, coefw_cld_cv
 !
 !outputs:
@@ -2811 +2813 @@
       REAL tls(nloc, nd), tps(nloc, nd)
       REAL qcondc(nloc, nd) ! cld
+      REAL qtc(nloc,nd), sigt(nloc,nd) ! cld
       REAL wd(nloc) ! gust
       REAL cbmf(nloc)
@@ -2830 +2833 @@
       REAL qcond(nloc, nd), nqcond(nloc, nd), wa(nloc, nd) ! cld
       REAL siga(nloc, nd), sax(nloc, nd), mac(nloc, nd) ! cld
-
+      REAL sument(nloc), sigment(nloc,nd), qtment(nloc,nd) ! cld
+      REAL qnk(nloc)
       REAL sumdq !jyg
 ! 
@@ -2861 +2865 @@
       qcondc(il, i) = 0.0 ! cld
       qcond(il, i) = 0.0 ! cld
+      qtc(il, i) = 0.0 ! cld
+      qtment(il, i) = 0.0 ! cld
+      sigment(il, i) = 0.0 ! cld
+      sigt(il, i) = 0.0 ! cld
       nqcond(il, i) = 0.0 ! cld
     END DO
@@ -3234 +3242 @@
 ! (saturated updrafts resulting from mixing)                                   ! cld
           qcond(il, i) = qcond(il, i) + (elij(il,k,i)-awat(il))                ! cld
-          nqcond(il, i) = nqcond(il, i) + 1. ! cld
+          qtment(il, i) = qtment(il, i) + qent(il,k,i)                         ! cld
+          nqcond(il, i) = nqcond(il, i) + 1.                                   ! cld
         END IF ! i
       END DO
@@ -3310 +3319 @@
 ! (saturated downdrafts resulting from mixing)                                 ! cld
           qcond(il, i) = qcond(il, i) + elij(il, k, i)                         ! cld
-          nqcond(il, i) = nqcond(il, i) + 1. ! cld
+          qtment(il, i) = qent(il,k,i) + qtment(il,i)                          ! cld
+          nqcond(il, i) = nqcond(il, i) + 1.                                   ! cld
         END IF ! cld
       END DO ! cld
@@ -3319 +3329 @@
       IF (i<=inb(il) .AND. nent(il,i)==0 .AND. iflag(il)<=1) THEN              ! cld
         qcond(il, i) = qcond(il, i) + (1.-ep(il,i))*clw(il, i)                 ! cld
+        qtment(il, i) = qent(il,k,i) + qtment(il,i)                          ! cld
         nqcond(il, i) = nqcond(il, i) + 1.                                     ! cld
       END IF                                                                   ! cld
@@ -3326 +3337 @@
       IF (i<=inb(il) .AND. nqcond(il,i)/=0 .AND. iflag(il)<=1) THEN            ! cld
         qcond(il, i) = qcond(il, i)/nqcond(il, i)                              ! cld
+        qtment(il, i) = qtment(il,i)/nqcond(il, i)                             ! cld
       END IF                                                                   ! cld
     END DO
@@ -3788 +3800 @@
       END IF                                                         ! cld
     END DO                                                           ! cld
-  END DO                                                             ! cld
-
-  DO i = 1, nl                                                       ! cld
+  END DO  
+                                                           ! cld
+  DO i = 1, nl  
+
+! 14/01/15 AJ je remets les parties manquantes cf JYG
+! Initialize sument to 0
+
+    DO il = 1,ncum
+     sument(il) = 0.
+    ENDDO
+
+! Sum mixed mass fluxes in sument
+
+    DO k = 1,nl
+      DO il = 1,ncum
+        IF  (k<=inb(il) .AND. i<=inb(il) .AND. iflag(il)<=1) THEN   ! cld
+          sument(il) =sument(il) + abs(ment(il,k,i))
+        ENDIF
+      ENDDO     ! il
+    ENDDO       ! k
+
+! 14/01/15 AJ delta n'a rien à faire là...                                                 
     DO il = 1, ncum                                                  ! cld
       IF (wa(il,i)>0.0 .AND. iflag(il)<=1) &                         ! cld
-        siga(il, i) = mac(il, i)/wa(il, i) &                         ! cld
-        *rrd*tvp(il, i)/p(il, i)/100./delta                          ! cld
+        siga(il, i) = mac(il, i)/(coefw_cld_cv*wa(il, i)) &          ! cld
+        *rrd*tvp(il, i)/p(il, i)/100.                                ! cld
+
       siga(il, i) = min(siga(il,i), 1.0)                             ! cld
-! IM cf. FH                                                          
+
+! IM cf. FH 
+! 14/01/15 AJ ne correspond pas à ce qui a été codé par JYG et SB           
+                                                         
       IF (iflag_clw==0) THEN                                         ! cld
         qcondc(il, i) = siga(il, i)*clw(il, i)*(1.-ep(il,i)) &       ! cld
           +(1.-siga(il,i))*qcond(il, i)                              ! cld
+
+
+        sigment(il,i)=sument(il)*tau_cld_cv/(ph(il,i)-ph(il,i+1))    ! cld
+        sigment(il, i) = min(1.e-4+sigment(il,i), 1.0 - siga(il,i))  ! cld
+        qtc(il, i) = (siga(il,i)*qnk(il)+sigment(il,i)*qtment(il,i)) & ! cld
+                     /(siga(il,i)+sigment(il,i))                     ! cld
+        sigt(il,i) = sigment(il, i) + siga(il, i)
+
+!        qtc(il, i) = siga(il,i)*qnk(il)+(1.-siga(il,i))*qtment(il,i) ! cld
+      print*,'BIGAUSSIAN CONV',siga(il,i),sigment(il,i),qtc(il,i)  
+                
       ELSE IF (iflag_clw==1) THEN                                    ! cld
         qcondc(il, i) = qcond(il, i)                                 ! cld
+        qtc(il,i) = qtment(il,i)                                     ! cld
       END IF                                                         ! cld