This repository contains code for performing the experiments and processing the model output data for the publication "A Fully Implicit Second Order Method for Viscous Free Surface Stokes Flow - Application to Glacier Simulations" by Josefin Ahlkrona, André Löfgren and A. Clara J. Henry.The simulations were performed in Biceps as well as Elmer/Ice. Biceps is a prognostic two-dimensional full-Stokes ice-sheet solver written by André Löfgren. The latest version can be accesed at https://github.com/andrelofgrenSU/biceps, accessed: 2025-04-23, 2025. For the purpose of the simulations of the publication "A Fully Implicit Second Order Method for Viscous Free Surface Stokes Flow - Application to Glacier Simulations" by Josefin Ahlkrona, André Löfgren and A. Clara J. Henry, the files biceps_demo.cpp, p_stokes_fem.cpp, free_surface_fem.cpp and corresponding header files were edited. The edited version as used in this paper can be found at https://zenodo.org/records/17371013 To run a single simulaitons, use biceps_demo.cpp. To produce the convergence experiments of the paper, the bash-script convergence.sh can be used. convergence.sh depends on run_convergence.sh, run_reference.sh and plot_convergence, which can all be found in this repository and should be placed in the path biceps/demos/cppElmer/Ice is a finite element software under open source licensing. Elmer was first licensed under GPL (GNU General Public License, v. 2.1) license in 2008. Since 2012 the ElmerSolver library has been published under LGPL license (GNU Lesser General Public License, v. 2.0). The Elmer revision used in this study can be found at https://zenodo.org/records/7892181. The extra files needed to perform the simulations in "A Fully Implicit Second Order Method for Viscous Free Surface Stokes Flow - Application to Glacier Simulations" by Josefin Ahlkrona, André Löfgren and A. Clara J. Henry are included in this repository, namely IncompressibleNSVecMod.F90 which should be placed in the path elmerfem/tree/devel/fem/src/modules The finite element mesh was generated with gmsh based on the file basemsh.geo which can be found in this repository. And the experiments can be performed by running Forward.sif which depends on the .dat files included in this repository. All output from the simulations can be found in the zip-file.