Sea Ice Dynamics in ICON-O - Numerical evaluation
Description
This data set includes the numerical data and the relevant source code accompanying the manuscript " Discretization of sea ice dynamics in the tangential plane on the sphere by a CD-grid-type finite element". In addition to the raw data, Python and shell scripts are deposited to process the data. The data set is used to evaluate the new discretization of the sea ice dynamics in the ocean model ICON-O. The approach is based on a CD-grid like approximation of the sea ice momentum equation, that is realized on an icosahedral grid. We numerically demonstrate the first order convergence in space by comparing the approximation to an analytic solution (4.1/ball.gp). We further show that a stabilization of the CD-grid on the icosahedral gird is required (4.1/vtk). To analyze the robustness of the CD-grid approach in ICON-O we solve the sea ice momentum equation on the sphere (4.2/ice_drift_sphere00002.vtk) and compare the results to a simulation conducted in two-dimensional planar geometry (4.2/ice_drift_planar00001.vtk). Finally, we compare the simulated sea ice drift from ICON-O with observations (4.3/Fig5.png). README.txt outlines how the routines are realized in ICON-O.
Files
Steps to reproduce
Exemplary shell commands -Example 4.1: ./tex2pdf.sh ball python3 ./3ddat2vtk.py ice_analytic_no_stabilization00002.dat -Example 4.2: python3 ./3ddat2vtk.py ice_drift_sphere00002.dat python3 ./dat2vtk.py ice_drift_planar00001.dat Example 4.3: see Mehlmann_and_Gutjahr_2022_Fig_5.pdf