We extract here the most important information contained in this file :
Eulerian grid maximum number of nodes (x/y) :
parameter (mnx1=1601,mny1=121,mn1=mnx1*mny1) :
Lagrangian grid/cloud maximum number of nodes ( x/y) :
note that here an additional number of particles up to 20 times the grid size is
allowed for future injection in the Lagrangian cloud arrays
parameter (mn2inject=mn1*20,mn2=mnx2*mny2+mn2inject)
parameter (mne=2*mnx2*mny2)
Material colors :
NOTE THAT SEVERAL 'COLORS' MAY CORRESPOND TO THE SAME MATERIAL (PHYSICAL) SET
parameter (MAXCOLORS = 300)
parameter (MAX_MATSETS = 30) ! maximum number of materialsets
parameter (MAXITERS=1000)
Linear Solver (BLKFCT) :
THE FOLLOWING ARRAYS ARE NEEDED BY THE LINEAR SPARSE CHOLESKY SOLVER BLKFCT
 NMAX should be greater than the expected number of equations in the Stokes system (2*mnx1*mny1)
 ANZMAX should be bigger than the number of nonzeros in the matrix (18*mnx1*mny1)
 LNZMAX is the critical number and hard to guess exactly.
It should be bigger than the (a priori unknown)
number of nonzeros in the Cholesky factor L.
The value here is set to a number compatible with most runs
by users of SOPALE in the Dalhousie Geodynamics Group, on Dalhousie's IBM sp4 (p690.ucis.dal.ca)
Would you exceed this number, the program would stop and BLKFCT would return the expected number that
should be used in the file : blkfct.out
PARAMETER ( NMAX = 210000,
& ANZMAX = 3700000,
& SUBMAX = MAX(ANZMAX,100000),
& LNZMAX = 20000000,
& TMPMAX = 200000,
& IWMAX = 7*NMAX + 3 )
Remarks :
Some variables never used : e.g.
parameter (mnxt1=201,mnyt1=101)
