This is the code used for the paper "Dynamics of the global meridional ice flow of Europa's icy shell", Yosef Ashkenazy, Roiy Sayag and Eli Tziperman, Nature Astronomy, 2017.

Different experiments shown in the paper:

1) no ocean, no convection

2) with ocean, no convection

3) no ocean, with convection

4) with ocean, with convection

correspond to different flags used in the code, using the variables
p.flag_ocean and p.simple_convection

The supplementary information to the above paper contains model experiments with salinity, activated using p.flag_salt

The file strain-rate-Nimmo-et-al-2007.txt is a required input file for this code, and the output/ directory contains an example output from a short test run of the model.