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

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

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

revision 146 by guez, Tue Jun 16 17:27:33 2015 UTC revision 147 by guez, Wed Jun 17 14:20:14 2015 UTC
# Line 30  contains Line 30  contains
30      ! Local:      ! Local:
31      real rlonm025(iim + 1), rlonp025(iim + 1), d_rlonv(iim)      real rlonm025(iim + 1), rlonp025(iim + 1), d_rlonv(iim)
32      REAL dzoom, step      REAL dzoom, step
33      DOUBLE PRECISION, dimension(0:nmax):: xtild, fhyp, G, Xf      DOUBLE PRECISION, dimension(0:nmax):: xtild, fhyp, G, Xf, ffdx
34      DOUBLE PRECISION ffdx, beta      DOUBLE PRECISION beta
35      INTEGER i, is2      INTEGER i, is2
36      DOUBLE PRECISION xxpr(nmax - 1), xmoy(nmax), fxm(nmax)      DOUBLE PRECISION xmoy(nmax), fxm(nmax)
37    
38      !----------------------------------------------------------------------      !----------------------------------------------------------------------
39    
# Line 57  contains Line 57  contains
57         forall (i = 1:nmax) xmoy(i) = 0.5d0 * (xtild(i-1) + xtild(i))         forall (i = 1:nmax) xmoy(i) = 0.5d0 * (xtild(i-1) + xtild(i))
58    
59         ! Compute fhyp:         ! Compute fhyp:
60         fhyp(1:nmax - 1) = tanh_cautious(taux * (dzoom / 2. &         fhyp(1:nmax - 1) = tanh_cautious(taux * (dzoom / 2d0 &
61              - xtild(1:nmax - 1)), xtild(1:nmax - 1) &              - xtild(1:nmax - 1)), xtild(1:nmax - 1) &
62              * (pi_d - xtild(1:nmax - 1)))              * (pi_d - xtild(1:nmax - 1)))
63         fhyp(0) = 1d0         fhyp(0) = 1d0
64         fhyp(nmax) = -1d0         fhyp(nmax) = -1d0
65    
66         fxm = tanh_cautious(taux * (dzoom / 2. - xmoy), xmoy * (pi_d - xmoy))         fxm = tanh_cautious(taux * (dzoom / 2d0 - xmoy), xmoy * (pi_d - xmoy))
67    
68         ! Calcul de beta         ! Compute \int_0 ^{\tilde x} F:
69    
70         ffdx = 0.         ffdx(0) = 0d0
71    
72         DO i = 1, nmax         DO i = 1, nmax
73            ffdx = ffdx + fxm(i) * (xtild(i) - xtild(i-1))            ffdx(i) = ffdx(i - 1) + fxm(i) * (xtild(i) - xtild(i-1))
74         END DO         END DO
75    
76         print *, "ffdx = ", ffdx         print *, "ffdx(nmax) = ", ffdx(nmax)
77         beta = (pi_d - grossismx * ffdx) / (pi_d - ffdx)         beta = (pi_d - grossismx * ffdx(nmax)) / (pi_d - ffdx(nmax))
78         print *, "beta = ", beta         print *, "beta = ", beta
79    
80         IF (2. * beta - grossismx <= 0.) THEN         IF (2d0 * beta - grossismx <= 0d0) THEN
81            print *, 'Bad choice of grossismx, taux, dzoomx.'            print *, 'Bad choice of grossismx, taux, dzoomx.'
82            print *, 'Decrease dzoomx or grossismx.'            print *, 'Decrease dzoomx or grossismx.'
83            STOP 1            STOP 1
# Line 85  contains Line 85  contains
85    
86         G = beta + (grossismx - beta) * fhyp         G = beta + (grossismx - beta) * fhyp
87    
88         ! Calcul de Xf         Xf(:nmax - 1) = beta * xtild(:nmax - 1) + (grossismx - beta) &
89                * ffdx(:nmax - 1)
        xxpr = beta + (grossismx - beta) * fxm(:nmax - 1)  
        Xf(0) = 0d0  
   
        DO i = 1, nmax - 1  
           Xf(i) = Xf(i-1) + xxpr(i) * (xtild(i) - xtild(i-1))  
        END DO  
   
90         Xf(nmax) = pi_d         Xf(nmax) = pi_d
91    
92         call invert_zoom_x(xf, xtild, G, rlonm025(:iim), xprimm025(:iim), &         call invert_zoom_x(xf, xtild, G, rlonm025(:iim), xprimm025(:iim), &
Legend:
Removed from v.146  
changed lines
  Added in v.147
  ViewVC Help
Powered by ViewVC 1.1.21