Sea Ice Dynamics in ICON-O - Numerical evaluation
Description
This data set includes the numerical data and the relevant source code accompanying the manuscript " Modeling sea ice dynamics in the tangential plane on the sphere". 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 demonstrate that the CD-grid discretization of the sea ice dynamics in ICON-O is less diffusive than the A-grid discretization(4.3/extent_ICON.gp) and allows a larger flow through narrow channels (4.3/transport_ICON.gp). We visualize the ice thickness and concentration of the A-grid (4.3/Agrid_ice_dynamics01392.vtk) and CD-grid type approach (4.3/Agrid_ice_dynamics01392.vtk) after two month. 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: ./tex2pdf.sh transport_ICON python3 ./3ddat2vtk.py CDgrid_ice_dynamics01392.dat