!
! The input file for running the PARIS code
! Parameters are read using a namelist statement.
! Blank lines and commented lines are ignored.
! Location of the parameters in the list is not important.

&PARAMETERS
                   
! name of the namelist
!=================================================================================================
! General parameters

TwoPhase  = T        ! TwoPhase: Is the flow two-phase?
DoVOF 	  = T
DoFront   = F
GetPropertiesFromFront = F ! T: uses Front-Tracking data to compute mu, rho and Surface Tension. 
		           ! F: uses VOF data to compute mu and rho. 

ZeroReynolds = F
Implicit = IMPTEMP     
hypre    = F       ! T: uses hyprepackage, F: uses SOR solver 
restart  = F       ! T: start the domain from a previous simulation
restartFront = F   ! T: start the frontfrom  a previous simulation
restartAverages = F
nBackup = 10000    ! number of time steps between  backups are kept.

NSTEP   = 50000  ! maximum number of time steps
EndTime = 1.0      ! When to stop simulation

MaxDt   = 0.1      ! maximum size of time step
dtFlag  = 1        ! 1: fixed dt;  2: fixed CFL
dt      = DTTEMP   ! dt in case of dtFlag=1
CFL     = 0.9
MAXERROR = MAXERRTEMP    ! Residual for Poisson solver
MAXERRORVOL = 1d-4
                   ! Numerical parameters

ITIME_SCHEME = 1   ! time scheme: 1:first order Euler, 2: second order quasi Crank-Nicolson
                   ! Should use first order Euler for Stokes flow ! ! ! 
MAXIT    = 1000
BETA     = 1.3     ! parameters for linear solver

!=================================================================================================
! Output parameters

termout  = 20
ICOut    = T        ! output initial condition
NOUT     = 200      ! write the solution to file every nout time steps
output_format = 2   ! 1:tecplot 2:vtk
out_path = 'out'    ! name of the output directory

nstats   = 10      ! number of time steps between flow statistics calculations

!=================================================================================================
!Grid parameters
npx     = NPXTEMP
npy     = NPYTEMP
npz     = NPZTEMP
                   ! number of processors in x,y,z direction
                   ! Total number of processes should be at least npx*npy*npz+1 since one process
                   ! is dedicated to the front. Any extra processes will be idle.
Nx       = NXTEMP
Ny       = NXTEMP
Nz       = NXTEMP
Ng       = 2
                   ! total grid size in x,y,z direction and number of ghost cells

read_x   = F       ! read the grid file for x-grid; If true xLength and xform are neglected
read_y   = F
read_z   = F

x_file   = 'xh.dat'	! input file for xh (Nx+1 points)
y_file   = 'yh.dat'
z_file   = 'zh.dat'

xform    = 0.0		!1.0
yform    = 0.0		!1.0
zform    = 0.0       	! non-uniformity of the grid if not reading an input file
                     	! 0:uniform; +:clustered at center; -:clustered near ends

XLENGTH  = 1.0
YLENGTH  = 1.0
ZLENGTH  = 1.0
                   ! domain size in x,y,z direction

!=================================================================================================
! Flow parameters

WallVel(1,1)  = 1. ! can set inflow velocities here
WallShear(4,1) = 0.0

U_init =  1d0
GX      = 0.0
GY      = 0.0
GZ      = 0.0
                   ! Components of the gravity in x,y,z direction

dPdX      = 0.0
dPdY      = 0.0
dPdZ      = 0.0
                    ! Px, Py, Pz: pressure gradients in case of pressure driven channel flow
                    ! Px = (P(xLength)-P(0))/xLength
BDRY_COND = 3  1  1  4  1  1
                    !Type of boundary condition in x,y,z direction: 0:wall  1:periodic  2:shear
		    !  3: inflow  4:outflow
                    !x- y- z- x+ y+ z+

RHO1     = 0.03591954022 
MU1      = 0.10321706962e-5
                    ! rho1, mu1 : density and viscosity of the matrix phase

RHO2     = 1.0
MU2      = 2.87356321839e-5
SIGMA    = 3.0e-5   ! rho2, mu2, sigma : properties of the drop and surface tension
                   

NumBubble = 0       ! number of bubbles

MaxPoint = 2000000
MaxElem  = 4000000
MaxFront = 100
amin     = 0.32
amax     = 0.96
aspmax   = 1.54

smooth   = T        !smooth the interface
nsmooth  = 10       !every nsmooth time steps
nregrid  = 10       !regrid the front every nregrid time steps

BUOYANCYCASE = 0
                   ! BuoyancyCase : determines what density will be subtracted from the gravity 
                   ! body force.
                   ! 0: rro=0,  1: rro=rho1,  2: rro=rho2,  3: rro=average(rho)
!=================================================================================================

/
! end of the namelist