The code used for Qin et al. (2021), Interactions between gas slug ascent and exchange flow in the conduit of persistently active volcanoes.

# JGR-Interactions between gas slug ascent and exchange flow in the conduit of 
persistently active volcanoes

Including the code used in Qin et al (2021): Interactions between gas
slug ascent and exchange flow in the conduit of persistently active volcanoes.

This repository includes the codes used to simulate the ascenting slug in 
exchange flow. 

Code Folder:

The code folder includes the codes used to reproduce the analogue experiment, 
which is developed using FORTRAN. Specifics on the development of the model and 
implementation are included in Qin et al (2021). To run this code, the reader 
should copy the whole folder into a single directory in his/her local computer 
or a server.

•    UMFPACK Library

Our codes use UMF library to solve the large sparse matrix. The UMFPACK library, 
written by Timothy Davis, is part of the SuiteSparse package, which is available
from http://www.cise.ufl.edu/research/sparse/ . The library files used by the 
authors have been included in the CODE/lib. To run a code, save all library 
files and the code in a same directory. However, these library files might not 
work in readers' local computer. Therefore, we suggest the readers to download 
and install their own library in their local computer. The UMFPACK use guide is 
in: https://fossies.org/linux/SuiteSparse/UMFPACK/Doc/UMFPACK_UserGuide.pdf.

•    Input

Our codes are suitable for a wide range of parameters. The input of geometry and
properties of fluids can be typed manually by readers in the “data input” area 
at the beginning of each code.

•    Output

We output the properties of multiphase, velocity and pressure profile frequently. 
Besides, in order to monitor the progress of the simulations, we output a runtime 
text file that prints residuals of the mechanical solvers and the position of 
the gas slug.

•    Run

To begin running a code, readers need to assign the name of the code in 
code/comp_2d2, then type and run “./comp_2d2” in a linux terminal to compile
the fortran code executable, which will be named "a.out". After the executable 
file was formed, type and run “./a.out”. The run time depends on the spatial and
temporal resolution. Simulation output files will be created in the same 
directory as the code is run.