Data and code for "A comparison of the accuracy and stability of numerical schemes used to solve the equations derived from the shallow-shelf approximation"

This repository contains the scripts necessary to get and plot the data used in "A comparison of the accuracy and stability of numerical schemes used to solve the equations derived from the shallow-shelf approximation" The repository is organized into two folders, PlottingScripts and ScriptsToObtainTheData. All scripts used to plot the figures in "A comparison of the accuracy and stability of numerical schemes used to solve the equations derived from the shallow-shelf approximation" are found in the PlottingScripts folder with the name of the figure the script create round in the title. The relevant data to plot the scripts is also stored in that folder. The scripts used to obtain the data are found in the ScriptsToObtainTheData folder. The code to get the data used to plot figure 3 is found in getFig3Data. The code to run the finite element model is found in the FiniteElementData folder. In order to run the finite element model users must install Úa from https://ghilmarg.github.io/UaSource/. The path to Úa must then be added to runTransient.m, getBetaSingle.m and getGammaSingle.m. To perform a single transient run call runTransient.m with a choice of Theta, q, p and CFL number. The value of beta can be changed by altering the MeshN. The dispersion ratio can be obtained from getGammaSingle.m and the phase velocity from getBetaSingle.m and then getPhaseVelMutiple.m.