109 |
END IF |
END IF |
110 |
END DO |
END DO |
111 |
|
|
112 |
IF (.not. (MINval(xvrai(:iim)) >= - pi_d - 1d-5 & |
IF (ik == 1 .and. (MINval(xvrai(:iim)) < - pi_d - 0.1d0 & |
113 |
.and. MAXval(xvrai(:iim)) <= pi_d + 1d-5)) THEN |
.or. MAXval(xvrai(:iim)) > pi_d + 0.1d0)) THEN |
114 |
IF (xzoom <= 0.) THEN |
IF (xzoom <= 0.) THEN |
115 |
IF (ik == 1) THEN |
i = 1 |
|
i = 1 |
|
116 |
|
|
117 |
do while (xvrai(i) < - pi_d .and. i < iim) |
do while (xvrai(i) < - pi_d .and. i < iim) |
118 |
i = i + 1 |
i = i + 1 |
119 |
end do |
end do |
120 |
|
|
121 |
if (xvrai(i) < - pi_d) then |
if (xvrai(i) < - pi_d) then |
122 |
print *, 'Xvrai plus petit que - pi !' |
print *, 'Xvrai plus petit que - pi !' |
123 |
STOP 1 |
STOP 1 |
124 |
end if |
end if |
125 |
|
|
126 |
is2 = i |
is2 = i |
127 |
END IF |
ELSE |
128 |
|
i = iim |
129 |
|
|
130 |
|
do while (xvrai(i) > pi_d .and. i > 1) |
131 |
|
i = i - 1 |
132 |
|
end do |
133 |
|
|
134 |
|
if (xvrai(i) > pi_d) then |
135 |
|
print *, 'Xvrai plus grand que pi !' |
136 |
|
STOP 1 |
137 |
|
end if |
138 |
|
|
139 |
|
is2 = i |
140 |
|
END IF |
141 |
|
END IF |
142 |
|
|
143 |
|
if (is2 /= 0) then |
144 |
|
IF (xzoom <= 0.) THEN |
145 |
IF (is2 /= 1) THEN |
IF (is2 /= 1) THEN |
146 |
DO ii = is2, iim |
DO ii = is2, iim |
147 |
xlon(ii-is2 + 1) = xvrai(ii) |
xlon(ii-is2 + 1) = xvrai(ii) |
152 |
xprimm(ii + iim-is2 + 1) = xxprim(ii) |
xprimm(ii + iim-is2 + 1) = xxprim(ii) |
153 |
END DO |
END DO |
154 |
END IF |
END IF |
155 |
ELSE |
else |
|
IF (ik == 1) THEN |
|
|
i = iim |
|
|
|
|
|
do while (xvrai(i) > pi_d .and. i > 1) |
|
|
i = i - 1 |
|
|
end do |
|
|
|
|
|
if (xvrai(i) > pi_d) then |
|
|
print *, 'Xvrai plus grand que pi !' |
|
|
STOP 1 |
|
|
end if |
|
|
|
|
|
is2 = i |
|
|
END IF |
|
|
|
|
156 |
idif = iim -is2 |
idif = iim -is2 |
157 |
|
|
158 |
DO ii = 1, is2 |
DO ii = 1, is2 |
164 |
xlon(ii) = xvrai(ii + is2) - twopi_d |
xlon(ii) = xvrai(ii + is2) - twopi_d |
165 |
xprimm(ii) = xxprim(ii + is2) |
xprimm(ii) = xxprim(ii + is2) |
166 |
END DO |
END DO |
167 |
END IF |
end IF |
168 |
END IF |
end if |
169 |
|
|
170 |
xlon(iim + 1) = xlon(1) + twopi |
xlon(iim + 1) = xlon(1) + twopi |
171 |
xprimm(iim + 1) = xprimm(1) |
xprimm(iim + 1) = xprimm(1) |