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 <em>newmodel_2a.py</em> 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<em> results</em> 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 <em>comparison </em>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).