# SAFIM: Semi-Analytical Floating Iceberg Model
# v1.0.0 DOI: 10.5281/zenodo.3562012
# v1.0.1 DOI: 10.5281/zenodo.3583022 (improved indentation)
# v1.0.2 DOI: 10.5281/zenodo.3783774 (corrected bug in force storage for Stormer-Verlet integration in time).

This is a Python implementation of Semi-Analytical Floating Iceberg Model.
It models capsize of unstable icebergs in free ocean (without constraints),
the model includes gravity, hydrostatic water pressure, added mass effect, and pressure drag.
The integration in time is made by Störmer-Verlet method.

Three files are provided:
(1) "safim.py" - implementation of the capsizing iceberg interacting with open water.
(2) "Example_of_simulation.ipynb" and (3) "Example_of_simulation.py" -> Python notebook and identical python executable (for Linux) to run simulations.
To run under Windows OS, simply remove the first line "#!/usr/bin/env python".

Example_of_simulation produces two figures:
"Force_x.pdf" - evolution of the horizontal force in time and
"Coord_z.pdf" - evolution of the vertical position
and stores the simulation results in
"Results.npz"
These results are obtained for iceberg of height 800 m, aspect ratio 0.25 initially tilted by 1 degree and being in hydrostatically equilibrated vertically-wise position.
All the parameters are listed in "Example_of_simulation.ipynb" or "Example_of_simulation.py".