Changeset 619 for LMDZ4/trunk/libf/phylmd/cv_driver.F

Timestamp:

Apr 15, 2005, 2:36:17 PM (19 years ago)

Author:

lmdzadmin

Message:

Rajout convection Kerry Emanuel pour traceurs- MAF+JYG

File:

• LMDZ4/trunk/libf/phylmd/cv_driver.F (modified) (14 diffs)

LMDZ4/trunk/libf/phylmd/cv_driver.F

@@ -5 +5 @@
      &                   t1,q1,qs1,u1,v1,tra1,
      &                   p1,ph1,iflag1,ft1,fq1,fu1,fv1,ftra1,
-     &                   precip1,
+     &                   precip1,VPrecip1,
      &                   cbmf1,sig1,w01,
-     &                   delt,Ma1,upwd1,dnwd1,dnwd01,qcondc1,wd1,cape1)
+     &                   icb1,inb1,
+     &                   delt,Ma1,upwd1,dnwd1,dnwd01,qcondc1,wd1,cape1,
+     &                   da1,phi1,mp1)
 C
       implicit none
@@ -37 +39 @@
 C      ftra1         Real           Output       tracor tend
 C      precip1       Real           Output       precipitation
+C      VPrecip1      Real           Output       vertical profile of precipitations
 C      cbmf1         Real           Output       cloud base mass flux
 C      sig1          Real           In/Out       section adiabatic updraft
@@ -84 +87 @@
       real precip1(len)
       real cbmf1(len)
+      real VPrecip1(len,nd+1)
       real Ma1(len,nd)
       real upwd1(len,nd)
@@ -93 +97 @@
       real cape1(len)     
 
+      real da1(len,nd),phi1(len,nd,nd),mp1(len,nd)
+      real da(len,nd),phi(len,nd,nd),mp(len,nd)
       real tra1(len,nd,ntra)
       real ftra1(len,nd,ntra)
@@ -185 +191 @@
 !  precip: Scalar convective precipitation rate (mm/day).
 !
+!  VPrecip: Vertical profile of convective precipitation (kg/m2/s).
+!
 !  wd:   A convective downdraft velocity scale. For use in surface
 !        flux parameterizations. See convect.ps file for details.
@@ -213 +221 @@
       integer nk1(klon)
       integer icb1(klon)
+      integer inb1(klon)
       integer icbs1(klon)
 
@@ -267 +276 @@
       real ments(nloc,klev,klev), qents(nloc,klev,klev)
       real sij(nloc,klev,klev), elij(nloc,klev,klev)
-      real mp(nloc,klev), qp(nloc,klev), up(nloc,klev), vp(nloc,klev)
+      real qp(nloc,klev), up(nloc,klev), vp(nloc,klev)
       real wt(nloc,klev), water(nloc,klev), evap(nloc,klev)
       real b(nloc,klev), ft(nloc,klev), fq(nloc,klev)
@@ -275 +284 @@
       real tps(nloc,klev), qprime(nloc), tprime(nloc)
       real precip(nloc)
+      real VPrecip(nloc,klev+1)
       real tra(nloc,klev,ntra), trap(nloc,klev,ntra)
       real ftra(nloc,klev,ntra), traent(nloc,klev,klev,ntra)
@@ -324 +334 @@
          clw(i,k)=0.0    
          gz1(i,k) = 0.
-
+         VPrecip1(i,k) = 0.
          Ma1(i,k)=0.0
          upwd1(i,k)=0.0
@@ -346 +356 @@
         wd1(i)=0.0
         cape1(i)=0.0
+        VPrecip1(i,nd+1)=0.0
  60   continue
 
@@ -552 +563 @@
      :                     ,ment,qent,uent,vent,nent,elij,traent,sig
      :                     ,tv,tvp
-     o                     ,iflag,precip,ft,fq,fu,fv,ftra
+     o                     ,iflag,precip,VPrecip,ft,fq,fu,fv,ftra
      o                     ,upwd,dnwd,dnwd0,ma,mike,tls,tps,qcondc,wd)
       endif
@@ -568 +579 @@
 
 !^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
+! --- passive tracers
+!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
+
+      if (iflag_con.eq.3) then
+       CALL cv3_tracer(nloc,len,ncum,nd,nd,
+     :                  ment,sij,da,phi)
+      endif
+
+!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
 ! --- UNCOMPRESS THE FIELDS
 !^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
-
-
+c set iflag1 =42 for non convective points 
+      do  i=1,len
+        iflag1(i)=42
+      end do
+c
       if (iflag_con.eq.3) then
        CALL cv3_uncompress(nloc,len,ncum,nd,ntra,idcum
      :          ,iflag
-     :          ,precip,sig,w0
+     :          ,precip,VPrecip,sig,w0
      :          ,ft,fq,fu,fv,ftra
+     :          ,inb 
      :          ,Ma,upwd,dnwd,dnwd0,qcondc,wd,cape
+     :          ,da,phi,mp
      o          ,iflag1
-     o          ,precip1,sig1,w01
+     o          ,precip1,VPrecip1,sig1,w01
      o          ,ft1,fq1,fu1,fv1,ftra1
-     o          ,Ma1,upwd1,dnwd1,dnwd01,qcondc1,wd1,cape1 )
+     o          ,inb1
+     o          ,Ma1,upwd1,dnwd1,dnwd01,qcondc1,wd1,cape1 
+     o          ,da1,phi1,mp1)
       endif
 
@@ -611 +638 @@
 c -- si .TRUE., on rend la gravite plus explicite et eventuellement
 c differente de 10.0 dans convect3: 
-      cvflag_grav = .FALSE.
+      cvflag_grav = .TRUE.
 
       return
@@ -654 +681 @@
 c ori      t0  = RTT
        t0  = 273.15 ! convect3 (RTT=273.16)
-       grav= 10.    ! implicitely or explicitely used in convect3
+c maf       grav= 10.    ! implicitely or explicitely used in convect3
+       grav= g    ! implicitely or explicitely used in convect3
       endif