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

Diff of /trunk/dyn3d/comgeom.f

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

revision 82 by guez, Wed Mar 5 14:57:53 2014 UTC revision 97 by guez, Fri Apr 25 14:58:31 2014 UTC
# Line 166  contains Line 166  contains
166      use conf_gcm_m, ONLY : fxyhypb, ysinus      use conf_gcm_m, ONLY : fxyhypb, ysinus
167      use fxy_m, only: fxy      use fxy_m, only: fxy
168      use fxyhyper_m, only: fxyhyper      use fxyhyper_m, only: fxyhyper
169        use fxysinus_m, only: fxysinus
170      use jumble, only: new_unit      use jumble, only: new_unit
171      use nr_util, only: pi      use nr_util, only: pi
172      USE paramet_m, ONLY : iip1, jjp1      USE paramet_m, ONLY : iip1, jjp1
173      USE serre, ONLY : alphax, alphay, clat, clon, dzoomx, dzoomy, grossismx, &      USE serre, ONLY : alphax, alphay, clat, clon, dzoomx, dzoomy, grossismx, &
174           grossismy, pxo, pyo, taux, tauy, transx, transy           grossismy, pxo, pyo, taux, tauy, transx, transy
175      ! Modifies pxo, pyo, transx, transy      ! Modifiés pxo, pyo, transx, transy
   
     ! Variables locales  
176    
177        ! Local:
178      INTEGER i, j, itmax, itmay, iter, unit      INTEGER i, j, itmax, itmay, iter, unit
179      REAL cvu(iip1, jjp1), cuv(iip1, jjm)      REAL cvu(iip1, jjp1), cuv(iip1, jjm)
180      REAL ai14, ai23, airez, un4rad2      REAL ai14, ai23, airez, un4rad2
# Line 216  contains Line 216  contains
216      print *, "gamdi_grot = ", gamdi_grot      print *, "gamdi_grot = ", gamdi_grot
217      print *, "gamdi_h = ", gamdi_h      print *, "gamdi_h = ", gamdi_h
218    
219      IF (.NOT. fxyhypb) THEN      IF (fxyhypb) THEN
220           ! Utilisation de fxyhyper, f(x, y) à dérivée tangente hyperbolique
221           print *, 'inigeom: Y = latitude, dérivée tangente hyperbolique'
222           CALL fxyhyper(clat, grossismy, dzoomy, tauy, clon, grossismx, dzoomx, &
223                taux, rlatu, yprimu, rlatv, yprimv, rlatu1, yprimu1, rlatu2, &
224                yprimu2, rlonu, xprimu, rlonv, xprimv, rlonm025, xprimm025, &
225                rlonp025, xprimp025)
226        ELSE
227         IF (ysinus) THEN         IF (ysinus) THEN
228            print *, ' Inigeom, Y = Sinus (Latitude) '            print *, 'inigeom: Y = sin(latitude)'
229            ! utilisation de f(x, y) avec y = sinus de la latitude            ! Utilisation de f(x, y) avec y = sinus de la latitude
230            CALL fxysinus(rlatu, yprimu, rlatv, yprimv, rlatu1, yprimu1, &            CALL fxysinus(rlatu, yprimu, rlatv, yprimv, rlatu1, yprimu1, &
231                 rlatu2, yprimu2, rlonu, xprimu, rlonv, xprimv, rlonm025, &                 rlatu2, yprimu2, rlonu, xprimu, rlonv, xprimv, rlonm025, &
232                 xprimm025, rlonp025, xprimp025)                 xprimm025, rlonp025, xprimp025)
# Line 268  contains Line 275  contains
275                 yprimu2, rlonu, xprimu, rlonv, xprimv, rlonm025, xprimm025, &                 yprimu2, rlonu, xprimu, rlonv, xprimv, rlonm025, xprimm025, &
276                 rlonp025, xprimp025)                 rlonp025, xprimp025)
277         END IF         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)  
278      END IF      END IF
279    
280      rlatu(1) = pi / 2.      rlatu(1) = pi / 2.
Legend:
Removed from v.82  
changed lines
  Added in v.97
  ViewVC Help
Powered by ViewVC 1.1.21