The Fast Lagrangian Analysis of Continua (FLAC, Cundall, 1989) technique is a hybrid finite-element/finite-difference method that explicitly solves mass, momentum and heat conservation in a gridded domain. 


==============================
Compile
==============================

The folder /src contains the FORTRAN files that need to be compiled to produce an executable “acflac_tracer_sed_diff_track”.

$ cd src
$ bash make_acflac

This will generate a binary "acflac_tracer_sed_diff_track".

Supported compilers:
gfortran -- GNU (default), v4.2 or later, otherwise openmp is disabled
ifort    -- Intel
pgf90    -- PGI, the only compiler for openacc
f90      -- Sun Studio (less tested)


==============================
Run
==============================

The folder /input contains the input files for each model shown in Figure 4, 5, 6.
The corresponding parameter input file “drezina_**.inp” should be renamed as “drezina.inp”.

$ cd ../input/inputa
$ ../../src/acflac_tracer_sed_diff_track drezinp.inp

Running a simulation requires the corresponding parameter input file “drezina.inp” and "hydrother.inp" in the same file.

==============================
Map	
==============================

When a simulation is finished, copy all the file in Map file to the input file by running:

$ cp ../../Map/* ./;

Then run Map.m file in MATLAB, which will generate a series of VTU files (*.vtu). The VTU files can be visualized in Paraview to create the snapshots shown in Figure 4 and Figure 6.


==============================
==============================

If you have any questions, please send email to the first author. Ming Chen (chenming123@sjtu.edu.cn)