Changeset 928

Timestamp:

Mar 18, 2008, 6:04:15 PM (17 years ago)

Author:

lmdzadmin

Message:

Bug variable d_tr_cv non-definie, cas iflag_con different de 2,3 FH
IM

File:

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

LMDZ4/trunk/libf/phylmd/phytrac.F

@@ -642 +642 @@
 c======================================================================
 c     print*,'Avant convection'
-c      do it=1,nqmax
-c         WRITE(itn,'(i2)') it
+       do it=1,nqmax
+          WRITE(itn,'(i2)') it
 c        call diagtracphy(tr_seri(:,:,it),paprs,'Avant conv'//itn)
-c      enddo
+       enddo
 
       if (convection) then
@@ -654 +654 @@
       IF ( conv_flg(it) == 0 ) CYCLE
 #endif
-      if (iflag_con.eq.2) then
+      if (iflag_con.lt.2) then
+       d_tr_cv=0.
+      else 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(1,1,it))
-      else if (iflag_con.eq.3) then
+      else
 c KE
       call cvltr(pdtphys, da, phi, mp, paprs,pplay, tr_seri(1,1,it),
@@ -664 +666 @@
       endif
 
-      DO k = 1, nlev
-      DO i = 1, klon
+       DO k = 1, nlev
+       DO i = 1, klon
          tr_seri(i,k,it) = tr_seri(i,k,it) + d_tr_cv(i,k,it)
-      ENDDO
-      ENDDO
+       ENDDO
+       ENDDO
 #ifdef INCA
       CALL minmaxqfi(tr_seri(1,1,it),0.,1.e33,'convection it = '