31 |
|
|
32 |
USE dimens_m, ONLY : iim, jjm |
USE dimens_m, ONLY : iim, jjm |
33 |
USE dynetat0_m, ONLY : rlatu, rlatv, xprimu, grossismx |
USE dynetat0_m, ONLY : rlatu, rlatv, xprimu, grossismx |
34 |
use inifgn_m, only: inifgn, eignfnu, eignfnv |
use inifgn_m, only: inifgn |
35 |
use jumble, only: new_unit |
use jumble, only: new_unit |
36 |
use nr_util, only: pi |
use nr_util, only: pi |
37 |
|
|
44 |
REAL colat0 ! > 0 |
REAL colat0 ! > 0 |
45 |
REAL eignft(iim, iim), coff |
REAL eignft(iim, iim), coff |
46 |
|
|
47 |
|
real eignfnu(iim, iim), eignfnv(iim, iim) |
48 |
|
! eigenfunctions of the discrete laplacian |
49 |
|
|
50 |
! Filtering coefficients (lamda_max * cos(rlat) / lamda): |
! Filtering coefficients (lamda_max * cos(rlat) / lamda): |
51 |
real coefilu(iim, jjm), coefilv(iim, jjm) |
real coefilu(iim, jjm), coefilv(iim, jjm) |
52 |
real coefilu2(iim, jjm), coefilv2(iim, jjm) |
real coefilu2(iim, jjm), coefilv2(iim, jjm) |
59 |
|
|
60 |
print *, "Call sequence information: inifilr" |
print *, "Call sequence information: inifilr" |
61 |
|
|
62 |
CALL inifgn(eignvl) |
CALL inifgn(eignvl, eignfnu, eignfnv) |
63 |
|
|
64 |
! compute eigenvalues and eigenfunctions |
! compute eigenvalues and eigenfunctions |
65 |
! compute the filtering coefficients for scalar lines and |
! compute the filtering coefficients for scalar lines and |