TEXAS A AND M UNIV COLLEGE STATION DEPT OF OCEANOGRAPHY
This report describes a finite difference solution to a coupled pair of steady state non-linear differential equations that arise in a model of the dissolution of silicious tests on the sea floor Schink et al, 1975. An iterative procedure is used to treat the non-linear coupling between the equation describing solid phase concentration. An analytical solution for the case where solid phase concentration is constant with depth is used to test for the finite difference solution. The FORTRAN computer program used to perform the calculations is presented in an appendix.