This repository contains the MATLAB code accompanying the paper "Constraints on magma source conditions from hydraulic fracture models and seismic observations of dyke tip deceleration".

# Constraints on magma source conditions from hydraulic fracture models and seismic observations of dyke tip deceleration

This repository contains the accompanying code.

## Repository Structure

### `ElasticIntegrals/`

Contains the elastic integrals for KGD (plane strain) and fixed-height fractures. 

### `ExtraFuncs/`

Folder containing extra functions used in the code

### `SelfSimilar/`

Solutions for a self-similar fracture profile - both KGD and PKN under different forcings. 

## Getting Started

### Main Simulation
To run the primary fixed-height simulation:
1. Add all files in to your MATLAB path
2. Execute: `ComputeWidthOverTimeNonDim_ILSA.m`