Sia fluxlim is a numerical method enhancement module for shallow ice approximate flow code, which implements the MUSCL super bee flux limiter algorithm. Historical background : This module was developed by open-source contributors to enhance the numerical performance of shallow ice approximation flow code, particularly in terms of stability and accuracy when dealing with shock waves and discontinuous solutions. The development of this module responds to the demand for more stable and accurate numerical methods in ice flow simulation. Technical features : • Implemented MUSCL super bee flux limiter algorithm to improve numerical stability • Enhanced ability to handle shock waves and discontinuous solutions • Improved the accuracy of numerical simulation of ice sheet flow • Stable performance in challenging ice flow scenarios • Easy integration with existing shallow ice approximation codes • High computational efficiency without significantly increasing computational costs Core functions : • Provide numerical method enhancement for shallow ice approximate flow codes • Improve the numerical stability of ice flow simulation • Dealing with shock waves and discontinuous solutions in ice flow • Improve the accuracy of numerical simulation of ice sheet flow • Provide robust performance in challenging ice flow scenarios • Seamless integration with existing code Application case : • Flow simulation of the edge region of the Greenland ice sheet • Simulation of rapidly changing glacier flows • Simulation of regions with significant gradients in ice flow • Improvement of numerical method for shallow ice approximate flow code • Accuracy evaluation of numerical simulation of ice sheet flow • Comparative study of different flux limiter algorithms Limitations : • Only applicable to shallow ice approximate flow codes, not suitable for complete Stokes equation solvers • The handling of certain extreme flow scenarios still needs improvement • Integration with certain specific code may require additional modifications • The selection of parameters for the flux limiter has a certain impact on the simulation results • Unable to solve the physical limitations of shallow ice approximation itself Input parameters : • Geometric data of ice sheet (thickness, surface elevation, etc.) • Physical parameters of ice (rheological parameters, etc.) • Numerical method parameters (flux limiter coefficients, etc.) • Boundary conditions and initial conditions • Simulate time steps and duration Output result : • Ice sheet velocity field and thickness variation • Numerical stability and convergence indicators • Comparison of simulation results under different flux limiter parameters • The processing effect of shock waves and discontinuous solutions