78 |
lat = jjp1 |
lat = jjp1 |
79 |
niv = llm |
niv = llm |
80 |
|
|
|
C *** Test de passage d'arguments ****** |
|
|
|
|
|
c DO 399 l = 1, llm |
|
|
c DO 399 j = 1, jjp1 |
|
|
c DO 399 i = 1, iip1 |
|
|
c IF (S0(i,j,l,ntra) .lt. 0. ) THEN |
|
|
c PRINT*,'S0(',i,j,l,')=',S0(i,j,l,ntra) |
|
|
c print*, 'sx(',i,j,l,')=',sx(i,j,l,ntra) |
|
|
c print*, 'sy(',i,j,l,')=',sy(i,j,l,ntra) |
|
|
c print*, 'sz(',i,j,l,')=',sz(i,j,l,ntra) |
|
|
c PRINT*, 'AIE !! debut ADVZ - pbl arg. passage dans ADVZ' |
|
|
c STOP |
|
|
c ENDIF |
|
|
399 CONTINUE |
|
|
|
|
|
C----------------------------------------------------------------- |
|
81 |
C *** Test : diag de la qqtite totale de traceur |
C *** Test : diag de la qqtite totale de traceur |
82 |
C dans l'atmosphere avant l'advection en z |
C dans l'atmosphere avant l'advection en z |
83 |
sqi = 0. |
sqi = 0. |
251 |
C fin de la boucle principale sur les latitudes |
C fin de la boucle principale sur les latitudes |
252 |
C |
C |
253 |
1 CONTINUE |
1 CONTINUE |
|
C |
|
|
C------------------------------------------------------------- |
|
|
C |
|
|
C ----------- AA Test en fin de ADVX ------ Controle des S* |
|
|
|
|
|
c DO 9999 l = 1, llm |
|
|
c DO 9999 j = 1, jjp1 |
|
|
c DO 9999 i = 1, iip1 |
|
|
c IF (S0(i,j,l,ntra).lt.0..and.LIMIT) THEN |
|
|
c PRINT*, '-------------------' |
|
|
c PRINT*, 'En fin de ADVZ' |
|
|
c PRINT*,'S0(',i,j,l,')=',S0(i,j,l,ntra) |
|
|
c print*, 'sx(',i,j,l,')=',sx(i,j,l,ntra) |
|
|
c print*, 'sy(',i,j,l,')=',sy(i,j,l,ntra) |
|
|
c print*, 'sz(',i,j,l,')=',sz(i,j,l,ntra) |
|
|
c WRITE (*,*) 'On arrete !! - pbl en fin de ADVZ1' |
|
|
c STOP |
|
|
c ENDIF |
|
|
9999 CONTINUE |
|
254 |
|
|
255 |
C *** ------------------- bouclage cyclique en X ------------ |
C *** ------------------- bouclage cyclique en X ------------ |
256 |
|
|