one can run the code as follows:

mpirun -np <nb_of_processors> ./douar

In this case, the program reads the 'input.txt' file and uses the user-defined boundary conditions in define_bc.f90.

On the other hand, on can also run the code as follows:

mpirun -np <nb_of_processors> ./douar input.txt_3Dp
mpirun -np <nb_of_processors> ./douar input.txt_2Dp
mpirun -np <nb_of_processors> ./douar input.txt_fold

these three input files are in the main directory, and if used, their corresponding boundary conditions are defined in 
define_bc.f90 and are automatically chosen by the code.

input.txt_3Dp  : square punch    (plasticity, one surface, linear viscosity)
input.txt_2Dp  : pseudo2D punch  (plasticity, one surface, linear viscosity)
input.txt_fold : folding problem (no plasticity, three surfaces, nonlinear viscosity)