7 |
! jfilts index of the first line filtered in SH |
! jfilts index of the first line filtered in SH |
8 |
|
|
9 |
! North: |
! North: |
10 |
|
|
11 |
real, allocatable:: matriceun(:, :, :), matrinvn(:, :, :) |
real, allocatable:: matriceun(:, :, :), matrinvn(:, :, :) |
12 |
! (iim, iim, 2:jfiltnu) |
! (iim, iim, 2:jfiltnu) |
13 |
|
|
14 |
real, allocatable:: matricevn(:, :, :) ! (iim, iim, jfiltnv) |
real, allocatable:: matricevn(:, :, :) ! (iim, iim, jfiltnv) |
15 |
|
|
16 |
! South: |
! South: |
17 |
|
|
18 |
real, allocatable:: matriceus(:, :, :), matrinvs(:, :, :) |
real, allocatable:: matriceus(:, :, :), matrinvs(:, :, :) |
19 |
! (iim, iim, jfiltsu:jjm) |
! (iim, iim, jfiltsu:jjm) |
20 |
|
|
27 |
! From filtrez/inifilr.F, version 1.1.1.1 2004/05/19 12:53:09 |
! From filtrez/inifilr.F, version 1.1.1.1 2004/05/19 12:53:09 |
28 |
! H. Upadhyaya, O. Sharma |
! H. Upadhyaya, O. Sharma |
29 |
|
|
30 |
! This routine computes the eigenfunctions of the laplacian on the |
! This routine computes the eigenvectors of the laplacian on the |
31 |
! stretched grid, and the filtering coefficients. The modes are |
! stretched grid, and the filtering coefficients. The modes are |
32 |
! filtered from modfrst to iim. |
! filtered from modfrst to iim. |
33 |
|
|
47 |
REAL eignft(iim, iim), coff |
REAL eignft(iim, iim), coff |
48 |
|
|
49 |
real eignfnu(iim, iim), eignfnv(iim, iim) |
real eignfnu(iim, iim), eignfnv(iim, iim) |
50 |
! eigenfunctions of the discrete laplacian |
! eigenvectors of the discrete laplacian |
51 |
|
|
52 |
! Filtering coefficients (lamda_max * cos(rlat) / lamda): |
! Filtering coefficients (lamda_max * cos(rlat) / lamda): |
53 |
real coefilu(iim, jjm), coefilv(iim, jjm) |
real coefilu(iim, jjm), coefilv(iim, jjm) |
63 |
|
|
64 |
CALL inifgn(eignvl, eignfnu, eignfnv) |
CALL inifgn(eignvl, eignfnu, eignfnv) |
65 |
|
|
66 |
! compute eigenvalues and eigenfunctions |
! compute eigenvalues and eigenvectors |
67 |
! compute the filtering coefficients for scalar lines and |
! compute the filtering coefficients for scalar lines and |
68 |
! meridional wind v-lines |
! meridional wind v-lines |
69 |
! we filter all those latitude lines where coefil < 1 |
! we filter all those latitude lines where coefil < 1 |