# SOLVER CONFIGURATION
#--------------------- 
fem.io.partitioning: rank
# number of process that prints output (-1 no output) 
 fem.verboserank: 0
#verbosity for non-linear solver
 fem.solver.newton.verbose: false
# true for iteration output, false means no output
 fem.solver.verbose: false

fem.solver.newton.linabstol: 5e-7
fem.solver.newton.tolerance: 5e-7

#type of time discretisation (Backward Euler ==1)
phaseflow.theta: 1.

# parameters from time dependent problems 
fem.timeprovider.starttime: 0
phaseflow.endtime: 3600
phaseflow.timestep: 60e0
fem.timeprovider.factor: 1.
fem.timeprovider.updatestep: 1

fem.preconditioning: false

##################################


# DOMAIN CONFIGURATION #########
#---------------------

# macro grid file (3dgrid.dgf)
fem.io.macroGridFile_3d: ../data/3dgrid.dgf

#FIRST LENS BND  #######
phaseflow.coord_x_lens0: 0.5
phaseflow.coord_x_lens1: 0.75

phaseflow.coord_y_lens0: 0.25
phaseflow.coord_y_lens1: 0.75

phaseflow.coord_z_lens0: 0.6
phaseflow.coord_z_lens1: 0.8
##################################

#SECOND LENS BND  #######
phaseflow.coord_x_secondlens0: 0.25
phaseflow.coord_x_secondlens1: 0.75

phaseflow.coord_y_secondlens0: 0.25
phaseflow.coord_y_secondlens1: 0.75

phaseflow.coord_z_secondlens0: 0.2
phaseflow.coord_z_secondlens1: 0.4
##################################

#DIRICHLET BND
phaseflow.coord_x_Dirich_bnd0:0
phaseflow.coord_x_Dirich_bnd1:1

phaseflow.coord_y_Dirich_bnd0:0
phaseflow.coord_y_Dirich_bnd1:1

phaseflow.coord_z_Dirich_bnd0:0
phaseflow.coord_z_Dirich_bnd1:0
##################################


#NEUMMAN BND
phaseflow.coord_x_Neum_bnd0: 0.384375
phaseflow.coord_x_Neum_bnd1: 0.51625

phaseflow.coord_y_Neum_bnd0:0.384375
phaseflow.coord_y_Neum_bnd1: 0.51625

phaseflow.coord_z_Neum_bnd0:0.
phaseflow.coord_z_Neum_bnd1: 1
##################################



# GENERAL #######################
#--------

#### Parameters for output ######
phaseflow.computeEOC: 0
# prefix data files 
fem.io.datafileprefix: 3dtwophase
# save every i-th step
fem.io.savestep: 60.0e00

# specify directory for data output (is created if not exists)
fem.prefix: ../Output3D/Deg1/MonomialBasis

fem.dofmanager.memoryfactor: 1.1
# output format: binary | vtk-cell | vtk-vertex (projected in DG case) | none 
fem.io.outputformat: vtk-vertex
fem.io.path: ./
fem.io.savecount: 0

# upwinding of advective term
phaseflow.with_upw: 0

#slope limiter
phaseflow.with_limiter: 0

#local adaptivity
phaseflow.local_adapt: 1

# tolerance for estimator 
phaseflow.tolerance: 2.5

#number of global refinements before start of simulation 
phaseflow.level: 0

#number of level of refinement
phasefield.maxlevel:3

#DG penalty
phaseflow.penaltypress:0.5e-1
phaseflow.penaltysat: 0.5e-2

#DG method NIPG=-1 SIPG=1 IIPG=0
phaseflow.DGeps: 1

##################################


# PHYSICAL PARAMETER #############
#--------------

phaseflow.gravity_value: -9.81

# fluids properties
phaseflow.vicositynonwetting: 0.9e-3
phaseflow.vicositywetting: 1e-3
phaseflow.densitynonwetting: 1.460e3
phaseflow.densitywetting: 1e3

# Domain  
phaseflow.porosity: 0.4
phaseflow.entrypressure: 755
phaseflow.domain1absolutepermeability: 6.64e-11
phaseflow.inhomogeneitytheta: 2.7
phaseflow.residualnonwetting: 0
phaseflow.residualwetting: 0.12

# Lens  
phaseflow.lensabsolutepermeability: 6.64e-15
phaseflow.lensentrypressure: 5000
phaseflow.lensinhomogeneitytheta: 2.
phaseflow.lensporosity: 0.39
phaseflow.lensresidualwetting: 0.1

# SecondLens  
phaseflow.secondlensabsolutepermeability: 6.64e-16
phaseflow.secondlensentrypressure: 5000
phaseflow.secondlensinhomogeneitytheta: 2.
phaseflow.secondlensporosity: 0.39
phaseflow.secondlensresidualwetting: 0.1
##################################


 