* Please note: first time you run PhreFlow, file 'inverse' will be created.
* This takes several hours, depending on your machine! Once this file is 
* created, do not delete it. File 'inverse' contains generic inverse matrix
* for all spheroidal inhomogeneities, regardless of their shape, size and 
* conductivity. Every time PhreFlow is executed, PhreFlow will search for 
* this file in the current directory. If the file exists, it will not be
* recreated. Hence, you need to create this file only once for all problems 
* you will ever solve with PhreFlow. 
*
* File 'solution' contains solution for a particular problem. This file is
* different for different problems. You should save this file if you want 
* to produce various grids, track various particles, etc based on the 
* same solution.
*
precision          4                              controls numbers of singularities (1-5, 5 is best)
spheroidprecision  5                              controls precision of spheroids (1-5, 5 is best)
relax              0.8                            relaxation coefficient for strengths of singularities
toler              1.d-7                          criterion for termination of iterations
maximumiterations  40                             for each time step
solvegridtrack                                    solve and grid and track                                    
budget                                            get budgets for every time step                                          
surface                                           compute the phreatic surface elevation
window             0. 0. 50. 100. 400 2           xcenter, ycenter, zcenter, size, number of nodes, coordinate axis to exclude (e.g. 3 means make x-y grid) 
gridstep 0                                        determines when to produce grids (see 'Window' command) and phreatic surface ('Surface' command)
2dpotential        0. 0. 50.                      These are coefficients of Taylor series used to represent 2D potential for lateral boundary conditions; a0, Real(a1), Imag(a1), Real(a2), Imag(a2), up to n=10 (2 through 10 are zero here).
domain             0. 0. 0. 50.                   base elevation, xcenter, ycenter, radius
conductivity       1.                             background conductivity
ImagedExtraction                                  well  
                   -35. 0. 15. 35. 0. 25.  20.    x1, y1, z1, x2, y2, z2,   extraction for stress period 1, for stress period 2, ...
&                                                 denotes end of element
ImagedExtraction                                  well  
                   0. -35. 19. 0. 40.  19. 20.    x1, y1, z1, x2, y2, z2,   extraction for stress period 1, for stress period 2, ...
&                                                 denotes end of element
recharge           0.0                            recharge for stress period 1, for stress period 2, ...
ImagedSpheroid                                    inhomogeneity
0. 0. 7.    30.  2.   0.1 0.1 1.     0.   2       x, y, z, long semiaxis, short semiaxis, x, y, z components of the rotational axis, conductivity, index (1=prolate, 2=oblate, 3 =sphere)
&                                                 denotes end of element
ImagedSpheroid                                    inhomogeneity
0. 0. 30.   30.  2.  -0.1 -0.1 1.    0.   2       x, y, z, long semiaxis, short semiaxis, x, y, z components of the rotational axis, conductivity, index (1=prolate, 2=oblate, 3 =sphere)
&                                                 denotes end of element
particle -50 -50 0 0  0 50 0. 1 1 10              xmin, xmax, ymin, ymax, zmin, zmax, release time, nx, ny, nz (specify nx*ny*nz particles)
particle  50  50 0 0  0 50 0. 1 1 10              xmin, xmax, ymin, ymax, zmin, zmax, release time, nx, ny, nz (specify nx*ny*nz particles)
time 300. 1000                                    total tracking time, printout resolution
end                                               of the file