--- trunk/Sources/dyn3d/comgeom.f	2015/05/26 17:46:03	139
+++ trunk/Sources/dyn3d/comgeom.f	2016/03/11 18:47:26	178
@@ -111,11 +111,10 @@
     ! Calcul des élongations cuij1, ..., cuij4, cvij1, ..., cvij4 aux mêmes
     ! endroits que les aires aireij1_2d, ..., aireij4_2d.
 
-    ! 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.
+    ! 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
@@ -125,9 +124,6 @@
     ! u(covariant) = cu_2d * cu_2d * u(contravariant)
     ! v(covariant) = cv_2d * cv_2d * v(contravariant)
 
-    ! On a l'application (x(X), y(Y)) avec - im / 2 + 1 <= X <= im / 2
-    ! et - jm / 2 <= Y <= jm / 2
-
     ! x est la longitude du point en radians.
     ! y est la latitude du point en radians.
     ! 
@@ -139,23 +135,16 @@
     ! dépendant de j uniquement, sera ici indicé aussi en i pour un
     ! adressage plus facile en ij.
 
-    ! cvu et cv_2d sont respectivement les valeurs de
-    ! cv_2d aux points u et v.
-
-    ! cu_2d, cuv, cuscal, cuz sont respectivement les valeurs de cu_2d
-    ! aux points u, v, scalaires, et z. Cf. "inigeom.txt".
+    ! cv_2d est aux points v. cu_2d est aux points u. Cf. "inigeom.txt".
 
     USE comconst, ONLY : g, omeg, rad
     USE comdissnew, ONLY : coefdis, nitergdiv, nitergrot, niterh
     use dynetat0_m, only: xprimp025, xprimm025, rlatu1, rlatu2, rlatu, rlatv, &
-         yprimu1, yprimu2, rlonu, rlonv
-    use jumble, only: new_unit
-    use nr_util, only: pi
+         yprimu1, yprimu2
     USE paramet_m, ONLY : iip1, jjp1
 
     ! Local:
-    INTEGER i, j, unit
-    REAL cvu(iip1, jjp1), cuv(iip1, jjm)
+    INTEGER i, j
     REAL ai14, ai23, airez, un4rad2
     REAL coslatm, coslatp, radclatm, radclatp
     REAL, dimension(iip1, jjp1):: cuij1, cuij2, cuij3, cuij4 ! in m
@@ -163,9 +152,6 @@
     REAL gamdi_gdiv, gamdi_grot, gamdi_h
     real, dimension(iim + 1, jjm + 1):: aireij1_2d, aireij2_2d, aireij3_2d, &
          aireij4_2d ! in m2
-    real airuscv2_2d(iim + 1, jjm) 
-    real airvscu2_2d(iim + 1, jjm), aiuscv2gam_2d(iim + 1, jjm) 
-    real aivscu2gam_2d(iim + 1, jjm)
 
     !------------------------------------------------------------------
 
@@ -335,16 +321,12 @@
        unsairz_gam_2d(iip1, j) = unsairz_gam_2d(1, j)
     END DO
 
-    ! Calcul des élongations cu_2d, cv_2d, cvu 
+    ! Calcul des élongations cu_2d, cv_2d
 
     DO j = 1, jjm
        DO i = 1, iim
           cv_2d(i, j) = 0.5 * &
                (cvij2(i, j) + cvij3(i, j) + cvij1(i, j + 1) + cvij4(i, j + 1))
-          cvu(i, j) = 0.5 * (cvij1(i, j) + cvij4(i, j) + cvij2(i, j) &
-               + cvij3(i, j))
-          cuv(i, j) = 0.5 * (cuij2(i, j) + cuij3(i, j) + cuij1(i, j + 1) &
-               + cuij4(i, j + 1))
           unscv2_2d(i, j) = 1. / cv_2d(i, j)**2
        END DO
        DO i = 1, iim
@@ -355,9 +337,7 @@
           cvscuvgam_2d(i, j) = cvsurcuv_2d(i, j)**(-gamdi_grot)
        END DO
        cv_2d(iip1, j) = cv_2d(1, j)
-       cvu(iip1, j) = cvu(1, j)
        unscv2_2d(iip1, j) = unscv2_2d(1, j)
-       cuv(iip1, j) = cuv(1, j)
        cuvsurcv_2d(iip1, j) = cuvsurcv_2d(1, j)
        cvsurcuv_2d(iip1, j) = cvsurcuv_2d(1, j)
        cuvscvgam1_2d(iip1, j) = cuvscvgam1_2d(1, j)
@@ -389,29 +369,9 @@
 
     cu_2d(:, 1) = 0.
     unscu2_2d(:, 1) = 0.
-    cvu(:, 1) = 0.
 
     cu_2d(:, jjp1) = 0.
     unscu2_2d(:, jjp1) = 0.
-    cvu(:, jjp1) = 0.
-
-    DO j = 1, jjm
-       DO i = 1, iim
-          airvscu2_2d(i, j) = airev_2d(i, j) / (cuv(i, j) * cuv(i, j))
-          aivscu2gam_2d(i, j) = airvscu2_2d(i, j)**(-gamdi_grot)
-       END DO
-       airvscu2_2d(iip1, j) = airvscu2_2d(1, j)
-       aivscu2gam_2d(iip1, j) = aivscu2gam_2d(1, j)
-    END DO
-
-    DO j = 2, jjm
-       DO i = 1, iim
-          airuscv2_2d(i, j) = aireu_2d(i, j) / (cvu(i, j) * cvu(i, j))
-          aiuscv2gam_2d(i, j) = airuscv2_2d(i, j)**(-gamdi_grot)
-       END DO
-       airuscv2_2d(iip1, j) = airuscv2_2d(1, j)
-       aiuscv2gam_2d(iip1, j) = aiuscv2gam_2d(1, j)
-    END DO
 
     ! Calcul des aires aux pôles :
 
@@ -443,14 +403,6 @@
        constang_2d(iip1, j) = constang_2d(1, j)
     END DO
 
-    call new_unit(unit)
-    open(unit, file="longitude_latitude.txt", status="replace", action="write")
-    write(unit, fmt=*) '"longitudes at V points (degrees)"', rlonv * 180. / pi
-    write(unit, fmt=*) '"latitudes at V points (degrees)"', rlatv * 180. / pi
-    write(unit, fmt=*) '"longitudes at U points (degrees)"', rlonu * 180. / pi
-    write(unit, fmt=*) '"latitudes at U points (degrees)"', rlatu * 180. / pi
-    close(unit)
-
   END SUBROUTINE inigeom
 
 end module comgeom