/[lmdze]/trunk/filtrez/filtreg_scal.f
ViewVC logotype

Diff of /trunk/filtrez/filtreg_scal.f

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

revision 30 by guez, Thu Apr 1 09:07:28 2010 UTC revision 32 by guez, Tue Apr 6 17:52:58 2010 UTC
# Line 13  contains Line 13  contains
13      ! pour l'opérateur filtre.      ! pour l'opérateur filtre.
14    
15      USE dimens_m, ONLY : iim, jjm      USE dimens_m, ONLY : iim, jjm
16      USE parafilt, ONLY : nfilun, nfilus, nfilvn, nfilvs      USE parafilt, ONLY: matriceun, matriceus, matricevn, matricevs, matrinvn, &
17             matrinvs
18      USE coefils, ONLY : jfiltnu, jfiltnv, jfiltsu, jfiltsv, sddu, sddv, &      USE coefils, ONLY : jfiltnu, jfiltnv, jfiltsu, jfiltsv, sddu, sddv, &
19           unsddu, unsddv           unsddu, unsddv
20    
# Line 48  contains Line 49  contains
49    
50      INTEGER jdfil1, jdfil2, jffil1, jffil2, jdfil, jffil      INTEGER jdfil1, jdfil2, jffil1, jffil2, jdfil, jffil
51      INTEGER i, j, l, k      INTEGER i, j, l, k
     REAL matriceun, matriceus, matricevn, matricevs, matrinvn, matrinvs  
     COMMON /matrfil/matriceun(iim, iim, nfilun), matriceus(iim, iim, nfilus), &  
          matricevn(iim, iim, nfilvn), matricevs(iim, iim, nfilvs), &  
          matrinvn(iim, iim, nfilun), matrinvs(iim, iim, nfilus)  
52      REAL eignq(iim), sdd1(iim), sdd2(iim)      REAL eignq(iim), sdd1(iim), sdd2(iim)
53      INTEGER hemisph      INTEGER hemisph
54    
# Line 156  contains Line 153  contains
153                     END DO                     END DO
154                     DO i = 1, iim                     DO i = 1, iim
155                        DO k = 1, iim                        DO k = 1, iim
156                           eignq(k) = eignq(k) + matriceun(k, i, j)*champ(i, j, l)                           eignq(k) = eignq(k) + matriceun(k, i, j) &
157                                  * champ(i, j, l)
158                        END DO                        END DO
159                     END DO                     END DO
160                  ELSE                  ELSE
# Line 165  contains Line 163  contains
163                     END DO                     END DO
164                     DO i = 1, iim                     DO i = 1, iim
165                        DO k = 1, iim                        DO k = 1, iim
166                           eignq(k) = eignq(k) + matricevn(k, i, j)*champ(i, j, l)                           eignq(k) = eignq(k) + matricevn(k, i, j) &
167                                  * champ(i, j, l)
168                        END DO                        END DO
169                     END DO                     END DO
170                  END IF                  END IF
Legend:
Removed from v.30  
changed lines
  Added in v.32
  ViewVC Help
Powered by ViewVC 1.1.21