-------------------------------------------------------------------------- 
                   Lithosphere and Mantle Evolution Model                   
     Compiled: Date: Jan 28 2021 - Time: 17:35:55 	    
-------------------------------------------------------------------------- 
        STAGGERED-GRID FINITE DIFFERENCE CANONICAL IMPLEMENTATION           
-------------------------------------------------------------------------- 
Parsing input file : ./t20_FSSA/RTI_FSSA.dat 
   Adding PETSc option: -snes_type ksponly
   Adding PETSc option: -objects_dump
Finished parsing input file : ./t20_FSSA/RTI_FSSA.dat 
--------------------------------------------------------------------------
Scaling parameters:
   Temperature : 1300. [C/K] 
   Length      : 100000. [m] 
   Viscosity   : 1e+18 [Pa*s] 
   Stress      : 1e+08 [Pa] 
--------------------------------------------------------------------------
Time stepping parameters:
   Simulation end time          : 5.5 [Myr] 
   Maximum number of steps      : 20 
   Time step                    : 0.01 [Myr] 
   Minimum time step            : 1e-05 [Myr] 
   Maximum time step            : 1. [Myr] 
   Time step increase factor    : 0.1 
   CFL criterion                : 0.75 
   CFLMAX (fixed time steps)    : 0.75 
   Output time step             : 0.25 [Myr] 
   Output every [n] steps       : 100 
   Output [n] initial steps     : 1 
--------------------------------------------------------------------------
--------------------------------------------------------------------------
Material parameters: 
--------------------------------------------------------------------------
   Phase ID : 0     --   Lithosphere 
   (dens)   : rho = 3300. [kg/m^3]  
   (diff)   : eta = 1e+21 [Pa*s]  Bd = 5e-22 [1/Pa/s]  

   Phase ID : 1     --   Asthenosphere 
   (dens)   : rho = 3200. [kg/m^3]  
   (diff)   : eta = 1e+20 [Pa*s]  Bd = 5e-21 [1/Pa/s]  

   Phase ID : 2     --   StickyAir 
   (diff)   : eta = 1e+16 [Pa*s]  Bd = 5e-17 [1/Pa/s]  

--------------------------------------------------------------------------
--------------------------------------------------------------------------
Grid parameters:
   Total number of cpu                  : 1 
   Processor grid  [nx, ny, nz]         : [1, 1, 1]
   Fine grid cells [nx, ny, nz]         : [50, 1, 100]
   Number of cells                      :  5000
   Number of faces                      :  20150
   Maximum cell aspect ratio            :  5.00000
   Lower coordinate bounds [bx, by, bz] : [-250., -1., -500.]
   Upper coordinate bounds [ex, ey, ez] : [250., 1., 100.]
--------------------------------------------------------------------------
Free surface parameters: 
   Sticky air phase ID       : 2 
   Initial surface level     : 0. [km] 
   Erosion model             : none
   Sedimentation model       : none
   Correct marker phases     @ 
--------------------------------------------------------------------------
Boundary condition parameters: 
   No-slip boundary mask [lt rt ft bk bm tp]  : 0 0 0 0 1 0 
--------------------------------------------------------------------------
Solution parameters & controls:
   Gravity [gx, gy, gz]                    : [0., 0., -10.] [m/s^2] 
   Surface stabilization (FSSA)            :  1. 
   Use lithostatic pressure for creep      @ 
   Minimum viscosity                       : 1e+16 [Pa*s] 
   Maximum viscosity                       : 1e+23 [Pa*s] 
   Max. melt fraction (viscosity, density) : 0.15    
   Rheology iteration number               : 25  
   Rheology iteration tolerance            : 1e-06    
   Ground water level type                 : none 
--------------------------------------------------------------------------
Advection parameters:
   Advection scheme              : Euler 1-st order (basic implementation)
   Marker setup scheme           : geometric primitives
   Velocity interpolation scheme : STAG (linear)
   Marker control type           : subgrid 
   Markers per cell [nx, ny, nz] : [5, 5, 5] 
   Marker distribution type      : uniform
   Background phase ID           : 0 
--------------------------------------------------------------------------
Reading geometric primitives ... done (0.0198972 sec)
--------------------------------------------------------------------------
Output parameters:
   Output file name                        : RT_FSSA 
   Write .pvd file                         : yes 
   Phase                                   @ 
   Density                                 @ 
   Total effective viscosity               @ 
   Creep effective viscosity               @ 
   Velocity                                @ 
   Pressure                                @ 
--------------------------------------------------------------------------
Surface output parameters:
   Write .pvd file : yes 
   Velocity        @ 
   Topography      @ 
   Amplitude       @ 
--------------------------------------------------------------------------
AVD output parameters:
   Write .pvd file       : yes 
   AVD refinement factor : 3 
--------------------------------------------------------------------------
Preconditioner parameters: 
   Matrix type                   : monolithic
   Preconditioner type           : user-defined
--------------------------------------------------------------------------
Solver parameters specified: 
   Outermost Krylov solver       : gmres 
   Solver type                   : serial direct/lu 
   Solver package                : petsc 
--------------------------------------------------------------------------
Saving output ... [mpiexec@Boris-Mac] Sending Ctrl-C to processes as requested
[mpiexec@Boris-Mac] Press Ctrl-C again to force abort
