/[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 167 by guez, Mon Aug 24 16:30:33 2015 UTC revision 170 by guez, Mon Sep 14 17:25:15 2015 UTC
# Line 23  contains Line 23  contains
23      ! Computes the eigenvalues and eigenvectors of the discrete analog      ! Computes the eigenvalues and eigenvectors of the discrete analog
24      ! of the second derivative with respect to longitude.      ! of the second derivative with respect to longitude.
25    
     use acc_m, only: acc  
26      USE dimens_m, ONLY: iim      USE dimens_m, ONLY: iim
27      USE dynetat0_m, ONLY: xprimu, xprimv      USE dynetat0_m, ONLY: xprimu, xprimv
28      use jumble, only: new_unit      use jumble, only: new_unit
# Line 62  contains Line 61  contains
61    
62      delta = matmul(deriv_v, deriv_u) ! second derivative at v longitudes      delta = matmul(deriv_v, deriv_u) ! second derivative at v longitudes
63      CALL jacobi(delta, eignval_v, eignfnv)      CALL jacobi(delta, eignval_v, eignfnv)
     CALL acc(eignfnv)  
64      CALL eigsrt(eignval_v, eignfnv)      CALL eigsrt(eignval_v, eignfnv)
65    
66      delta = matmul(deriv_u, deriv_v) ! second derivative at u longitudes      delta = matmul(deriv_u, deriv_v) ! second derivative at u longitudes
67      CALL jacobi(delta, eignval_u, eignfnu)      CALL jacobi(delta, eignval_u, eignfnu)
     CALL acc(eignfnu)  
68      CALL eigsrt(eignval_u, eignfnu)      CALL eigsrt(eignval_u, eignfnu)
69    
70      call new_unit(unit)      call new_unit(unit)
Legend:
Removed from v.167  
changed lines
  Added in v.170
  ViewVC Help
Powered by ViewVC 1.1.21