source: trunk/LMDZ.PLUTO.old/libf/phypluto/dtridgl.F @ 3232

      SUBROUTINE DTRIDGL(L,AF,BF,CF,DF,XK)

C  GCM2.0  Feb 2003

C     DOUBLE PRECISION VERSION OF TRIDGL

!      PARAMETER (NMAX=201)
      IMPLICIT REAL*8  (A-H,O-Z)
      DIMENSION AF(L),BF(L),CF(L),DF(L),XK(L)
!      DIMENSION AS(NMAX),DS(NMAX)
      DIMENSION AS(2*L),DS(2*L)

C*    THIS SUBROUTINE SOLVES A SYSTEM OF TRIDIAGIONAL MATRIX
C*    EQUATIONS. THE FORM OF THE EQUATIONS ARE:
C*    A(I)*X(I-1) + B(I)*X(I) + C(I)*X(I+1) = D(I)
C*    WHERE I=1,L  LESS THAN 103.
C* ..............REVIEWED -CP........

C======================================================================C

      AS(L) = AF(L)/BF(L)
      DS(L) = DF(L)/BF(L)

      DO I=2,L
        !print*, 'BF=',BF(L+1-I)
        !print*, 'AS=',AS(L+2-I)
        !print*, 'CF=',CF(L+1-I)
        !!! modif TB22:
        !BF(L+1-I)=max(BF(L+1-I),1e-20)
        X         = 1./(BF(L+1-I) - CF(L+1-I)*AS(L+2-I))
        AS(L+1-I) = AF(L+1-I)*X
        DS(L+1-I) = (DF(L+1-I)-CF(L+1-I)*DS(L+2-I))*X
      END DO
 
      XK(1)=DS(1)
      DO I=2,L
        XKB   = XK(I-1)
        XK(I) = DS(I)-AS(I)*XKB
      END DO
      !print*, 'TB22 dtrigl done'

      RETURN
      END