57 |
|
|
58 |
REAL dlatu(jjm) |
REAL dlatu(jjm) |
59 |
REAL rlamda(2:iim) ! > 0, in descending order |
REAL rlamda(2:iim) ! > 0, in descending order |
60 |
real eignvl(iim) ! eigenvalues sorted in descending order (<= 0) |
real eignvl(iim) ! eigenvalues (<= 0) sorted in descending order |
61 |
INTEGER j, unit |
INTEGER j, unit |
62 |
REAL colat0 ! > 0 |
REAL colat0 ! > 0 |
63 |
integer j1 ! index of smallest positive latitude |
integer j1 ! index of smallest positive latitude |
77 |
PRINT *, 'colat0 = ', colat0 |
PRINT *, 'colat0 = ', colat0 |
78 |
|
|
79 |
rlamda = iim / (pi * colat0 / grossismx) / sqrt(- eignvl(2: iim)) |
rlamda = iim / (pi * colat0 / grossismx) / sqrt(- eignvl(2: iim)) |
80 |
|
print *, "1 / rlamda(iim) = ", 1. / rlamda(iim) |
81 |
call new_unit(unit) |
call new_unit(unit) |
82 |
open(unit, file = "modfrst.csv", status = "replace", action = "write") |
open(unit, file = "modfrst.csv", status = "replace", action = "write") |
83 |
write(unit, fmt = *) '"rlat (degrees)" modfrst' ! title line |
write(unit, fmt = *) '"rlat (degrees)" modfrst' ! title line |
84 |
|
|
|
! D\'etermination de jfilt[ns][uv] : |
|
|
|
|
85 |
j1 = ifirstloc(rlatu <= 0.) |
j1 = ifirstloc(rlatu <= 0.) |
86 |
|
|
87 |
call inifilr_hemisph(rlatu(j1 - 1:2:- 1), colat0, rlamda, unit, eignfnv, & |
call inifilr_hemisph(rlatu(j1 - 1:2:- 1), colat0, rlamda, unit, eignfnv, & |