B.1.7 I6LMDA For the a matrix computed in 15CMAT, I6LMDA computes the eigenvalues using Scientific Subroutine Package subprograms HSBG and ATEI6. The largest eigenvalue in magnitude is then found. B.1.8 I7MPSP The subprogram I7MPSP proposes modifications to solution procedures by re-ordering the functions and variables in the FVIM and specifying that certain of them belong to Newton-Raphson solution procedures and others to Gauss-Seidel solution procedures. The functions and variables to be included in a Newton-Raphson partition are determined from an examination of the modified Jacobian (coefficient matrix). Non-output elements which are much larger than the output terms in their row and column are found. The function type containing these terms is moved to the Newton-Raphson partition, along with the variable type if it is a legal output; otherwise another variable type which is a legal output is moved. Only elements which are a certain percent of the largest element (user specified) are treated in this manner. B.2 Level 2 Subprograms B.2.1 J1FIX STOSIM (3.2), J1CTWE (5.2), RSTSIM (5.11), WRTDIR (6.2) J1FIX is invoked by ALWET when a variable proposed as a torn variable is found to be a two-way edge complement of a variable which was torn because of that two-way edge. The SIM, as it existed when the tear decision for the two-way edge was made, is restored using RSTSIM. The two-way edge links are traced using the subprogram