6 |
|
|
7 |
private iim |
private iim |
8 |
|
|
9 |
real sddu(iim), sddv(iim) ! SQRT(dx) |
real sddu(iim), sddv(iim) ! SQRT(dx / di) |
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) |
26 |
use nr_util, only: pi |
use nr_util, only: pi |
27 |
use numer_rec_95, only: jacobi, eigsrt |
use numer_rec_95, only: jacobi, eigsrt |
28 |
|
|
29 |
real, intent(out):: dv(:) ! (iim) |
real, intent(out):: dv(:) ! (iim) eigenvalues sorted in descending order |
30 |
|
|
31 |
! Local: |
! Local: |
32 |
REAL vec(iim, iim), vec1(iim, iim) |
REAL vec(iim, iim), vec1(iim, iim) |
35 |
|
|
36 |
!---------------------------------------------------------------- |
!---------------------------------------------------------------- |
37 |
|
|
38 |
|
print *, "Call sequence information: inifgn" |
39 |
|
|
40 |
sddv = sqrt(xprimv(:iim)) |
sddv = sqrt(xprimv(:iim)) |
41 |
sddu = sqrt(xprimu(:iim)) |
sddu = sqrt(xprimu(:iim)) |
42 |
unsddu = 1. / sddu |
unsddu = 1. / sddu |