Appendix B - Computer Program Listing

? THIS PROGRAM SIMULATES METHANE FLOWRATES THROUGH THE GOAL ? SEAM BY THE GAUSS-SEIDEL ITERATIVE TECHNIQUE. THE MODEL IS A TWO ? DIMENSIONAL SYSTEM. THE ITERATIVE FORM OF THE GAUSS-SEIDEL IS ? OBTAINED FROM A PARTIAL DIFFERENTIAL EQUATION OF FLOW BY CENTRAL ? FINITE-DIFFERENCING APPROXIMATION SCHEME ? NPRESS IS THE NUMBER OF SPECIFIED PRESSURES EXCLUDING ? RESERVOIR AND FACE PRESSURES ? NQUAN IS THE NUMBER OF SPECIFIED FLOW QUANTITY.... EXCEPT FLOW ? OUT OF THE FACE. THIS- QUANTITY CAN BE POSITIVE(INJECTION) OR ? NEGATIVE ( PRODUCTION) ? XPERM(I,J) I S THE DIRECTIONAL PERMEABILITY IN THE X-DIRECTION AT THE POINT (I,J) ? YPERM(I,J) IS THE DIRECTIONAL PERMEABILITY IN THE Y-DIRECTION ? AT THE POINT (I,J) ? RESEP' IS THE RESERVOIR PRESSURE ? FACEP IS THE FACE PRESSURE. THIS PRESSURE IS ASSUMED TO BE ? APPROXIMATELY ATMOSPHERIC ? LENGTH, WIDTH, AND THICK ARE THE ZONE OR MEDIUM DIMENSIONS * PUHL I S A DUMMY NAME FOR VISCOSITY ? DMAXI IS VARIABLE NAME FOR INITIAL MAXIMUM ERROR OF ITERATION ? DMAXF IS A VARIABLE NAME FOR FINAL MAXIMUM ERROR OF ITERATION ? NITER IS A VARIABLE NAME FOR STORAGE OF NUMBER OF ITERATIONS ? RAD IS A VARIABLE NAME FOR THE SPECTRAL RADIUS OF THE GAUSS¬ ? SEIDEL MATRI-X. THIS PARAMETER ALSO MEASURES THE RATE OF ? CONVERGENCE OF THE ITERATION ? KOUNT IS AN ITERATION COUNTER ? MAXIT IS THE MAXIMUM NUMBER OF ITERATION SPECIFIED BY THE USER # EPS IS THE MAXIMUM RESIDUAL ALLOWABLE ERROR-FOR CONVERGENCE