/[lmdze]/trunk/Sources/dyn3d/comgeom.f

Diff of /trunk/Sources/dyn3d/comgeom.f

Parent Directory | Revision Log | View Patch Patch

-revision 139 by guez,
Tue May 26 17:46:03 2015 UTC
+revision 178 by guez,
Fri Mar 11 18:47:26 2016 UTC
 Line 111 
 contains
      ! 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
+     ! Calcul des coefficients cu_2d, cv_2d, 1. / cu_2d**2, 1. /
-     ! coefficients cu_2d, cv_2d, 1. / cu_2d**2, 1. / cv_2d**2. Les
+     ! cv_2d**2. Les coefficients cu_2d et cv_2d permettent de passer
-     ! coefficients cu_2d et cv_2d permettent de passer des vitesses
+     ! des vitesses naturelles aux vitesses covariantes et
-     ! naturelles aux vitesses covariantes et contravariantes, ou
+     ! contravariantes, ou vice-versa.
-     ! vice-versa.
      ! On a :
      ! u(covariant) = cu_2d * u(naturel), u(contravariant) = u(naturel) / cu_2d
-Line 125 
 contains
+Line 124 
 contains
      ! 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.
      !
-Line 139 
 contains
+Line 135 
 contains
      ! 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 est aux points v. cu_2d est aux points u. Cf. "inigeom.txt".
-     ! 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".
      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
+          yprimu1, yprimu2
-     use jumble, only: new_unit
-     use nr_util, only: pi
      USE paramet_m, ONLY : iip1, jjp1
      ! Local:
-     INTEGER i, j, unit
+     INTEGER i, j
-     REAL cvu(iip1, jjp1), cuv(iip1, jjm)
      REAL ai14, ai23, airez, un4rad2
      REAL coslatm, coslatp, radclatm, radclatp
      REAL, dimension(iip1, jjp1):: cuij1, cuij2, cuij3, cuij4 ! in m
-Line 163 
 contains
+Line 152 
 contains
      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)
      !------------------------------------------------------------------
-Line 335 
 contains
+Line 321 
 contains
         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
-Line 355 
 contains
+Line 337 
 contains
            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)
-Line 389 
 contains
+Line 369 
 contains
      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 :
-Line 443 
 contains
+Line 403 
 contains
         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

 Legend:



Removed from v.139
 


changed lines


 
Added in v.178
 Legend:



Removed from v.139
 


changed lines


 
Added in v.178
-Removed from v.139
+Added in v.178

	ViewVC Help
Powered by ViewVC 1.1.21