--- trunk/dyn3d/comgeom.f	2015/01/07 14:34:57	119
+++ trunk/dyn3d/comgeom.f	2015/02/06 15:00:28	125
@@ -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
@@ -190,18 +190,20 @@
 
     PRINT *, 'Call sequence information: inigeom'
 
-    IF (nitergdiv/=2) THEN
-       gamdi_gdiv = coefdis / (real(nitergdiv)-2.)
+    IF (nitergdiv /= 2) THEN
+       gamdi_gdiv = coefdis / (nitergdiv - 2)
     ELSE
        gamdi_gdiv = 0.
     END IF
-    IF (nitergrot/=2) THEN
-       gamdi_grot = coefdis / (real(nitergrot)-2.)
+
+    IF (nitergrot /= 2) THEN
+       gamdi_grot = coefdis / (nitergrot - 2)
     ELSE
        gamdi_grot = 0.
     END IF
-    IF (niterh/=2) THEN
-       gamdi_h = coefdis / (real(niterh)-2.)
+
+    IF (niterh /= 2) THEN
+       gamdi_h = coefdis / (niterh - 2)
     ELSE
        gamdi_h = 0.
     END IF
@@ -210,7 +212,6 @@
     print *, "gamdi_grot = ", gamdi_grot
     print *, "gamdi_h = ", gamdi_h
 
-    print *, 'inigeom: Y = latitude, dérivée tangente hyperbolique'
     CALL fyhyp(rlatu, yprimu, rlatv, rlatu2, yprimu2, rlatu1, yprimu1)
     CALL fxhyp(xprimm025, rlonv, xprimv, rlonu, xprimu, xprimp025)