322
Appendix B
B.6 Source Programs
C
C
Program: MCB2D.FOR
C*****************************************************************
C
C
C
Main program
C*****************************************************************
C
C
C
C
101
102
103
104
105
106
107
108
109
110
C
111
C
DIMENSION A(80000), N(10000), CHA(70)
IMPLICIT REAL*8 (A-H, O-Z)
MA=80000
MN=10000
OPEN(l,FILE='DATAN')
OPEN(2,FILE='DATAG')
OPEN(3,FILE='DATAI')
OPEN(4,FILE='DATAB')
OPEN(5,FILE='DATAP')
OPEN(6,FILE='RESUL')
OPEN(7,FILE='DATAW')
OPEN(8,FILE='DATAS')
OPEN(9,FILE='DATAO')
OPEN(10,FILE='DATAC')
OPEN(ll,FILE='DRAWC')
OPEN(12,FILE='DRAWT')
WRITE (* , 101)
WRITE (* , 102)
WRITE (* , 103)
WRITE (* , 104)
WRITE (* , 105)
WRITE (* , 106)
WRITE (* , 107)
WRITE (* , 108)
WRITE (* , 109)
WRITE (* , 110)
FORMAT (11111111)
FORMAT (15X, 1*************************************************')
FORMAT (15X, '*
*')
FORMAT (15X, '*
PROGRAM: MCB2D
version 2.0
*')
FORMAT (15X, '* For simulating coupled 2-D groundwater flow *')
FORMAT (15X, '* and mass transport in saturated zones
*')
FORMAT (15X, '*
*')
FORMAT (l5X, ' *
By Prof. Ne-Zheng Sun
*' )
format (15x, '*
*')
FORMAT (15X, 1*************************************************')
WRITE (*,111)
FORMAT (/,5X, 'Please wait ..... ',
&
1,5X, 'The program is reading control integers')
C Reading control integers
C
CALL JUMP(l, 2)
READ (1, 2) CHA