25 |
|
|
26 |
REAL, intent(in):: dzooma ! distance totale de la zone du zoom |
REAL, intent(in):: dzooma ! distance totale de la zone du zoom |
27 |
|
|
28 |
REAL, intent(in):: tau |
REAL, intent(in):: tau |
29 |
! raideur de la transition de l'intérieur à l'extérieur du zoom |
! raideur de la transition de l'intérieur à l'extérieur du zoom |
30 |
|
|
31 |
! arguments de sortie |
! arguments de sortie |
58 |
INTEGER is2 |
INTEGER is2 |
59 |
SAVE is2 |
SAVE is2 |
60 |
|
|
|
DOUBLE PRECISION heavyside |
|
|
|
|
61 |
!---------------------------------------------------------------------- |
!---------------------------------------------------------------------- |
62 |
|
|
63 |
pi = 2. * ASIN(1.) |
pi = 2. * ASIN(1.) |
80 |
STOP 1 |
STOP 1 |
81 |
ELSE |
ELSE |
82 |
dzoom = dzooma * pi/180. |
dzoom = dzooma * pi/180. |
83 |
ENDIF |
END IF |
84 |
|
|
85 |
print *, ' xzoom(rad), grossism, tau, dzoom (rad):' |
print *, ' xzoom(rad), grossism, tau, dzoom (rad):' |
86 |
print *, xzoom, grossism, tau, dzoom |
print *, xzoom, grossism, tau, dzoom |
87 |
|
|
88 |
DO i = 0, nmax2 |
DO i = 0, nmax2 |
89 |
xtild(i) = - pi + FLOAT(i) * depi /nmax2 |
xtild(i) = - pi + REAL(i) * depi /nmax2 |
90 |
ENDDO |
ENDDO |
91 |
|
|
92 |
DO i = nmax, nmax2 |
DO i = nmax, nmax2 |
107 |
ELSE |
ELSE |
108 |
fhyp (i) = TANH (fa/fb) |
fhyp (i) = TANH (fa/fb) |
109 |
ENDIF |
ENDIF |
110 |
ENDIF |
END IF |
111 |
|
|
112 |
IF (xtild(i) == 0.) fhyp(i) = 1. |
IF (xtild(i) == 0.) fhyp(i) = 1. |
113 |
IF (xtild(i) == pi) fhyp(i) = -1. |
IF (xtild(i) == pi) fhyp(i) = -1. |
114 |
ENDDO |
END DO |
115 |
|
|
116 |
! Calcul de beta |
! Calcul de beta |
117 |
|
|
150 |
print *, 'Attention ! La valeur beta calculée dans fxhyp est mauvaise.' |
print *, 'Attention ! La valeur beta calculée dans fxhyp est mauvaise.' |
151 |
print *, 'Modifier les valeurs de grossismx, tau ou dzoomx et relancer.' |
print *, 'Modifier les valeurs de grossismx, tau ou dzoomx et relancer.' |
152 |
STOP 1 |
STOP 1 |
153 |
ENDIF |
END IF |
154 |
|
|
155 |
! calcul de Xprimt |
! calcul de Xprimt |
156 |
|
|
157 |
DO i = nmax, nmax2 |
DO i = nmax, nmax2 |
158 |
Xprimt(i) = beta + (grossism - beta) * fhyp(i) |
Xprimt(i) = beta + (grossism - beta) * fhyp(i) |
159 |
ENDDO |
END DO |
160 |
|
|
161 |
DO i = nmax + 1, nmax2 |
DO i = nmax + 1, nmax2 |
162 |
Xprimt(nmax2 - i) = Xprimt(i) |
Xprimt(nmax2 - i) = Xprimt(i) |
163 |
ENDDO |
END DO |
164 |
|
|
165 |
! Calcul de Xf |
! Calcul de Xf |
166 |
|
|
218 |
ENDIF |
ENDIF |
219 |
|
|
220 |
DO i = ii1, ii2 |
DO i = ii1, ii2 |
221 |
xlon2 = - pi + (FLOAT(i) + xuv - decalx) * depi / FLOAT(iim) |
xlon2 = - pi + (REAL(i) + xuv - decalx) * depi / REAL(iim) |
222 |
Xfi = xlon2 |
Xfi = xlon2 |
223 |
|
|
224 |
it = nmax2 |
it = nmax2 |
264 |
end if |
end if |
265 |
|
|
266 |
|
|
267 |
xxprim(i) = depi/ (FLOAT(iim) * Xprimin) |
xxprim(i) = depi/ (REAL(iim) * Xprimin) |
268 |
xvrai(i) = xi + xzoom |
xvrai(i) = xi + xzoom |
269 |
end DO |
end DO |
270 |
|
|