/[lmdze]/trunk/dyn3d/advtrac.f
ViewVC logotype

Diff of /trunk/dyn3d/advtrac.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 91 by guez, Wed Mar 26 17:18:58 2014 UTC revision 104 by guez, Thu Sep 4 10:05:52 2014 UTC
# Line 32  contains Line 32  contains
32      REAL, save:: massem(ip1jmp1, llm)      REAL, save:: massem(ip1jmp1, llm)
33      real zdp(ip1jmp1)      real zdp(ip1jmp1)
34      REAL pbarug(ip1jmp1, llm), pbarvg(ip1jm, llm), wg(ip1jmp1, llm)      REAL pbarug(ip1jmp1, llm), pbarvg(ip1jm, llm), wg(ip1jmp1, llm)
     REAL cpuadv(nqmx)  
35    
36      INTEGER:: iadvtr = 0      INTEGER:: iadvtr = 0
37      INTEGER ij, l, iq      INTEGER ij, l, iq
# Line 124  contains Line 123  contains
123               ! Schema de Frederic Hourdin               ! Schema de Frederic Hourdin
124            ELSE IF (iadv(iq)==12) THEN            ELSE IF (iadv(iq)==12) THEN
125               ! Pas de temps adaptatif               ! Pas de temps adaptatif
126               CALL adaptdt(iadv(iq), dtbon, n, pbarug, massem)               CALL adaptdt(dtbon, n, pbarug, massem)
127               IF (n>1) THEN               IF (n>1) THEN
128                  WRITE (*, *) 'WARNING horizontal dt=', dtbon, 'dtvr=', dtvr, &                  WRITE (*, *) 'WARNING horizontal dt=', dtbon, 'dtvr=', dtvr, &
129                       'n=', n                       'n=', n
# Line 134  contains Line 133  contains
133               END DO               END DO
134            ELSE IF (iadv(iq)==13) THEN            ELSE IF (iadv(iq)==13) THEN
135               ! Pas de temps adaptatif               ! Pas de temps adaptatif
136               CALL adaptdt(iadv(iq), dtbon, n, pbarug, massem)               CALL adaptdt(dtbon, n, pbarug, massem)
137               IF (n>1) THEN               IF (n>1) THEN
138                  WRITE (*, *) 'WARNING horizontal dt=', dtbon, 'dtvr=', dtvr, &                  WRITE (*, *) 'WARNING horizontal dt=', dtbon, 'dtvr=', dtvr, &
139                       'n=', n                       'n=', n
# Line 148  contains Line 147  contains
147               ! Schema de Prather               ! Schema de Prather
148            ELSE IF (iadv(iq)==30) THEN            ELSE IF (iadv(iq)==30) THEN
149               ! Pas de temps adaptatif               ! Pas de temps adaptatif
150               CALL adaptdt(iadv(iq), dtbon, n, pbarug, massem)               CALL adaptdt(dtbon, n, pbarug, massem)
151               IF (n>1) THEN               IF (n>1) THEN
152                  WRITE (*, *) 'WARNING horizontal dt=', dtbon, 'dtvr=', dtvr, &                  WRITE (*, *) 'WARNING horizontal dt=', dtbon, 'dtvr=', dtvr, &
153                       'n=', n                       'n=', n
# Line 158  contains Line 157  contains
157            ELSE IF (iadv(iq)==11 .OR. (iadv(iq)>=16 .AND. iadv(iq)<=18)) THEN            ELSE IF (iadv(iq)==11 .OR. (iadv(iq)>=16 .AND. iadv(iq)<=18)) THEN
158               ! Test sur le flux horizontal               ! Test sur le flux horizontal
159               ! Pas de temps adaptatif               ! Pas de temps adaptatif
160               CALL adaptdt(iadv(iq), dtbon, n, pbarug, massem)               CALL adaptdt(dtbon, n, pbarug, massem)
161               IF (n>1) THEN               IF (n>1) THEN
162                  WRITE (*, *) 'WARNING horizontal dt=', dtbon, 'dtvr=', dtvr, &                  WRITE (*, *) 'WARNING horizontal dt=', dtbon, 'dtvr=', dtvr, &
163                       'n=', n                       'n=', n
Legend:
Removed from v.91  
changed lines
  Added in v.104
  ViewVC Help
Powered by ViewVC 1.1.21