Changeset 33 for LMDZ.3.3


Ignore:
Timestamp:
Feb 7, 2000, 4:43:34 PM (25 years ago)
Author:
lmdz
Message:

Toutes references aux traceurs ont ete eliminees par L. Li
LF

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ.3.3/trunk/libf/phylmd/ajsec.F

    r2 r33  
    1       SUBROUTINE ajsec(paprs, pplay, t,q, d_t,d_q, itr, tr, d_tr)
     1      SUBROUTINE ajsec(paprs, pplay, t,q, d_t,d_q)
    22      IMPLICIT none
    33c======================================================================
     
    1616      REAL t(klon,klev), q(klon,klev)
    1717      REAL d_t(klon,klev), d_q(klon,klev)
    18 c
    19       INTEGER it, itr
    20       REAL tr(klon,klev,nbtr), d_tr(klon,klev,nbtr)
    21       REAL ztr(klon,klev,nbtr), trm(nbtr)
    2218c
    2319      INTEGER limbas, limhau ! les couches a ajuster
     
    4541      ENDDO
    4642      ENDDO
    47       IF (itr.GE.1) THEN
    48       DO it = 1, itr
    49       DO k = 1, klev
    50       DO i = 1, klon
    51          d_tr(i,k,it) = 0.0
    52       ENDDO
    53       ENDDO
    54       ENDDO
    55       ENDIF
    5643c------------------------------------- detection des profils a modifier
    5744      DO k = limbas, limhau
     
    6249      ENDDO
    6350      ENDDO
    64 c
    65       IF (itr.GE.1) THEN
    66       DO it = 1, itr
    67       DO k = limbas, limhau
    68       DO i = 1, klon
    69          ztr(i,k,it) = tr(i,k,it)
    70       ENDDO
    71       ENDDO
    72       ENDDO
    73       ENDIF
    7451c
    7552      DO k = limbas, limhau
     
    10279              hm = zh(i,k2)
    10380              qm = zq(i,k2)
    104               IF (itr.GE.1) THEN
    105               DO it = 1, itr
    106                 trm(it) = ztr(i,k2,it)
    107               ENDDO
    108               ENDIF
    10981 8020         CONTINUE
    11082                sm = sm +zpkdp(i,k)
    11183                hm = hm +zpkdp(i,k) * (zh(i,k)-hm) / sm
    11284                qm = qm +zpkdp(i,k) * (zq(i,k)-qm) / sm
    113                 IF (itr.GE.1) THEN
    114                 DO it = 1, itr
    115                 trm(it) = trm(it) +
    116      .                    zpkdp(i,k) * (ztr(i,k,it)-trm(it)) / sm
    117                 ENDDO
    118                 ENDIF
    11985                down = .FALSE.
    12086                IF (k1 .ne. limbas) THEN
     
    137103                zq(i,k) = qm
    138104              ENDDO
    139               IF (itr.GE.1) THEN
    140               DO it = 1, itr
    141               DO k = k1, k2
    142                 ztr(i,k,it) = trm(it)
    143               ENDDO
    144               ENDDO
    145               ENDIF
    146105              k2 = k2 + 1
    147106            ENDIF
     
    173132         d_q(i,k) = 0.0
    174133      ENDDO
    175       ENDDO
    176       ENDIF
    177 c
    178       IF (itr.GE.1) THEN
    179       DO it = 1, itr
    180 c
    181       DO k = limbas, limhau
    182       DO i = 1, klon
    183          d_tr(i,k,it) = ztr(i,k,it) - tr(i,k,it)
    184       ENDDO
    185       ENDDO
    186 c
    187       IF (limbas.GT.1) THEN
    188       DO k = 1, limbas-1
    189       DO i = 1, klon
    190          d_tr(i,k,it) = 0.0
    191       ENDDO
    192       ENDDO
    193       ENDIF
    194 c
    195       IF (limhau.LT.klev) THEN
    196       DO k = limhau+1, klev
    197       DO i = 1, klon
    198          d_tr(i,k,it) = 0.0
    199       ENDDO
    200       ENDDO
    201       ENDIF
    202 c
    203134      ENDDO
    204135      ENDIF
Note: See TracChangeset for help on using the changeset viewer.