Changeset 619 for LMDZ4/trunk/libf/phylmd/concvl.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/concvl.F (modified) (6 diffs)

LMDZ4/trunk/libf/phylmd/concvl.F

@@ -7 +7 @@
      .             upwd,dnwd,dnwdbis,Ma,cape,tvp,iflag,
      .             pbase,bbase,dtvpdt1,dtvpdq1,dplcldt,dplcldr,
-     .             qcondc,wd)
+     .             qcondc,wd,
+     .             pmflxr,pmflxs,
+     .             da,phi,mp)
  
 c
@@ -56 +58 @@
        INTEGER ntra
        REAL work1(klon,klev),work2(klon,klev)
+       REAL pmflxr(klon,klev+1),pmflxs(klon,klev+1)
 c
        REAL d_t(klon,klev),d_q(klon,klev),d_u(klon,klev),d_v(klon,klev)
@@ -65 +68 @@
        REAL upwd(klon,klev),dnwd(klon,klev),dnwdbis(klon,klev)
        REAL Ma(klon,klev),cape(klon),tvp(klon,klev)
+       real da(klon,klev),phi(klon,klev,klev),mp(klon,klev)
        INTEGER iflag(klon)
        REAL rflag(klon)
@@ -101 +105 @@
          DO i=1,klon
          em_ph(i,k) = paprs(i,k) / 100.0
+         pmflxs(i,k)=0.
       ENDDO
       ENDDO
@@ -145 +150 @@
      $              em_p,em_ph,iflag,
      $              d_t,d_q,d_u,d_v,d_tra,rain,
-     $              cbmf,work1,work2,
-     $              dtime,Ma,upwd,dnwd,dnwdbis,qcondc,wd,cape)
+     $              pmflxr,cbmf,work1,work2,
+     $              kbas,ktop,
+     $              dtime,Ma,upwd,dnwd,dnwdbis,qcondc,wd,cape,
+     $              da,phi,mp)
 
 C------------------------------------------------------------------
@@ -163 +170 @@
         ENDDO
       ENDDO
- 
+       DO itra = 1,ntra
+        DO k = 1, klev
+         DO i = 1, klon
+            d_tra(i,k,itra) =dtime*d_tra(i,k,itra) 
+         ENDDO
+        ENDDO
+       ENDDO
 c les traceurs ne sont pas mis dans cette version de convect4:
       if (iflag_con.eq.4) then