C
C
B.6. Source Pro grams
339
READ (10, 10)
PT, DTF(2), DTL(2), ALF(2), EPS(2)
END IF
RETURN
END
c****************************************************·***************
C
C
Subroutine FDB forms boundary conditions for each time step
C
C********************************************************************
C
C
C
C
C
C
C
C
C
C
C
C
C
SUBROUTINE FDB (TSUM,HB1,CB1,HO,CO, STEPB,NHBP, NCBP,
NP,NHB,NCB,NTB,IHC)
IMPLICIT REAL*8 (A-H, O-Z)
DIMENSION HB1(NHB,NTB) ,CB1(NCB,NTB) ,HO(NP) ,CO(NP) ,
1
STEPB(NTB) ,NHBP(NHB) ,NCBP(NCB)
IF (NHB .GT. 0) THEN
DO K=l, NTB
IF (STEPB(K) .GT. TSUM) THEN
DO J=l, NHB
A=HB1(J, K-1)
T1=STEPB(K-1)
B=HB1 (J, K)
T2=STEPB(K)
CALL LIP(TSUM, A, B, Tl, T2, C)
HO (NHBP (J» =C
END DO
IF (IHC .EQ. 2
.AND. NCB .GT. 0) THEN
DO J=l, NCB
A=CB1(J, K-1)
T1=STEPB(K-1)
B=CB1 (J, K)
T2=STEPB(K)
CALL LIP(TSUM, A, B, Tl, T2, C)
CO (NCBP (J) ) =C
END DO
END IF
RETURN
END IF
END DO
END IF
RETURN
END
SUBROUTINE LIP
IMPLICIT REAL*8
(T, A, B, Tl,
(A-H, O-Z)
C=A+(B-A)*(T-T1)/(T2-T1)
RETURN
END
T2, C)