Sources/dyn3d/fyhyp.f

module fyhyp_m

  IMPLICIT NONE

contains

  SUBROUTINE fyhyp(rlatu, rlatv, rlatu2, yprimu2, rlatu1, yprimu1)

    ! From LMDZ4/libf/dyn3d/fyhyp.F, version 1.2, 2005/06/03 09:11:32

    ! Author: P. Le Van, from analysis by R. Sadourny 

    ! Calcule les latitudes et dérivées dans la grille du GCM pour une
    ! fonction f(y) à dérivée tangente hyperbolique.

    ! Il vaut mieux avoir : grossismy * dzoom < pi / 2

    use coefpoly_m, only: coefpoly, a0, a1, a2, a3
    USE dimens_m, only: jjm
    use dynetat0_m, only: clat, grossismy, dzoomy, tauy
    use heavyside_m, only: heavyside

    REAL, intent(out):: rlatu(:) ! (jjm + 1)
    REAL, intent(out):: rlatv(:) ! (jjm)
    real, intent(out):: rlatu2(:), yprimu2(:), rlatu1(:), yprimu1(:) ! (jjm)

    ! Local: 

    INTEGER, PARAMETER:: nmax=30000, nmax2=2*nmax
    REAL dzoom ! distance totale de la zone du zoom (en radians)
    DOUBLE PRECISION ylat(jjm + 1), yprim(jjm + 1)
    DOUBLE PRECISION yuv
    DOUBLE PRECISION, save:: yt(0:nmax2)
    DOUBLE PRECISION fhyp(0:nmax2), beta
    DOUBLE PRECISION, save:: ytprim(0:nmax2)
    DOUBLE PRECISION fxm(0:nmax2)
    DOUBLE PRECISION, save:: yf(0:nmax2)
    DOUBLE PRECISION yypr(0:nmax2)
    DOUBLE PRECISION yvrai(jjm + 1), yprimm(jjm + 1), ylatt(jjm + 1)
    DOUBLE PRECISION pi, pis2, epsilon, pisjm
    DOUBLE PRECISION yo1, yi, ylon2, ymoy, yprimin
    DOUBLE PRECISION yfi, yf1, ffdy
    DOUBLE PRECISION ypn
    DOUBLE PRECISION, save::deply, y00

    INTEGER i, j, it, ik, iter, jlat, jjpn
    INTEGER, save:: jpn
    DOUBLE PRECISION yi2, heavyy0, heavyy0m
    DOUBLE PRECISION fa(0:nmax2), fb(0:nmax2)
    REAL y0min, y0max

    !-------------------------------------------------------------------

    print *, "Call sequence information: fyhyp"

    pi = 2.*asin(1.)
    pis2 = pi/2.
    pisjm = pi/real(jjm)
    epsilon = 1e-3
    dzoom = dzoomy*pi
    print *, 'yzoom(rad), grossismy, tauy, dzoom (rad):'
    print *, clat, grossismy, tauy, dzoom

    DO i = 0, nmax2
       yt(i) = -pis2 + real(i)*pi/nmax2
    END DO

    heavyy0m = heavyside(-clat)
    heavyy0 = heavyside(clat)
    y0min = 2.*clat*heavyy0m - pis2
    y0max = 2.*clat*heavyy0 + pis2

    fa = 999.999
    fb = 999.999

    DO i = 0, nmax2
       IF (yt(i)<clat) THEN
          fa(i) = tauy*(yt(i)-clat + dzoom/2.)
          fb(i) = (yt(i)-2.*clat*heavyy0m + pis2)*(clat-yt(i))
       ELSE IF (yt(i)>clat) THEN
          fa(i) = tauy*(clat-yt(i) + dzoom/2.)
          fb(i) = (2.*clat*heavyy0-yt(i) + pis2)*(yt(i)-clat)
       END IF

       IF (200.*fb(i)<-fa(i)) THEN
          fhyp(i) = -1.
       ELSE IF (200.*fb(i)<fa(i)) THEN
          fhyp(i) = 1.
       ELSE
          fhyp(i) = tanh(fa(i)/fb(i))
       END IF

       IF (yt(i)==clat) fhyp(i) = 1.
       IF (yt(i)==y0min .OR. yt(i)==y0max) fhyp(i) = -1.
    END DO

    ! Calcul de beta 

    ffdy = 0.

    DO i = 1, nmax2
       ymoy = 0.5*(yt(i-1) + yt(i))
       IF (ymoy<clat) THEN
          fa(i) = tauy*(ymoy-clat + dzoom/2.)
          fb(i) = (ymoy-2.*clat*heavyy0m + pis2)*(clat-ymoy)
       ELSE IF (ymoy>clat) THEN
          fa(i) = tauy*(clat-ymoy + dzoom/2.)
          fb(i) = (2.*clat*heavyy0-ymoy + pis2)*(ymoy-clat)
       END IF

       IF (200.*fb(i)<-fa(i)) THEN
          fxm(i) = -1.
       ELSE IF (200.*fb(i)<fa(i)) THEN
          fxm(i) = 1.
       ELSE
          fxm(i) = tanh(fa(i)/fb(i))
       END IF
       IF (ymoy==clat) fxm(i) = 1.
       IF (ymoy==y0min .OR. yt(i)==y0max) fxm(i) = -1.
       ffdy = ffdy + fxm(i)*(yt(i)-yt(i-1))
    END DO

    beta = (grossismy*ffdy-pi)/(ffdy-pi)

    IF (2. * beta - grossismy <= 0.) THEN
       print *, 'Attention ! La valeur beta calculee dans la routine fyhyp ' &
            // 'est mauvaise. Modifier les valeurs de grossismy, tauy ou ' &
            // 'dzoomy et relancer.'
       STOP 1
    END IF

    ! calcul de Ytprim 

    DO i = 0, nmax2
       ytprim(i) = beta + (grossismy-beta)*fhyp(i)
    END DO

    ! Calcul de Yf 

    yf(0) = -pis2
    DO i = 1, nmax2
       yypr(i) = beta + (grossismy-beta)*fxm(i)
    END DO

    DO i = 1, nmax2
       yf(i) = yf(i-1) + yypr(i)*(yt(i)-yt(i-1))
    END DO

    ! yuv = 0. si calcul des latitudes aux pts. U 
    ! yuv = 0.5 si calcul des latitudes aux pts. V 

    loop_ik: DO ik = 1, 4
       IF (ik==1) THEN
          yuv = 0.
          jlat = jjm + 1
       ELSE IF (ik==2) THEN
          yuv = 0.5
          jlat = jjm
       ELSE IF (ik==3) THEN
          yuv = 0.25
          jlat = jjm
       ELSE IF (ik==4) THEN
          yuv = 0.75
          jlat = jjm
       END IF

       yo1 = 0.
       DO j = 1, jlat
          yo1 = 0.
          ylon2 = -pis2 + pisjm*(real(j) + yuv-1.)
          yfi = ylon2

          it = nmax2
          DO while (it >= 1 .and. yfi < yf(it))
             it = it - 1
          END DO

          yi = yt(it)
          IF (it==nmax2) THEN
             it = nmax2 - 1
             yf(it + 1) = pis2
          END IF

          ! Interpolation entre yi(it) et yi(it + 1) pour avoir Y(yi)
          ! et Y'(yi) 

          CALL coefpoly(yf(it), yf(it + 1), ytprim(it), ytprim(it + 1), &
               yt(it), yt(it + 1))

          yf1 = yf(it)
          yprimin = a1 + 2.*a2*yi + 3.*a3*yi*yi

          iter = 1
          DO
             yi = yi - (yf1-yfi)/yprimin
             IF (abs(yi-yo1)<=epsilon .or. iter == 300) exit
             yo1 = yi
             yi2 = yi*yi
             yf1 = a0 + a1*yi + a2*yi2 + a3*yi2*yi
             yprimin = a1 + 2.*a2*yi + 3.*a3*yi2
          END DO
          if (abs(yi-yo1) > epsilon) then
             print *, 'Pas de solution.', j, ylon2
             STOP 1
          end if

          yprimin = a1 + 2.*a2*yi + 3.*a3*yi*yi
          yprim(j) = pi/(jjm*yprimin)
          yvrai(j) = yi
       END DO

       DO j = 1, jlat - 1
          IF (yvrai(j + 1)<yvrai(j)) THEN
             print *, 'Problème avec rlat(', j + 1, ') plus petit que rlat(', &
                  j, ')'
             STOP 1
          END IF
       END DO

       print *, 'Reorganisation des latitudes pour avoir entre - pi/2 et pi/2'

       IF (ik==1) THEN
          ypn = pis2
          DO j = jjm + 1, 1, -1
             IF (yvrai(j)<=ypn) exit
          END DO

          jpn = j
          y00 = yvrai(jpn)
          deply = pis2 - y00
       END IF

       DO j = 1, jjm + 1 - jpn
          ylatt(j) = -pis2 - y00 + yvrai(jpn + j-1)
          yprimm(j) = yprim(jpn + j-1)
       END DO

       jjpn = jpn
       IF (jlat==jjm) jjpn = jpn - 1

       DO j = 1, jjpn
          ylatt(j + jjm + 1-jpn) = yvrai(j) + deply
          yprimm(j + jjm + 1-jpn) = yprim(j)
       END DO

       ! Fin de la reorganisation

       DO j = 1, jlat
          ylat(j) = ylatt(jlat + 1-j)
          yprim(j) = yprimm(jlat + 1-j)
       END DO

       DO j = 1, jlat
          yvrai(j) = ylat(j)*180./pi
       END DO

       IF (ik==1) THEN
          DO j = 1, jjm + 1
             rlatu(j) = ylat(j)
          END DO
       ELSE IF (ik==2) THEN
          DO j = 1, jjm
             rlatv(j) = ylat(j)
          END DO
       ELSE IF (ik==3) THEN
          DO j = 1, jjm
             rlatu2(j) = ylat(j)
             yprimu2(j) = yprim(j)
          END DO
       ELSE IF (ik==4) THEN
          DO j = 1, jjm
             rlatu1(j) = ylat(j)
             yprimu1(j) = yprim(j)
          END DO
       END IF
    END DO loop_ik

    DO j = 1, jjm
       ylat(j) = rlatu(j) - rlatu(j + 1)
    END DO

    DO j = 1, jjm
       IF (rlatu1(j) <= rlatu2(j)) THEN
          print *, 'Attention ! rlatu1 < rlatu2 ', rlatu1(j), rlatu2(j), j
          STOP 13
       ENDIF

       IF (rlatu2(j) <= rlatu(j+1)) THEN
          print *, 'Attention ! rlatu2 < rlatup1 ', rlatu2(j), rlatu(j+1), j
          STOP 14
       ENDIF

       IF (rlatu(j) <= rlatu1(j)) THEN
          print *, ' Attention ! rlatu < rlatu1 ', rlatu(j), rlatu1(j), j
          STOP 15
       ENDIF

       IF (rlatv(j) <= rlatu2(j)) THEN
          print *, ' Attention ! rlatv < rlatu2 ', rlatv(j), rlatu2(j), j
          STOP 16
       ENDIF

       IF (rlatv(j) >= rlatu1(j)) THEN
          print *, ' Attention ! rlatv > rlatu1 ', rlatv(j), rlatu1(j), j
          STOP 17
       ENDIF

       IF (rlatv(j) >= rlatu(j)) THEN
          print *, ' Attention ! rlatv > rlatu ', rlatv(j), rlatu(j), j
          STOP 18
       ENDIF
    ENDDO

    print *, 'Latitudes'
    print 3, minval(ylat(:jjm)) *180d0/pi, maxval(ylat(:jjm))*180d0/pi

