326
Appendix B
C
C
C
C
CLOSE(8)
CLOSE(9)
CLOSE(10)
CLOSE(ll)
CLOSE(12)
STOP
END
SUBROUTINE JUMP (J,I)
DO K=l,I
READ (J, *)
END DO
RETURN
END
C*********************************************************************
C
C
Subroutine MCB2D solves coupled flow and mass transport problems
C
C********************************************************************
C
C
C
SUBROUTINE
1
2
3
4
5
6
7
8
1
2
MCB2D (X,Y,TR,BR,TH,PS,IJM,MA,MCPN,NPZON,NSZON,
HOO,COO,HB1,CB1,STEPB,NHBP,NCBP,
WRAT,WCON,WLX,WLY,LEW,STEPW,SRAT,SCON,STEPS,
DWR,DWC,DSR,DSC,
ZTK,ZSE,ZPR,ZAL,ZAT,ZAD,
DTK,DSE,DPR,DAL,DAT,DAD,PSE,PPR,PAD,
HO,H1,H2,CO,C1,C2,HFOO,HFO,CFOO,CFO,
DXX,DXY,DYY,VX,VY,CFA,CFB,IE,IEE,
HOB,COB,OBWX,OBWY,LEOB,OBT,
NP,NP1,NE,NB,NZON,IH,IHC,IST,ISO,IPR,
NW,NTW,NHB,NCB,NTB,NS,NTS,NOBW,NOBT)
IMPLICIT REAL*8 (A-H, O-Z)
DIMENSION X(NP) ,Y(NP),TR(NP) ,BR(NP) ,TH(NP) ,PS(NP),IJM(3,NE),
1
MA(NP1) ,MCPN(NB,NP) ,NPZON(NE) ,NSZON(NE) ,
2
HOO(NP),COO(NP) ,HB1(NHB,NTB),CB1(NCB,NTB) ,STEPB(NTB) ,
3
NHBP(NHB),NCBP(NCB),
4
WRAT(NW,NTW),WCON(NW,NTW),WLX(NW),WLY(NW),LEW(NW),
5
STEPW(NW),SRAT(NS,NTS),SCON(NS,NTS),STEPS(NTS),
6
DWR(NP) ,DWC(NP) ,DSR(NP),DSC(NP) ,
7
ZTK (NZON) , ZSE (NZON) , ZPR (NZON) , ZAL (NZON) , ZAT (NZON) ,
8
ZAD(NZON) ,DTK(NE) ,DSE(NE),DPR(NE) ,DAL(NE) ,DAT(NE) ,
1
DAD(NE) ,PSE(NP) ,PPR(NP) ,PAD(NP) ,
2
HO(NP),H1(NP),H2(NP),CO(NP),C1(NP),C2(NP),
3
HFOO(NP) ,HFO(NP) ,CFOO(NP) ,CFO(NP) ,CFA(IH) ,CFB(IH) ,
4
VX(NE) ,VY(NE) ,DXX(NE) ,DXY(NE) ,DYY(NE) ,IE(NP) ,IEE(NP) ,
5
HOB (NOBW,NOBT),COB(NOBW,NOBT),OBWX(NOBW),OBWY(NOBW),
6
LEOB(NOBW),OBT(NOBT),DTF(2),DTL(2),ALF(2),EPS(2)
C Reading initial and boundary conditions, sink/source terms,
C parameters, observations and computation control parameters.
C
C
WRITE (*, 101)
101 FORMAT (/,5X, 'The program is reading input data')
CALL DATAI (HOO,COO,NP,IHC)
CALL DATAB (HB1,CB1,NHBP,NCBP,STEPB,NHB,NCB,NTB,IHC)
CALL DATAP (ZTK,ZSE,ZPR,ZAL,ZAT,ZAD,DTK,DSE,DPR,DAL,DAT,