42 |
|
|
43 |
print *, "Call sequence information: fxhyp" |
print *, "Call sequence information: fxhyp" |
44 |
|
|
|
xzoom = clon * pi_d / 180d0 |
|
|
|
|
|
IF (grossismx == 1.) THEN |
|
|
decalx = 1d0 |
|
|
else |
|
|
decalx = 0.75d0 |
|
|
END IF |
|
|
|
|
45 |
dzoom = dzoomx * twopi_d |
dzoom = dzoomx * twopi_d |
46 |
xtild = arth(- pi_d, pi_d / nmax, 2 * nmax + 1) |
xtild = arth(- pi_d, pi_d / nmax, 2 * nmax + 1) |
47 |
|
|
106 |
print *, "beta = ", beta |
print *, "beta = ", beta |
107 |
|
|
108 |
IF (2. * beta - grossismx <= 0.) THEN |
IF (2. * beta - grossismx <= 0.) THEN |
109 |
print *, 'Attention ! La valeur beta calculée dans fxhyp est mauvaise.' |
print *, 'Bad choice of grossismx, taux, dzoomx.' |
110 |
print *, 'Modifier les valeurs de grossismx, taux ou dzoomx et relancer.' |
print *, 'Decrease dzoomx or grossismx.' |
111 |
STOP 1 |
STOP 1 |
112 |
END IF |
END IF |
113 |
|
|
145 |
|
|
146 |
Xf(2 * nmax) = pi_d |
Xf(2 * nmax) = pi_d |
147 |
|
|
148 |
|
IF (grossismx == 1.) THEN |
149 |
|
decalx = 1d0 |
150 |
|
else |
151 |
|
decalx = 0.75d0 |
152 |
|
END IF |
153 |
|
|
154 |
|
xzoom = clon * pi_d / 180d0 |
155 |
call fxhyp_loop_ik(1, decalx, xf, xtild, Xprimt, xzoom, rlonm025, & |
call fxhyp_loop_ik(1, decalx, xf, xtild, Xprimt, xzoom, rlonm025, & |
156 |
xprimm025, xuv = - 0.25d0) |
xprimm025, xuv = - 0.25d0) |
157 |
call fxhyp_loop_ik(2, decalx, xf, xtild, Xprimt, xzoom, rlonv, xprimv, & |
call fxhyp_loop_ik(2, decalx, xf, xtild, Xprimt, xzoom, rlonv, xprimv, & |
164 |
print * |
print * |
165 |
|
|
166 |
forall (i = 1: iim) xlon(i) = rlonv(i + 1) - rlonv(i) |
forall (i = 1: iim) xlon(i) = rlonv(i + 1) - rlonv(i) |
167 |
print *, "Minimum longitude step:", MINval(xlon) * 180. / pi_d, "°" |
print *, "Minimum longitude step:", MINval(xlon) * 180. / pi_d, "degrees" |
168 |
print *, "Maximum longitude step:", MAXval(xlon) * 180. / pi_d, "°" |
print *, "Maximum longitude step:", MAXval(xlon) * 180. / pi_d, "degrees" |
169 |
|
|
170 |
DO i = 1, iim + 1 |
DO i = 1, iim + 1 |
171 |
IF (rlonp025(i) < rlonv(i)) THEN |
IF (rlonp025(i) < rlonv(i)) THEN |