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

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

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

revision 134 by guez, Wed Apr 29 15:47:56 2015 UTC revision 144 by guez, Wed Jun 10 16:46:46 2015 UTC
# Line 10  contains Line 10  contains
10    
11      use coefpoly_m, only: coefpoly      use coefpoly_m, only: coefpoly
12      USE dimens_m, ONLY: iim      USE dimens_m, ONLY: iim
13        use dynetat0_m, only: clon
14      use nr_util, only: pi_d, twopi_d      use nr_util, only: pi_d, twopi_d
     use serre, only: clon  
15    
16      DOUBLE PRECISION, intent(in):: Xf(0:), xtild(0:), Xprimt(0:) ! (0:2 * nmax)      DOUBLE PRECISION, intent(in):: Xf(0:), xtild(0:), Xprimt(0:) ! (0:2 * nmax)
17      real, intent(out):: xlon(:), xprimm(:) ! (iim)      real, intent(out):: xlon(:), xprimm(:) ! (iim)
# Line 41  contains Line 41  contains
41         ! Calcul de Xf(xvrai(i))         ! Calcul de Xf(xvrai(i))
42    
43         xvrai(i) = xtild(it)         xvrai(i) = xtild(it)
44           IF (it == 2 * nmax) it = 2 * nmax -1
        IF (it == 2 * nmax) THEN  
           it = 2 * nmax -1  
        END IF  
   
45         CALL coefpoly(Xf(it), Xf(it + 1), Xprimt(it), Xprimt(it + 1), &         CALL coefpoly(Xf(it), Xf(it + 1), Xprimt(it), Xprimt(it + 1), &
46              xtild(it), xtild(it + 1), a0, a1, a2, a3)              xtild(it), xtild(it + 1), a0, a1, a2, a3)
47         Xf1 = Xf(it)         Xf1 = Xf(it)
Legend:
Removed from v.134  
changed lines
  Added in v.144
  ViewVC Help
Powered by ViewVC 1.1.21