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

Diff of /trunk/dyn3d/comgeom.f

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

trunk/dyn3d/comgeom.f90 revision 79 by guez, Fri Feb 28 17:52:47 2014 UTC trunk/dyn3d/comgeom.f revision 112 by guez, Thu Sep 18 13:36:51 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, pxo, pyo, transx, transy
174           grossismy, pxo, pyo, taux, tauy, transx, transy      ! Modifiés pxo, pyo, transx, transy
     ! Modifies pxo, pyo, transx, transy  
   
     ! Variables locales  
175    
176        ! Local:
177      INTEGER i, j, itmax, itmay, iter, unit      INTEGER i, j, itmax, itmay, iter, unit
178      REAL cvu(iip1, jjp1), cuv(iip1, jjm)      REAL cvu(iip1, jjp1), cuv(iip1, jjm)
179      REAL ai14, ai23, airez, un4rad2      REAL ai14, ai23, airez, un4rad2
# Line 216  contains Line 215  contains
215      print *, "gamdi_grot = ", gamdi_grot      print *, "gamdi_grot = ", gamdi_grot
216      print *, "gamdi_h = ", gamdi_h      print *, "gamdi_h = ", gamdi_h
217    
218      IF (.NOT. fxyhypb) THEN      IF (fxyhypb) THEN
219           print *, 'inigeom: Y = latitude, dérivée tangente hyperbolique'
220           CALL fxyhyper(rlatu, yprimu, rlatv, yprimv, rlatu1, yprimu1, rlatu2, &
221                yprimu2, rlonu, xprimu, rlonv, xprimv, rlonm025, xprimm025, &
222                rlonp025, xprimp025)
223        ELSE
224         IF (ysinus) THEN         IF (ysinus) THEN
225            print *, ' Inigeom, Y = Sinus (Latitude) '            print *, 'inigeom: Y = sin(latitude)'
226            ! utilisation de f(x, y) avec y = sinus de la latitude            ! Utilisation de f(x, y) avec y = sinus de la latitude
227            CALL fxysinus(rlatu, yprimu, rlatv, yprimv, rlatu1, yprimu1, &            CALL fxysinus(rlatu, yprimu, rlatv, yprimv, rlatu1, yprimu1, &
228                 rlatu2, yprimu2, rlonu, xprimu, rlonv, xprimv, rlonm025, &                 rlatu2, yprimu2, rlonu, xprimu, rlonv, xprimv, rlonm025, &
229                 xprimm025, rlonp025, xprimp025)                 xprimm025, rlonp025, xprimp025)
# Line 268  contains Line 272  contains
272                 yprimu2, rlonu, xprimu, rlonv, xprimv, rlonm025, xprimm025, &                 yprimu2, rlonu, xprimu, rlonv, xprimv, rlonm025, xprimm025, &
273                 rlonp025, xprimp025)                 rlonp025, xprimp025)
274         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)  
275      END IF      END IF
276    
277      rlatu(1) = pi / 2.      rlatu(1) = pi / 2.
Legend:
Removed from v.79  
changed lines
  Added in v.112
  ViewVC Help
Powered by ViewVC 1.1.21