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

Diff of /trunk/filtrez/inifgn.f

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

revision 131 by guez, Wed Jan 28 16:10:02 2015 UTC revision 132 by guez, Fri Mar 20 16:31:06 2015 UTC
# Line 10  contains Line 10  contains
10    
11      ! H.Upadyaya, O.Sharma      ! H.Upadyaya, O.Sharma
12    
     USE dimens_m, ONLY: iim  
     USE comgeom, ONLY: xprimu, xprimv  
13      USE coefils, ONLY: eignfnu, eignfnv, sddu, sddv, unsddu, unsddv      USE coefils, ONLY: eignfnu, eignfnv, sddu, sddv, unsddu, unsddv
14        USE comgeom, ONLY: xprimu, xprimv
15        USE dimens_m, ONLY: iim
16        use numer_rec_95, only: jacobi
17    
18      real dv(iim)      real, intent(out):: dv(iim)
19    
20      ! Local:      ! Local:
21      REAL vec(iim, iim), vec1(iim, iim)      REAL vec(iim, iim), vec1(iim, iim)
# Line 23  contains Line 24  contains
24      REAL pi      REAL pi
25      INTEGER i, j, k, imm1, nrot      INTEGER i, j, k, imm1, nrot
26    
27      EXTERNAL acc, jacobi      EXTERNAL acc
28    
29      !----------------------------------------------------------------      !----------------------------------------------------------------
30    
# Line 74  contains Line 75  contains
75         END DO         END DO
76      END DO      END DO
77    
78      CALL jacobi(vec, iim, iim, dv, eignfnv, nrot)      CALL jacobi(vec, dv, eignfnv, nrot)
79      CALL acc(eignfnv, d, iim)      CALL acc(eignfnv, d, iim)
80      CALL eigen_sort(dv, eignfnv, iim, iim)      CALL eigen_sort(dv, eignfnv, iim, iim)
81    
82      CALL jacobi(vec1, iim, iim, du, eignfnu, nrot)      CALL jacobi(vec1, du, eignfnu, nrot)
83      CALL acc(eignfnu, d, iim)      CALL acc(eignfnu, d, iim)
84      CALL eigen_sort(du, eignfnu, iim, iim)      CALL eigen_sort(du, eignfnu, iim, iim)
85    
Legend:
Removed from v.131  
changed lines
  Added in v.132
  ViewVC Help
Powered by ViewVC 1.1.21