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

Diff of /trunk/dyn3d/advz.f

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

revision 3 by guez, Wed Feb 27 13:16:39 2008 UTC revision 32 by guez, Tue Apr 6 17:52:58 2010 UTC
# Line 33  C  w : flux de masse en z en Pa.m2.s-1 Line 33  C  w : flux de masse en z en Pa.m2.s-1
33        INTEGER ntra        INTEGER ntra
34        PARAMETER (ntra = 1)        PARAMETER (ntra = 1)
35    
36        REAL dtz        REAL, intent(in):: dtz
37        REAL w ( iip1,jjp1,llm )        REAL w ( iip1,jjp1,llm )
38            
39  C  moments: SM  total mass in each grid box  C  moments: SM  total mass in each grid box
# Line 78  C Line 78  C
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.
# Line 267  C Line 251  C
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                
Legend:
Removed from v.3  
changed lines
  Added in v.32
  ViewVC Help
Powered by ViewVC 1.1.21