search

ESCOMP / CISM

Community Ice Sheet Model
GNU Lesser General Public License v3.0
6 stars 11 forks source link

Updated the Dukowicz slab test case #37

Open whlipscomb opened 3 years ago

whlipscomb commented 3 years ago

This commit modifies the run and plot scripts for the Dukowicz slab test case, as described in Section 5 of Dukowicz (2012, The Cryosphere). The test case consists of an ice slab of uniform thickness moving down an inclined plane by a combination of sliding and shearing. Analytic Stokes and first-order velocity solutions exist for all values of Glen's n >= 1, and can be compared to those simulated by CISM. The original scripts contained warnings that the test was not fully supported. With these changes, the test is supported for CISM users.

This commit also adds a new script to run the solver stability tests described in Section 3.4 of Robinson et al. (2021, submitted to The Cryosphere), and updates the README file.

In most cases, answers are unchanged. To support the slab test, Glen's n is now a user-configurable variable. This might lead to roundoff-level changes on some machines.

whlipscomb commented 3 months ago

@Katetc, this PR can be removed, since the slab test changes were brought into main as part of the basal_physics4 branch.