SUBROUTINE SINT(TIME, T, TSTAR, II, S, SI, NS, OSINT)

C     + + + PURPOSE + + +
C
C     SUBROUTINE SINT CALCULATES THE INTEGRAL OF S WRT 
C     FROM 0.0 TO TIME
C
C     CALLED FROM HDEPTH, PHI, PSIS
C     AUTHOR(S): D. FLANAGAN, J. ASCOUGH
C     VERSION: THIS MODULE TAKEN FROM ASCOUGH STANDALONE IRS CODE
C     DATE CODED:  3-28-2005
C     CODED BY: D. FLANAGAN
C
C     + + + PARAMETER DECLARATIONS + + +
C      
      INTEGER MXTIME
      PARAMETER (MXTIME = 1000)

C     + + + ARGUMENT DECLARATIONS + + +
C
      INTEGER II, NS
      REAL TSTAR, T(MXTIME), S(MXTIME), SI(MXTIME+1), OSINT
      DOUBLE PRECISION TIME
     
C     + + + END SPECIFICATIONS + + +
C
      IF (TIME.LT.TSTAR) THEN
         IF (TIME.GE.T(II+1)) GO TO 20
         IF (TIME.LT.T(II)) GO TO 30
C        
   10    OSINT = SI(II) + S(II) * (TIME-T(II))
C
         RETURN
C        
   20    II = II + 1
C
         IF (TIME.GE.T(II+1)) GO TO 20
         GO TO 10
C
   30    II = II - 1
C
         IF (TIME.LT.T(II)) GO TO 30
         GO TO 10
      END IF
C
      OSINT = SI(NS+1)
C     
      RETURN
      END