/[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 96 by guez, Fri Apr 4 11:30:34 2014 UTC 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
          grossismy, pxo, pyo, taux, tauy, transx, transy  
174      ! Modifiés pxo, pyo, transx, transy      ! Modifiés pxo, pyo, transx, transy
175    
176      ! Local:      ! Local:
# Line 215  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 267  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.96  
changed lines
  Added in v.112
  ViewVC Help
Powered by ViewVC 1.1.21