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

Timestamp:

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

Author:

lmdzadmin

Message:

Rajout convection Kerry Emanuel pour traceurs- MAF+JYG

File:

• LMDZ4/trunk/libf/phylmd/phytrac.F (modified) (9 diffs)

LMDZ4/trunk/libf/phylmd/phytrac.F

@@ -4 +4 @@
 c
 c
-      SUBROUTINE phytrac (rnpb,nstep,
+      SUBROUTINE phytrac (iflag_con,rnpb,nstep,
      I                    julien,gmtime,
      I                    debutphy,lafin,
@@ -146 +146 @@
       REAL flxmass_w(klon,klev)
 #endif
+      integer iflag_con
 
 cAA Rem : nbtr : nombre de vrais traceurs est defini dans dimphy.h
@@ -166 +167 @@
       REAL pen_d(nlon,nlev) ! flux entraine dans le panache descendant
       REAL pde_d(nlon,nlev) ! flux detraine dans le panache descendant
+c KE
+      real da(nlon,nlev),phi(nlon,nlev,nlev),mp(nlon,nlev)
+      REAL upwd(nlon,nlev)      ! saturated updraft mass flux
+      REAL dnwd(nlon,nlev)      ! saturated downdraft mass flux
+
 c
 c   Couche limite:
@@ -281 +287 @@
 c
       REAL d_tr(klon,klev), d_trs(klon) ! tendances de traceurs 
-      REAL d_tr_cl(klon,klev) ! tendance de traceurs  couche limite
-      REAL d_tr_cli(klon,klev,nbtr) ! tendance de traceurs  CL pour chq traceur
-      REAL d_tr_cv(klon,klev) ! tendance de traceurs  convection
-      REAL d_tr_cvi(klon,klev,nbtr) ! tendance de traceurs  conv pour chq traceur
+      REAL d_tr_cl(klon,klev,nbtr) ! tendance de traceurs  couche limite
+      REAL d_tr_cv(klon,klev,nbtr) ! tendance de traceurs  conv pour chq traceur
       REAL d_tr_th(klon,klev,nbtr) ! la tendance des thermiques
       REAL d_tr_dec(klon,klev,nbtr) ! la tendance de la decroissance 
@@ -528 +532 @@
 
 c Abder
-      if(nqmax.gt.2) aerosol(3)=.true.
+ctestmaf      if(nqmax.gt.2) aerosol(3)=.true.
 
        do i=1,nlon
@@ -548 +552 @@
 c======================================================================
 c     print*,'Avant convection'
-      do it=1,nqmax
-         WRITE(itn,'(i2)') it
+c      do it=1,nqmax
+c         WRITE(itn,'(i2)') it
 c        call diagtracphy(tr_seri(:,:,it),paprs,'Avant conv'//itn)
-      enddo
+c      enddo
 
       if (convection) then
@@ -560 +564 @@
       IF ( conv_flg(it) == 0 ) CYCLE
 #endif
+      if (iflag_con.eq.2) then
+c tiedke
       CALL nflxtr(pdtphys, pmfu, pmfd, pen_u, pde_u, pen_d, pde_d,
-     .            pplay, paprs, tr_seri(1,1,it), d_tr_cv)
+     .            pplay, paprs, tr_seri(1,1,it), d_tr_cv(1,1,it))
+      else if (iflag_con.eq.3) then
+c KE
+      call cvltr(pdtphys, da, phi, mp, paprs,pplay, tr_seri(1,1,it),
+     .           upwd,dnwd,d_tr_cv(1,1,it))
+      endif
+
       DO k = 1, nlev
       DO i = 1, klon
-         tr_seri(i,k,it) = tr_seri(i,k,it) + d_tr_cv(i,k)
-         d_tr_cvi(i,k,it)=d_tr_cv(i,k)
-c        print*,'en k i d_tr_cv=',k,i,d_tr_cv(i,k)
+         tr_seri(i,k,it) = tr_seri(i,k,it) + d_tr_cv(i,k,it)
       ENDDO
       ENDDO
-c      WRITE(itn,'(i1)') it
 #ifdef INCA
       CALL minmaxqfi(tr_seri(1,1,it),0.,1.e33,'convection it = '
@@ -676 +685 @@
      e                    tautr(it),vdeptr(it),
      e                    xlat,
-     s                    d_tr_cl,d_trs)
+     s                    d_tr_cl(1,1,it),d_trs)
           DO k = 1, nlev
             DO i = 1, klon
-              tr_seri(i,k,it) = tr_seri(i,k,it) + d_tr_cl(i,k)
-              d_tr_cli(i,k,it)=d_tr_cl(i,k)
+              tr_seri(i,k,it) = tr_seri(i,k,it) + d_tr_cl(i,k,it)
             ENDDO
           ENDDO
@@ -709 +717 @@
      s               tr_seri(1,1,it), source,
      e               paprs, pplay, delp,
-     s               d_tr )
+     s               d_tr_cl(1,1,it))
             DO k = 1, nlev
                DO i = 1, klon
-                  tr_seri(i,k,it) = tr_seri(i,k,it) + d_tr(i,k)
-                  d_tr_cli(i,k,it)=d_tr_cl(i,k)
+                  tr_seri(i,k,it) = tr_seri(i,k,it) + d_tr_cl(i,k,it)
                ENDDO
             ENDDO
-Cmaf provisoire suppression des prints
 Cmaf          WRITE(itn,'(i1)') it
 cmaf          CALL minmaxqfi(tr_seri(1,1,it),0.,1.e33,'cltracn it='//itn)