--- trunk/dyn3d/comgeom.f	2014/03/05 14:57:53	82
+++ trunk/dyn3d/comgeom.f	2015/01/28 16:10:02	121
@@ -124,11 +124,11 @@
     ! Calcul des élongations cuij1, ..., cuij4, cvij1, ..., cvij4 aux mêmes
     ! endroits que les aires aireij1_2d, ..., aireij4_2d.
 
-    ! Choix entre une fonction "f(y)" à dérivée sinusoïdale ou à
-    ! dérivée tangente hyperbolique. Calcul des coefficients cu_2d,
-    ! cv_2d, 1. / cu_2d**2, 1. / cv_2d**2. Les coefficients cu_2d et cv_2d
-    ! permettent de passer des vitesses naturelles aux vitesses
-    ! covariantes et contravariantes, ou vice-versa.
+    ! Fonction "f(y)" à dérivée tangente hyperbolique. Calcul des
+    ! coefficients cu_2d, cv_2d, 1. / cu_2d**2, 1. / cv_2d**2. Les
+    ! coefficients cu_2d et cv_2d permettent de passer des vitesses
+    ! naturelles aux vitesses covariantes et contravariantes, ou
+    ! vice-versa.
 
     ! On a :
     ! u(covariant) = cu_2d * u(naturel), u(contravariant) = u(naturel) / cu_2d
@@ -163,29 +163,23 @@
 
     USE comconst, ONLY : g, omeg, rad
     USE comdissnew, ONLY : coefdis, nitergdiv, nitergrot, niterh
-    use conf_gcm_m, ONLY : fxyhypb, ysinus
-    use fxy_m, only: fxy
-    use fxyhyper_m, only: fxyhyper
+    use fxhyp_m, only: fxhyp
+    use fyhyp_m, only: fyhyp
     use jumble, only: new_unit
     use nr_util, only: pi
     USE paramet_m, ONLY : iip1, jjp1
-    USE serre, ONLY : alphax, alphay, clat, clon, dzoomx, dzoomy, grossismx, &
-         grossismy, pxo, pyo, taux, tauy, transx, transy
-    ! Modifies pxo, pyo, transx, transy
 
-    ! Variables locales
-
-    INTEGER i, j, itmax, itmay, iter, unit
+    ! Local:
+    INTEGER i, j, unit
     REAL cvu(iip1, jjp1), cuv(iip1, jjm)
     REAL ai14, ai23, airez, un4rad2
-    REAL eps, x1, xo1, f, df, xdm, y1, yo1, ydm
     REAL coslatm, coslatp, radclatm, radclatp
     REAL, dimension(iip1, jjp1):: cuij1, cuij2, cuij3, cuij4 ! in m
     REAL, dimension(iip1, jjp1):: cvij1, cvij2, cvij3, cvij4 ! in m
     REAL rlatu1(jjm), yprimu1(jjm), rlatu2(jjm), yprimu2(jjm)
-    real yprimv(jjm), yprimu(jjp1)
+    real yprimu(jjp1)
     REAL gamdi_gdiv, gamdi_grot, gamdi_h
-    REAL rlonm025(iip1), xprimm025(iip1), rlonp025(iip1), xprimp025(iip1)
+    REAL xprimm025(iip1), xprimp025(iip1)
     real, dimension(iim + 1, jjm + 1):: aireij1_2d, aireij2_2d, aireij3_2d, &
          aireij4_2d ! in m2
     real airuscv2_2d(iim + 1, jjm) 
@@ -201,11 +195,13 @@
     ELSE
        gamdi_gdiv = 0.
     END IF
+
     IF (nitergrot/=2) THEN
        gamdi_grot = coefdis / (real(nitergrot)-2.)
     ELSE
        gamdi_grot = 0.
     END IF
+
     IF (niterh/=2) THEN
        gamdi_h = coefdis / (real(niterh)-2.)
     ELSE
@@ -216,66 +212,8 @@
     print *, "gamdi_grot = ", gamdi_grot
     print *, "gamdi_h = ", gamdi_h
 
-    IF (.NOT. fxyhypb) THEN
-       IF (ysinus) THEN
-          print *, ' Inigeom, Y = Sinus (Latitude) '
-          ! utilisation de f(x, y) avec y = sinus de la latitude 
-          CALL fxysinus(rlatu, yprimu, rlatv, yprimv, rlatu1, yprimu1, &
-               rlatu2, yprimu2, rlonu, xprimu, rlonv, xprimv, rlonm025, &
-               xprimm025, rlonp025, xprimp025)
-       ELSE
-          print *, 'Inigeom, Y = Latitude, der. sinusoid .'
-          ! utilisation de f(x, y) a tangente sinusoidale, y etant la latit
-
-          pxo = clon * pi / 180.
-          pyo = 2. * clat * pi / 180.
-
-          ! determination de transx (pour le zoom) par Newton-Raphson
-
-          itmax = 10
-          eps = .1E-7
-
-          xo1 = 0.
-          DO iter = 1, itmax
-             x1 = xo1
-             f = x1 + alphax * sin(x1-pxo)
-             df = 1. + alphax * cos(x1-pxo)
-             x1 = x1 - f / df
-             xdm = abs(x1-xo1)
-             IF (xdm<=eps) EXIT
-             xo1 = x1
-          END DO
-
-          transx = xo1
-
-          itmay = 10
-          eps = .1E-7
-
-          yo1 = 0.
-          DO iter = 1, itmay
-             y1 = yo1
-             f = y1 + alphay * sin(y1-pyo)
-             df = 1. + alphay * cos(y1-pyo)
-             y1 = y1 - f / df
-             ydm = abs(y1-yo1)
-             IF (ydm<=eps) EXIT
-             yo1 = y1
-          END DO
-
-          transy = yo1
-
-          CALL fxy(rlatu, yprimu, rlatv, yprimv, rlatu1, yprimu1, rlatu2, &
-               yprimu2, rlonu, xprimu, rlonv, xprimv, rlonm025, xprimm025, &
-               rlonp025, xprimp025)
-       END IF
-    ELSE
-       ! Utilisation de fxyhyper, f(x, y) à dérivée tangente hyperbolique
-       print *, 'Inigeom, Y = Latitude, dérivée tangente hyperbolique'
-       CALL fxyhyper(clat, grossismy, dzoomy, tauy, clon, grossismx, dzoomx, &
-            taux, rlatu, yprimu, rlatv, yprimv, rlatu1, yprimu1, rlatu2, &
-            yprimu2, rlonu, xprimu, rlonv, xprimv, rlonm025, xprimm025, &
-            rlonp025, xprimp025)
-    END IF
+    CALL fyhyp(rlatu, yprimu, rlatv, rlatu2, yprimu2, rlatu1, yprimu1)
+    CALL fxhyp(xprimm025, rlonv, xprimv, rlonu, xprimu, xprimp025)
 
     rlatu(1) = pi / 2.
     rlatu(jjp1) = -rlatu(1)