! ***********************************************************
! ********** Input file to a Diffpack simulator *************
! ***********************************************************

set tolerance = 1.0e-5   ! tolerance in outer (Gauss-Seidel) iteration loop
set max no of outer iterations = 10   ! for outer Gauss-Seidel loop

set gridfile = P=PreproBox|d=2 [0,1]x[0,1]|d=2 elm_tp=ElmTensorProd2 div=[4,4], grading=[1,1]

! physical parameters:
set power law exponent = { 1.0 & 0.5}  ! 4 test cases!!!
set alpha = {0.0 & 1.0}
set mu_0 = 1.0
set T_0 = 0.0
set compare to analytical solution = OFF
set initial guess = Boundary_layer  ! sinus_approx 1D_analytic

set redefine velocity boundary indicators = nb=3 names= fixed-wall moving-wall symmetry-plane 1=(1 2 3 4) 2=() 3=()
set redefine temperature boundary indicators = nb=4 names= x xx xxx xxxx 1=(1 2 3 4) 2=() 3=() 4=()

! **** momentum equation parameters: ****

set momentum pressure gradient = 1.0
sub momentum #1 NonLinEqSolver_prm   ! nonlinear solver in momentum eq.
set momentum #1 nonlinear iteration method = NewtonRaphson ! 
set momentum #1 max nonlinear iterations = 10
set momentum #1 max estimated nonlinear error = 1.0e-6
set momentum #1 nonlinear iteration stopping criterion = 1
ok
sub momentum #2 NonLinEqSolver_prm   ! nonlinear solver in momentum eq.
set momentum #2 nonlinear iteration method = SuccessiveSubst
set momentum #2 max nonlinear iterations = 10
set momentum #2 max estimated nonlinear error = 1.0e-6
set momentum #2 nonlinear iteration stopping criterion = 1
ok

! **** energy equation parameters: *****

sub energy #1 NonLinEqSolver_prm
set energy #1 nonlinear iteration method = NewtonRaphson 
set energy #1 max nonlinear iterations = 10
set energy #1 max estimated nonlinear error = 1.0e-6
set energy #1 nonlinear iteration stopping criterion = 1
ok
sub energy #2 NonLinEqSolver_prm
set energy #2 nonlinear iteration method = SuccessiveSubst
set energy #2 max nonlinear iterations = 10
set energy #2 max estimated nonlinear error = 1.0e-6
set energy #2 nonlinear iteration stopping criterion = 1
ok

! **** linear solvers (same system/solver used for both PDEs) ****

sub LinEqAdmFE
sub Matrix_prm                      ! submenu for matrix parameters
set matrix type = MatBand            ! matrix format, here: banded matrix
!set matrix type = MatSparse
ok
sub LinEqSolver_prm                 ! submenu for linear solvers
set basic method = GaussElim         ! Gaussian elimination is safe...
!set basic method = CGS
ok
sub Precond_prm                     ! submenu for preconditioning
set preconditioning type = PrecNone  ! no preconditioning, RILU: PrecRILU
set RILU relaxation parameter = 0.0  ! =0 gives ILU, =1 gives MILU etc
ok

sub ConvMonitorList_prm              ! termination criteria for linear solvers
sub Define ConvMonitor #1            ! criterion #1
set #1: convergence tolerance = 1.0e-4  ! set eps_r in  ||r|| < eps_r||b||
ok
ok
ok

sub SaveSimRes                       ! storage for visualization/plotting
set line1: start = (0.5,0)           ! plot unknowns along a line
set line1: stop = (0.5,1)            ! from (0.5,0) to (0.5,1)
set line1: resolution = 61           ! no of sampling points along the line
ok

set purpose = exercise for IMA Workshop on Software for PDEs

ok  ! quite menu and compute solution....