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

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

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 119 by guez, Wed Jan 7 14:34:57 2015 UTC revision 121 by guez, Wed Jan 28 16:10:02 2015 UTC
# Line 124  contains Line 124  contains
124      ! Calcul des élongations cuij1, ..., cuij4, cvij1, ..., cvij4 aux mêmes      ! Calcul des élongations cuij1, ..., cuij4, cvij1, ..., cvij4 aux mêmes
125      ! endroits que les aires aireij1_2d, ..., aireij4_2d.      ! endroits que les aires aireij1_2d, ..., aireij4_2d.
126    
127      ! Choix entre une fonction "f(y)" à dérivée sinusoïdale ou à      ! Fonction "f(y)" à dérivée tangente hyperbolique. Calcul des
128      ! dérivée tangente hyperbolique. Calcul des coefficients cu_2d,      ! coefficients cu_2d, cv_2d, 1. / cu_2d**2, 1. / cv_2d**2. Les
129      ! cv_2d, 1. / cu_2d**2, 1. / cv_2d**2. Les coefficients cu_2d et cv_2d      ! coefficients cu_2d et cv_2d permettent de passer des vitesses
130      ! permettent de passer des vitesses naturelles aux vitesses      ! naturelles aux vitesses covariantes et contravariantes, ou
131      ! covariantes et contravariantes, ou vice-versa.      ! vice-versa.
132    
133      ! On a :      ! On a :
134      ! u(covariant) = cu_2d * u(naturel), u(contravariant) = u(naturel) / cu_2d      ! u(covariant) = cu_2d * u(naturel), u(contravariant) = u(naturel) / cu_2d
# Line 195  contains Line 195  contains
195      ELSE      ELSE
196         gamdi_gdiv = 0.         gamdi_gdiv = 0.
197      END IF      END IF
198    
199      IF (nitergrot/=2) THEN      IF (nitergrot/=2) THEN
200         gamdi_grot = coefdis / (real(nitergrot)-2.)         gamdi_grot = coefdis / (real(nitergrot)-2.)
201      ELSE      ELSE
202         gamdi_grot = 0.         gamdi_grot = 0.
203      END IF      END IF
204    
205      IF (niterh/=2) THEN      IF (niterh/=2) THEN
206         gamdi_h = coefdis / (real(niterh)-2.)         gamdi_h = coefdis / (real(niterh)-2.)
207      ELSE      ELSE
# Line 210  contains Line 212  contains
212      print *, "gamdi_grot = ", gamdi_grot      print *, "gamdi_grot = ", gamdi_grot
213      print *, "gamdi_h = ", gamdi_h      print *, "gamdi_h = ", gamdi_h
214    
     print *, 'inigeom: Y = latitude, dérivée tangente hyperbolique'  
215      CALL fyhyp(rlatu, yprimu, rlatv, rlatu2, yprimu2, rlatu1, yprimu1)      CALL fyhyp(rlatu, yprimu, rlatv, rlatu2, yprimu2, rlatu1, yprimu1)
216      CALL fxhyp(xprimm025, rlonv, xprimv, rlonu, xprimu, xprimp025)      CALL fxhyp(xprimm025, rlonv, xprimv, rlonu, xprimu, xprimp025)
217    
Legend:
Removed from v.119  
changed lines
  Added in v.121
  ViewVC Help
Powered by ViewVC 1.1.21