Pre and post- processing files and scripts for the paper 'The effect of local shear stress on glacier sliding'

The folder contains:

• A set of .F90 to be compiled in order to run the simulations with Elmer Ice
• A python script called newmodel_2a.py for generating the complete analytical model shown in the appendix
• scripts for post-processing the results. The paths have to be edited. They read the .dat files, extract the main variables: u_b, \tau_b, N and compute C to produce the friction law plots. It also computes the ratio Af/As and produces the respective plots
• A compressed folder called results with the outputs of the main results. They include the results (.dat file) and the .sif files to run the simulations with FEM software Elmer/Ice. Their names are already included in a small subroutine in the postprocessing script
• A compressed folder called comparison with the outputs of the comparison between models. They include the results (.dat file) and the .sif files to run the simulations with FEM software Elmer/Ice, as well as the complete solution for the new models (in order to extract \tau_f for the sandpaper model we need the complete flow solution).