18 |
! 1., taux=0., clon=0.) est à - 180 degrés. |
! 1., taux=0., clon=0.) est à - 180 degrés. |
19 |
|
|
20 |
USE dimens_m, ONLY: iim |
USE dimens_m, ONLY: iim |
21 |
use fxhyp_loop_ik_m, only: fxhyp_loop_ik, nmax |
use dynetat0_m, only: clon, grossismx, dzoomx, taux |
22 |
|
use invert_zoom_x_m, only: invert_zoom_x, nmax |
23 |
use nr_util, only: pi, pi_d, twopi, twopi_d, arth |
use nr_util, only: pi, pi_d, twopi, twopi_d, arth |
24 |
use principal_cshift_m, only: principal_cshift |
use principal_cshift_m, only: principal_cshift |
|
use serre, only: clon, grossismx, dzoomx, taux |
|
25 |
|
|
26 |
REAL, intent(out):: xprimm025(:), rlonv(:), xprimv(:) ! (iim + 1) |
REAL, intent(out):: xprimm025(:), rlonv(:), xprimv(:) ! (iim + 1) |
27 |
real, intent(out):: rlonu(:), xprimu(:), xprimp025(:) ! (iim + 1) |
real, intent(out):: rlonu(:), xprimu(:), xprimp025(:) ! (iim + 1) |
53 |
rlonm025(:iim) = rlonv(:iim) - 0.25 * step |
rlonm025(:iim) = rlonv(:iim) - 0.25 * step |
54 |
rlonp025(:iim) = rlonv(:iim) + 0.25 * step |
rlonp025(:iim) = rlonv(:iim) + 0.25 * step |
55 |
rlonu(:iim) = rlonv(:iim) + 0.5 * step |
rlonu(:iim) = rlonv(:iim) + 0.5 * step |
56 |
else |
else test_grossismx |
57 |
dzoom = dzoomx * twopi_d |
dzoom = dzoomx * twopi_d |
58 |
xtild = arth(- pi_d, pi_d / nmax, 2 * nmax + 1) |
xtild = arth(- pi_d, pi_d / nmax, 2 * nmax + 1) |
59 |
|
|
157 |
|
|
158 |
Xf(2 * nmax) = pi_d |
Xf(2 * nmax) = pi_d |
159 |
|
|
160 |
call fxhyp_loop_ik(xf, xtild, Xprimt, rlonm025(:iim), xprimm025(:iim), & |
call invert_zoom_x(xf, xtild, Xprimt, rlonm025(:iim), xprimm025(:iim), & |
161 |
xuv = - 0.25d0) |
xuv = - 0.25d0) |
162 |
call fxhyp_loop_ik(xf, xtild, Xprimt, rlonv(:iim), xprimv(:iim), & |
call invert_zoom_x(xf, xtild, Xprimt, rlonv(:iim), xprimv(:iim), & |
163 |
xuv = 0d0) |
xuv = 0d0) |
164 |
call fxhyp_loop_ik(xf, xtild, Xprimt, rlonu(:iim), xprimu(:iim), & |
call invert_zoom_x(xf, xtild, Xprimt, rlonu(:iim), xprimu(:iim), & |
165 |
xuv = 0.5d0) |
xuv = 0.5d0) |
166 |
call fxhyp_loop_ik(xf, xtild, Xprimt, rlonp025(:iim), xprimp025(:iim), & |
call invert_zoom_x(xf, xtild, Xprimt, rlonp025(:iim), xprimp025(:iim), & |
167 |
xuv = 0.25d0) |
xuv = 0.25d0) |
168 |
end if test_grossismx |
end if test_grossismx |
169 |
|
|