SUBROUTINE TDA1(NTIME,R,S,XD,DF,MODB,T)
      IMPLICIT DOUBLE PRECISION (A-Z)
      INTEGER IDIM
      PARAMETER (IDIM=50)
      DIMENSION R(IDIM),S(IDIM),XD(IDIM),A(4),T(9,IDIM)
      INTEGER NTIME,N1,I,J,MODB,NCELLS,DF
C
C         SUBROUTINE COMPUTES COMPONENT L1(I) FOR THE TEST BETWEEN
C         ModelS D AND A (WITH POOLING).
C
      NCELLS=0
      N1=NTIME-1
      DO 10 I=1,4
   10   A(I)=.0D0
      DO 40 I=1,NTIME-1
        A(1)=A(1)+R(I)
        A(2)=A(2)+S(I)-R(I)
        A(3)=A(3)+(.1D1-XD(I))*S(I)
        A(4)=A(4)+XD(I)*S(I)
        IF(A(1).LT..2D1.OR.A(2).LT..2D1.OR.
     .     A(3).LT..2D1.OR.A(4).LT..2D1)GO TO 40
        IF(I.NE.NTIME-2)GO TO 20
        IF(R(N1).LT..2D1.OR.(S(N1)-R(N1)).LT..2D1.OR.
     .    S(N1)*(.1D1-XD(N1)).LT..2D1.OR.S(N1)*XD(N1).LT..2D1)GO TO 40
   20   T(MODB,I)=(A(2)+A(1))*(A(1)-A(3))**2 /A(4)/A(3)
        IF(T(MODB,I).GE..0D0)GO TO 25
        T(MODB,I)=.0D0
        GO TO 40
   25   NCELLS=NCELLS+1
        T(MODB,IDIM)=T(MODB,IDIM)+T(MODB,I)
        DO 30 J=1,4
   30     A(J)=.0D0
   40   CONTINUE
      DF=DF-(N1-NCELLS)
      RETURN
      END