/[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 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
# Line 215  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 267  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.96  
changed lines
  Added in v.97
  ViewVC Help
Powered by ViewVC 1.1.21