3   Format(1x, ' Au centre du zoom, la longueur de la maille est', &
         ' d environ ', f0.2, ' degres ', /, &
         ' alors que la maille en dehors de la zone du zoom est ', &
         "d'environ ", f0.2, ' degres ')

  END SUBROUTINE fyhyp

end module fyhyp_m
1	guez	97	module fyhyp_m
2	guez	3
3	guez	112	IMPLICIT NONE
4	guez	3
5	guez	97	contains
6	guez	3
7	guez	139	SUBROUTINE fyhyp(rlatu, rlatv, rlatu2, yprimu2, rlatu1, yprimu1)
8	guez	3
9	guez	97	! From LMDZ4/libf/dyn3d/fyhyp.F, version 1.2, 2005/06/03 09:11:32
10	guez	3
11	guez	97	! Author: P. Le Van, from analysis by R. Sadourny
12	guez	3
13	guez	97	! Calcule les latitudes et dérivées dans la grille du GCM pour une
14	guez	119	! fonction f(y) à dérivée tangente hyperbolique.
15	guez	3
16	guez	121	! Il vaut mieux avoir : grossismy * dzoom < pi / 2
17	guez	3
18	guez	149	use coefpoly_m, only: coefpoly, a0, a1, a2, a3
19	guez	97	USE dimens_m, only: jjm
20	guez	139	use dynetat0_m, only: clat, grossismy, dzoomy, tauy
21	guez	131	use heavyside_m, only: heavyside
22	guez	3
23	guez	212	REAL, intent(out):: rlatu(:) ! (jjm + 1)
24			REAL, intent(out):: rlatv(:) ! (jjm)
25			real, intent(out):: rlatu2(:), yprimu2(:), rlatu1(:), yprimu1(:) ! (jjm)
26	guez	3
27	guez	97	! Local:
28	guez	3
29	guez	97	INTEGER, PARAMETER:: nmax=30000, nmax2=2*nmax
30			REAL dzoom ! distance totale de la zone du zoom (en radians)
31	guez	118	DOUBLE PRECISION ylat(jjm + 1), yprim(jjm + 1)
32	guez	97	DOUBLE PRECISION yuv
33			DOUBLE PRECISION, save:: yt(0:nmax2)
34			DOUBLE PRECISION fhyp(0:nmax2), beta
35			DOUBLE PRECISION, save:: ytprim(0:nmax2)
36			DOUBLE PRECISION fxm(0:nmax2)
37			DOUBLE PRECISION, save:: yf(0:nmax2)
38			DOUBLE PRECISION yypr(0:nmax2)
39	guez	118	DOUBLE PRECISION yvrai(jjm + 1), yprimm(jjm + 1), ylatt(jjm + 1)
40	guez	131	DOUBLE PRECISION pi, pis2, epsilon, pisjm
41	guez	97	DOUBLE PRECISION yo1, yi, ylon2, ymoy, yprimin
42			DOUBLE PRECISION yfi, yf1, ffdy
43	guez	212	DOUBLE PRECISION ypn
44			DOUBLE PRECISION, save::deply, y00
45	guez	3
46	guez	212	INTEGER i, j, it, ik, iter, jlat, jjpn
47			INTEGER, save:: jpn
48	guez	149	DOUBLE PRECISION yi2, heavyy0, heavyy0m
49	guez	97	DOUBLE PRECISION fa(0:nmax2), fb(0:nmax2)
50			REAL y0min, y0max
51	guez	3
52	guez	97	!-------------------------------------------------------------------
53	guez	3
54	guez	120	print *, "Call sequence information: fyhyp"
55
56	guez	97	pi = 2.*asin(1.)
57			pis2 = pi/2.
58			pisjm = pi/real(jjm)
59			epsilon = 1e-3
60	guez	121	dzoom = dzoomy*pi
61	guez	119	print *, 'yzoom(rad), grossismy, tauy, dzoom (rad):'
62	guez	131	print *, clat, grossismy, tauy, dzoom
63	guez	81
64	guez	97	DO i = 0, nmax2
65			yt(i) = -pis2 + real(i)*pi/nmax2
66			END DO
67	guez	81
68	guez	131	heavyy0m = heavyside(-clat)
69			heavyy0 = heavyside(clat)
70			y0min = 2.clatheavyy0m - pis2
71			y0max = 2.clatheavyy0 + pis2
72	guez	81
73	guez	97	fa = 999.999
74			fb = 999.999
75	guez	81
76	guez	97	DO i = 0, nmax2
77	guez	131	IF (yt(i)<clat) THEN
78			fa(i) = tauy*(yt(i)-clat + dzoom/2.)
79			fb(i) = (yt(i)-2.clatheavyy0m + pis2)*(clat-yt(i))
80			ELSE IF (yt(i)>clat) THEN
81			fa(i) = tauy*(clat-yt(i) + dzoom/2.)
82			fb(i) = (2.clatheavyy0-yt(i) + pis2)*(yt(i)-clat)
83	guez	97	END IF
84	guez	81
85	guez	97	IF (200.*fb(i)<-fa(i)) THEN
86			fhyp(i) = -1.
87			ELSE IF (200.*fb(i)<fa(i)) THEN
88			fhyp(i) = 1.
89			ELSE
90			fhyp(i) = tanh(fa(i)/fb(i))
91			END IF
92	guez	81
93	guez	131	IF (yt(i)==clat) fhyp(i) = 1.
94	guez	97	IF (yt(i)==y0min .OR. yt(i)==y0max) fhyp(i) = -1.
95			END DO
96	guez	81
97	guez	97	! Calcul de beta
98	guez	81
99	guez	97	ffdy = 0.
100	guez	81
101	guez	97	DO i = 1, nmax2
102			ymoy = 0.5*(yt(i-1) + yt(i))
103	guez	131	IF (ymoy<clat) THEN
104			fa(i) = tauy*(ymoy-clat + dzoom/2.)
105			fb(i) = (ymoy-2.clatheavyy0m + pis2)*(clat-ymoy)
106			ELSE IF (ymoy>clat) THEN
107			fa(i) = tauy*(clat-ymoy + dzoom/2.)
108			fb(i) = (2.clatheavyy0-ymoy + pis2)*(ymoy-clat)
109	guez	97	END IF
110	guez	81
111	guez	97	IF (200.*fb(i)<-fa(i)) THEN
112			fxm(i) = -1.
113			ELSE IF (200.*fb(i)<fa(i)) THEN
114			fxm(i) = 1.
115			ELSE
116			fxm(i) = tanh(fa(i)/fb(i))
117			END IF
118	guez	131	IF (ymoy==clat) fxm(i) = 1.
119	guez	97	IF (ymoy==y0min .OR. yt(i)==y0max) fxm(i) = -1.
120			ffdy = ffdy + fxm(i)*(yt(i)-yt(i-1))
121			END DO
122	guez	81
123	guez	119	beta = (grossismy*ffdy-pi)/(ffdy-pi)
124	guez	81
125	guez	119	IF (2. * beta - grossismy <= 0.) THEN
126	guez	97	print *, 'Attention ! La valeur beta calculee dans la routine fyhyp ' &
127			// 'est mauvaise. Modifier les valeurs de grossismy, tauy ou ' &
128			// 'dzoomy et relancer.'
129			STOP 1
130	guez	81	END IF
131
132	guez	97	! calcul de Ytprim
133	guez	81
134	guez	97	DO i = 0, nmax2
135	guez	119	ytprim(i) = beta + (grossismy-beta)*fhyp(i)
136	guez	97	END DO
137	guez	81
138	guez	97	! Calcul de Yf
139	guez	81
140	guez	97	yf(0) = -pis2
141			DO i = 1, nmax2
142	guez	119	yypr(i) = beta + (grossismy-beta)*fxm(i)
143	guez	81	END DO
144	guez	3
145	guez	97	DO i = 1, nmax2
146			yf(i) = yf(i-1) + yypr(i)*(yt(i)-yt(i-1))
147	guez	81	END DO
148	guez	3
149	guez	97	! yuv = 0. si calcul des latitudes aux pts. U
150			! yuv = 0.5 si calcul des latitudes aux pts. V
151	guez	3
152	guez	97	loop_ik: DO ik = 1, 4
153			IF (ik==1) THEN
154			yuv = 0.
155			jlat = jjm + 1
156			ELSE IF (ik==2) THEN
157			yuv = 0.5
158			jlat = jjm
159			ELSE IF (ik==3) THEN
160			yuv = 0.25
161			jlat = jjm
162			ELSE IF (ik==4) THEN
163			yuv = 0.75
164			jlat = jjm
165			END IF
166	guez	3
167	guez	97	yo1 = 0.
168			DO j = 1, jlat
169			yo1 = 0.
170			ylon2 = -pis2 + pisjm*(real(j) + yuv-1.)
171			yfi = ylon2
172	guez	3
173	guez	97	it = nmax2
174			DO while (it >= 1 .and. yfi < yf(it))
175			it = it - 1
176			END DO
177	guez	3
178	guez	97	yi = yt(it)
179			IF (it==nmax2) THEN
180			it = nmax2 - 1
181			yf(it + 1) = pis2
182			END IF
183	guez	3
184	guez	97	! Interpolation entre yi(it) et yi(it + 1) pour avoir Y(yi)
185			! et Y'(yi)
186	guez	3
187	guez	97	CALL coefpoly(yf(it), yf(it + 1), ytprim(it), ytprim(it + 1), &
188	guez	149	yt(it), yt(it + 1))
189	guez	3
190	guez	97	yf1 = yf(it)
191			yprimin = a1 + 2.a2yi + 3.a3yi*yi
192	guez	3
193	guez	97	iter = 1
194			DO
195			yi = yi - (yf1-yfi)/yprimin
196			IF (abs(yi-yo1)<=epsilon .or. iter == 300) exit
197			yo1 = yi
198			yi2 = yi*yi
199			yf1 = a0 + a1yi + a2yi2 + a3yi2yi
200			yprimin = a1 + 2.a2yi + 3.a3yi2
201			END DO
202			if (abs(yi-yo1) > epsilon) then
203			print *, 'Pas de solution.', j, ylon2
204			STOP 1
205			end if
206	guez	3
207	guez	97	yprimin = a1 + 2.a2yi + 3.a3yi*yi
208			yprim(j) = pi/(jjm*yprimin)
209			yvrai(j) = yi
210			END DO
211	guez	3
212	guez	97	DO j = 1, jlat - 1
213			IF (yvrai(j + 1)<yvrai(j)) THEN
214			print *, 'Problème avec rlat(', j + 1, ') plus petit que rlat(', &
215			j, ')'
216			STOP 1
217			END IF
218			END DO
219	guez	3
220	guez	97	print *, 'Reorganisation des latitudes pour avoir entre - pi/2 et pi/2'
221	guez	3
222	guez	97	IF (ik==1) THEN
223			ypn = pis2
224	guez	118	DO j = jjm + 1, 1, -1
225	guez	97	IF (yvrai(j)<=ypn) exit
226			END DO
227	guez	3
228	guez	97	jpn = j
229			y00 = yvrai(jpn)
230			deply = pis2 - y00
231			END IF
232	guez	3
233	guez	97	DO j = 1, jjm + 1 - jpn
234			ylatt(j) = -pis2 - y00 + yvrai(jpn + j-1)
235			yprimm(j) = yprim(jpn + j-1)
236			END DO
237	guez	3
238	guez	97	jjpn = jpn
239			IF (jlat==jjm) jjpn = jpn - 1
240	guez	3
241	guez	97	DO j = 1, jjpn
242			ylatt(j + jjm + 1-jpn) = yvrai(j) + deply
243			yprimm(j + jjm + 1-jpn) = yprim(j)
244			END DO
245	guez	3
246	guez	97	! Fin de la reorganisation
247	guez	3
248	guez	97	DO j = 1, jlat
249			ylat(j) = ylatt(jlat + 1-j)
250			yprim(j) = yprimm(jlat + 1-j)
251			END DO
252	guez	81
253	guez	97	DO j = 1, jlat
254			yvrai(j) = ylat(j)*180./pi
255			END DO
256	guez	81
257	guez	97	IF (ik==1) THEN
258	guez	118	DO j = 1, jjm + 1
259	guez	119	rlatu(j) = ylat(j)
260	guez	97	END DO
261			ELSE IF (ik==2) THEN
262	guez	118	DO j = 1, jjm
263	guez	119	rlatv(j) = ylat(j)
264	guez	97	END DO
265			ELSE IF (ik==3) THEN
266	guez	118	DO j = 1, jjm
267	guez	97	rlatu2(j) = ylat(j)
268			yprimu2(j) = yprim(j)
269			END DO
270			ELSE IF (ik==4) THEN
271	guez	118	DO j = 1, jjm
272	guez	97	rlatu1(j) = ylat(j)
273			yprimu1(j) = yprim(j)
274			END DO
275			END IF
276			END DO loop_ik
277	guez	81
278	guez	97	DO j = 1, jjm
279	guez	119	ylat(j) = rlatu(j) - rlatu(j + 1)
280	guez	97	END DO
281	guez	81
282	guez	119	DO j = 1, jjm
283			IF (rlatu1(j) <= rlatu2(j)) THEN
284			print *, 'Attention ! rlatu1 < rlatu2 ', rlatu1(j), rlatu2(j), j
285			STOP 13
286			ENDIF
287
288			IF (rlatu2(j) <= rlatu(j+1)) THEN
289			print *, 'Attention ! rlatu2 < rlatup1 ', rlatu2(j), rlatu(j+1), j
290			STOP 14
291			ENDIF
292
293			IF (rlatu(j) <= rlatu1(j)) THEN
294			print *, ' Attention ! rlatu < rlatu1 ', rlatu(j), rlatu1(j), j
295			STOP 15
296			ENDIF
297
298			IF (rlatv(j) <= rlatu2(j)) THEN
299			print *, ' Attention ! rlatv < rlatu2 ', rlatv(j), rlatu2(j), j
300			STOP 16
301			ENDIF
302
303			IF (rlatv(j) >= rlatu1(j)) THEN
304			print *, ' Attention ! rlatv > rlatu1 ', rlatv(j), rlatu1(j), j
305			STOP 17
306			ENDIF
307
308			IF (rlatv(j) >= rlatu(j)) THEN
309			print *, ' Attention ! rlatv > rlatu ', rlatv(j), rlatu(j), j
310			STOP 18
311			ENDIF
312			ENDDO
313
314			print *, 'Latitudes'
315	guez	212	print 3, minval(ylat(:jjm)) 180d0/pi, maxval(ylat(:jjm))180d0/pi
316	guez	119
317			3 Format(1x, ' Au centre du zoom, la longueur de la maille est', &
318			' d environ ', f0.2, ' degres ', /, &
319			' alors que la maille en dehors de la zone du zoom est ', &
320	guez	120	"d'environ ", f0.2, ' degres ')
321	guez	119
322	guez	97	END SUBROUTINE fyhyp
323	guez	81
324	guez	97	end module fyhyp_m