33 |
DOUBLE PRECISION xtild(0:2 * nmax) |
DOUBLE PRECISION xtild(0:2 * nmax) |
34 |
DOUBLE PRECISION fhyp(nmax:2 * nmax), ffdx, beta, Xprimt(0:2 * nmax) |
DOUBLE PRECISION fhyp(nmax:2 * nmax), ffdx, beta, Xprimt(0:2 * nmax) |
35 |
DOUBLE PRECISION Xf(0:2 * nmax), xxpr(2 * nmax) |
DOUBLE PRECISION Xf(0:2 * nmax), xxpr(2 * nmax) |
36 |
DOUBLE PRECISION xzoom, fa, fb |
DOUBLE PRECISION fa, fb |
37 |
INTEGER i, is2 |
INTEGER i, is2 |
38 |
DOUBLE PRECISION xmoy, fxm |
DOUBLE PRECISION xmoy, fxm |
39 |
|
|
41 |
|
|
42 |
print *, "Call sequence information: fxhyp" |
print *, "Call sequence information: fxhyp" |
43 |
|
|
|
xzoom = clon * pi_d / 180d0 |
|
|
|
|
44 |
test_grossismx: if (grossismx == 1.) then |
test_grossismx: if (grossismx == 1.) then |
45 |
step = twopi / iim |
step = twopi / iim |
46 |
|
|
49 |
xprimv(:iim) = step |
xprimv(:iim) = step |
50 |
xprimu(:iim) = step |
xprimu(:iim) = step |
51 |
|
|
52 |
rlonv(:iim) = arth(- pi + clon * pi / 180., step, iim) |
rlonv(:iim) = arth(- pi + clon, step, iim) |
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 |
157 |
|
|
158 |
Xf(2 * nmax) = pi_d |
Xf(2 * nmax) = pi_d |
159 |
|
|
160 |
call fxhyp_loop_ik(xf, xtild, Xprimt, xzoom, rlonm025(:iim), & |
call fxhyp_loop_ik(xf, xtild, Xprimt, rlonm025(:iim), xprimm025(:iim), & |
161 |
xprimm025(:iim), xuv = - 0.25d0) |
xuv = - 0.25d0) |
162 |
call fxhyp_loop_ik(xf, xtild, Xprimt, xzoom, rlonv(:iim), & |
call fxhyp_loop_ik(xf, xtild, Xprimt, rlonv(:iim), xprimv(:iim), & |
163 |
xprimv(:iim), xuv = 0d0) |
xuv = 0d0) |
164 |
call fxhyp_loop_ik(xf, xtild, Xprimt, xzoom, rlonu(:iim), & |
call fxhyp_loop_ik(xf, xtild, Xprimt, rlonu(:iim), xprimu(:iim), & |
165 |
xprimu(:iim), xuv = 0.5d0) |
xuv = 0.5d0) |
166 |
call fxhyp_loop_ik(xf, xtild, Xprimt, xzoom, rlonp025(:iim), & |
call fxhyp_loop_ik(xf, xtild, Xprimt, rlonp025(:iim), xprimp025(:iim), & |
167 |
xprimp025(:iim), xuv = 0.25d0) |
xuv = 0.25d0) |
168 |
end if test_grossismx |
end if test_grossismx |
169 |
|
|
170 |
is2 = 0 |
is2 = 0 |