# RECONNECTION SIMULATION WITH PARTICLES TRACKING
# SIMULATION USED IN IVY’S PAPER

# DIRECTORIES WITHOUT  / at the end !!!!!
SaveDirName = data
RestartDirName = data

# set the maximum number of particles allocated  
NpMaxNpRatio = 2.0

# New flags:
Case              = GEM       # Select the case
PoissonCorrection = no       # Poisson correction
WriteMethod       = shdf5     # Output method [ phdf5, shdf5, pvtk]
SimulationName    = MagnetotailEngParticle  # Simulation name for the output

#  %%%%%%%%%%%%%%%%%%% Magnetic Reconnection %%%%%%%%%%%%%%%%%%
B0x = 0.01559
B0y = 0.0
B0z = 0.0
delta = 0.5
    
#  %%%%%%%%%%%%%%%%%%% TIME %%%%%%%%%%%%%%%%%%
dt = 0.125          # dt = time step   
ncycles = 18001     # cycles
th = 1.0           # th =   decentering parameter
c = 1.0

#  %%%%%%%%%%%%%%%%%%% SMOOTH %%%%%%%%%%%%%%%%%%
Smooth = 0.2      # Smoothing value (5-points stencil)

# %%%%%%%%%%%%%%%%%% BOX SIZE %%%%%%%%%%%%%%%
Lx = 120   # Lx = simulation box length - x direction    
Ly = 40   # Ly = simulation box length - y direction   
Lz = 10   # Lz = simulation box length - z direction
nxc = 3072  # nxc = number of cells - x direction        
nyc = 1024  # nyc = number of cells - y direction 
nzc = 1  # nzc = number of cells - z direction 

# %%%%%%%%%%%%%% MPI TOPOLOGY %%%%%%%%%%%%%%
# number of MPI subdomains in each direction
XLEN = 48
YLEN = 16
ZLEN = 1
# topology of subdomains in each dimension (1=true, 0=false)
PERIODICX = 1
PERIODICY = 0
PERIODICZ = 1

# %%%%%%%%%%%%%% PARTICLES %%%%%%%%%%%%%%%%%
#    ns = number of species
#    0 = electrons
#    1 = protons
#    ...
ns = 4

# qom = charge to mass ratio for different species */
qom = -100.0 1.0 -100.0 1.0

# Initial density (make sure that plasma is neutral)
rhoINIT = 1.0 1.0 0.1 0.1
rhoINJECT=0.0 0.0 0.0 0.0
# TrackParticleID[species] = 1=true, 0=false --> Assign ID to particles 
TrackParticleID= 1 1 1 1 
# npcelx = number of particles per cell - Direction X 
npcelx = 8 8 8 8 
# npcely = number of particles per cell - Direction Y */
npcely = 8 8 8 8 
# npcelz = number of particles per cell - Direction Z */
npcelz = 8 8 8 8


# uth = thermal velocity for different species - Direction X  */
uth  = 0.045 0.0100 0.045 0.0100
# vth = thermal velocity for different species - Direction Y */
vth  = 0.045 0.0100 0.045 0.0100
# wth = thermal velocity for different species - Direction Z */
wth  = 0.045 0.0100 0.045 0.0100
# u0 = drift velocity   - Direction X    */
u0 = 0.0 0.0 0.0 0.0
# v0 = drift velocity   - Direction Y    */
v0 = 0.0 0.0 0.0 0.0
# w0 = drift velocity   - Direction Z    */
w0 = 0.005196 -0.025980 0.0 0.0 

# &&&&&&&&&&&& boundary conditions &&&&&&&&&&&&&&&
# PHI Electrostatic Potential     
# 0,1 = Dirichilet boundary condition      
# 2   = Neumann boundary condition
bcPHIfaceXright = 1
bcPHIfaceXleft  = 1
bcPHIfaceYright = 1
bcPHIfaceYleft  = 1
bcPHIfaceZright = 1
bcPHIfaceZleft  = 1
    
#    EM field boundary condition
#    0 = perfect conductor
#    1 = magnetic mirror

    bcEMfaceXright = 0
    bcEMfaceXleft =  0
    bcEMfaceYright = 0
    bcEMfaceYleft =  0
    bcEMfaceZright = 0
    bcEMfaceZleft =  0
#    Particles Boundary condition
#    0 = exit
#    1 = perfect mirror
#    2 = riemission
#    Caveat: if your processor topology is set to be periodic in a direction automatically the boundary condition in that direction will be periodic*/
    bcPfaceXright = 1
    bcPfaceXleft =  1
    bcPfaceYright = 1
    bcPfaceYleft =  1
    bcPfaceZright = 1
    bcPfaceZleft =  1

#  print to video results */
    verbose = 1

# velocity of the injection from the wall 
    Vinj= 0.0
	
# CG solver stopping criterium tolerance 
   CGtol = 1E-3
# GMRES solver stopping criterium tolerance 
   GMREStol = 1E-3
# mover predictor corrector iteration 
   NiterMover  = 3
# Output for field 
   FieldOutputCycle= 100
   FieldOutputTag= Eall+Ball+Jsall+phi
   MomentsOutputTag =rhos+pressure
# Output for particles if 1 it doesnt save particles data
   ParticlesOutputCycle = 100
   ParticlesOutputTag= ID + position + velocity  
# restart cycle 
   RestartOutputCycle = 0
# record the last cyle, true by default
   CallFinalize = 0