C
RETURN
END
B.6. Source Programs
335
C********************************************************************
C
C
C
Subroutine DATAB reads boundary conditions
C********************************************************************
C
C
SUBROUTINE DATAB (HB1,CB1,NHBP,NCBP,STEPB,NHB,NCB,NTB,IHC)
IMPLICIT REAL*8 (A-H, O-Z)
DIMENSION HB1(NHB,NTB),CB1(NCB,NTB) ,NHBP(NHB),
1
NCBP(NCB),STEPB(NTB)
C Reading boundary conditions of the flow problem
C
IF (NHB .GT. 0) THEN
C
CALL JUMP (4, 2)
C
DO K=l, NTB
CALL JUMP (4, 3)
READ (4, 4) STEPB(K)
4 FORMAT (10X, Fl5.2)
CALL JUMP (4,1)
DO 1=1, NHB
READ (4, 10) NHBP(I) , HB1(I,K)
10 FORMAT (10x, IlO, F20.5)
END DO
END DO
ENDIF
C
C Reading boundary conditions of the transport problem
C
C
C
C
C
C
IF (IHC .EQ. 2) THEN
IF (NCB .GT. 0) THEN
CALL JUMP (4, 2)
DO K=l, NTB
CALL JUMP (4, 3)
READ (4, 4) STEPB(K)
CALL JUMP (4,1)
DO 1=1, NCB
READ (4, 20) NCBP(I), CB1(I,K)
20 FORMAT (10x, 110, F20.5)
END DO
END DO
ENDIF
ENDIF
RETURN
END
C***************************************************** ***************
C
C
C
Subroutine DATAP reads model parameters
C********************************************************************
C