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

Diff of /trunk/Sources/filtrez/inifgn.f

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

revision 142 by guez, Fri Jun 5 18:58:06 2015 UTC revision 143 by guez, Tue Jun 9 14:32:46 2015 UTC
# Line 10  module inifgn_m Line 10  module inifgn_m
10    real unsddu(iim), unsddv(iim)    real unsddu(iim), unsddv(iim)
11    
12    real eignfnu(iim, iim), eignfnv(iim, iim)    real eignfnu(iim, iim), eignfnv(iim, iim)
13    ! eignfn eigenfunctions of the discrete laplacian    ! eigenfunctions of the discrete laplacian
14    
15  contains  contains
16    
# Line 18  contains Line 18  contains
18    
19      ! From LMDZ4/libf/filtrez/inifgn.F, v 1.1.1.1 2004/05/19 12:53:09      ! From LMDZ4/libf/filtrez/inifgn.F, v 1.1.1.1 2004/05/19 12:53:09
20    
21      ! H.Upadyaya, O.Sharma      ! H. Upadyaya, O. Sharma
22    
23        use acc_m, only: acc
24      USE dimens_m, ONLY: iim      USE dimens_m, ONLY: iim
25      USE dynetat0_m, ONLY: xprimu, xprimv      USE dynetat0_m, ONLY: xprimu, xprimv
26      use nr_util, only: pi      use nr_util, only: pi
27      use numer_rec_95, only: jacobi      use numer_rec_95, only: jacobi, eigsrt
28    
29      real, intent(out):: dv(iim)      real, intent(out):: dv(:) ! (iim)
30    
31      ! Local:      ! Local:
32      REAL vec(iim, iim), vec1(iim, iim)      REAL vec(iim, iim), vec1(iim, iim)
33      REAL du(iim)      REAL du(iim)
     real d(iim)  
34      INTEGER i, j, k, nrot      INTEGER i, j, k, nrot
35    
     EXTERNAL acc  
   
36      !----------------------------------------------------------------      !----------------------------------------------------------------
37    
38      sddv = sqrt(xprimv(:iim))      sddv = sqrt(xprimv(:iim))
# Line 82  contains Line 80  contains
80      END DO      END DO
81    
82      CALL jacobi(vec, dv, eignfnv, nrot)      CALL jacobi(vec, dv, eignfnv, nrot)
83      CALL acc(eignfnv, d, iim)      CALL acc(eignfnv)
84      CALL eigen_sort(dv, eignfnv, iim, iim)      CALL eigsrt(dv, eignfnv)
85    
86      CALL jacobi(vec1, du, eignfnu, nrot)      CALL jacobi(vec1, du, eignfnu, nrot)
87      CALL acc(eignfnu, d, iim)      CALL acc(eignfnu)
88      CALL eigen_sort(du, eignfnu, iim, iim)      CALL eigsrt(du, eignfnu)
89    
90    END SUBROUTINE inifgn    END SUBROUTINE inifgn
91    
Legend:
Removed from v.142  
changed lines
  Added in v.143
  ViewVC Help
Powered by ViewVC 1.1.21