4. Test simulations

These are a number of computationally rather inexpensive test simulations, of which the run-specs header files are contained in the SICOPOLIS repository.

Run v5_vialov3d25
3-d version of the 2-d Vialov profile (Vialov [31]),
SIA, resolution 25 km, \(t=0\ldots{}100\,\mathrm{ka}\).
Similar to the EISMINT Phase 1 fixed-margin experiment (Huybrechts et al. [23]), but without thermodynamics. Instead, isothermal conditions with \(T=-10^{\circ}\mathrm{C}\) everywhere are assumed.
Run v5_emtp2sge25_expA
EISMINT Phase 2 Simplified Geometry Experiment A,
SIA, resolution 25 km, \(t=0\ldots{}200\,\mathrm{ka}\) (Payne et al. [29]).
The thermodynamics solver for this run is the one-layer melting-CTS enthalpy scheme (ENTM), while all other runs employ the polythermal two-layer scheme (POLY) (Greve and Blatter [18]).
Run v5_grl16_bm5_ss25ka
Greenland ice sheet, SIA, resolution 16 km,
short steady-state run (\(t=0\ldots{}25\,\mathrm{ka}\)) for modern climate conditions (unpublished).
Run v5_ant40_b2_ss25ka
Antarctic ice sheet without ice shelves, SIA, resolution 40 km,
short steady-state run (\(t=0\ldots{}25\,\mathrm{ka}\)) for modern climate conditions (unpublished).
Run v5_grl20_b2_paleo21
Greenland ice sheet, SIA, resolution 20 km,
\(t=-140\,\mathrm{ka}\ldots{}0\), basal sliding ramped up during the first 5 ka.
Modified, low-resolution version of the spin-up for ISMIP6 InitMIP (Greve et al. [19]).
Runs v5_grl10_b2_{paleo21, future21_ctrl, future21_asmb}
Greenland ice sheet, SIA, resolution 10 km,
\(t=-9\,\mathrm{ka}\ldots{}0\) for the paleo run, \(t=0\ldots{}100\,\mathrm{a}\) for the future runs.
10-km version of the spin-up and the schematic future climate runs for ISMIP6 InitMIP (Greve et al. [19]).
Runs v5_ant64_b2_{spinup09_init100a, spinup09_fixtopo, spinup09, future09_ctrl}
Antarctic ice sheet with hybrid shallow-ice–shelfy-stream dynamics
(Bernales et al. [3]) and ice shelves (SSA), resolution 64 km;
\(t=-140.1\ldots{}-140.0\,\mathrm{ka}\) for the init run without basal sliding (…_init100a),
\(t=-140\,\mathrm{ka}\ldots{}0\) for the run with almost fixed topography (…_fixtopo), basal sliding ramped up during the first 5 ka,
\(t=-0.5\,\mathrm{ka}\ldots{}0\) for the final, freely-evolving-topography part of the (…_spinup09),
\(t=0\ldots{}100\,\mathrm{a}\) for the constant-climate control run (…_future09_ctrl).
64-km version of the spin-up and the constant-climate control run for ISMIP6 InitMIP (Greve and Galton-Fenzi, pers. comm. 2017).
Runs v5_asf2_steady, v5_asf2_surge
Austfonna, SIA, resolution 2 km, \(t=0\ldots{}10\,\mathrm{ka}\).
Similar to Dunse et al. [8]’s Exp. 2 (steady fast flow) and Exp. 5 (surging-type flow), respectively.
Runs v5_nmars10_steady, v5_smars10_steady
North-/south-polar cap of Mars, SIA, resolution 10 km, \(t=-10\,\mathrm{Ma}\ldots{}0\).
Steady-state runs by Greve [16].
Run v5_nhem80_nt012_new
Northern hemisphere, SIA, resolution 80 km, \(t=-250\,\mathrm{ka}\ldots{}0\).
Similar to run nt012 by Greve et al. [22].
Run v5_heino50_st
ISMIP HEINO standard run ST, SIA, resolution 50 km, \(t=0\ldots{}200\,\mathrm{ka}\) (Calov et al. [7]).

Computing times:

Run

Model time

Time step

CPU time

v5_vialov3d25

100 ka

20 a

1.0 min

v5_emtp2sge25_expA

200 ka

20 a

3.9 min

v5_grl16_bm5_ss25ka

25 ka

5 a

9.7 min

v5_ant40_b2_ss25ka

25 ka

10 a

5.0 min

v5_grl20_b2_paleo21

140 ka

5 a

0.8 hrs

v5_grl10_b2_paleo21*

9 ka

1 a

1.0 hrs

v5_grl10_b2_future21_ctrl

100 a

1 a

0.9 min

v5_grl10_b2_future21_asmb

100 a

1 a

0.9 min

v5_ant64_b2_spinup09_init100a

100 a

2 / 10 a

4.1 sec

v5_ant64_b2_spinup09_fixtopo

140 ka

5 / 10 a

0.7 hrs

v5_ant64_b2_spinup09

500 a

2 / 10 a

0.5 min

v5_ant64_b2_future09_ctrl

100 a

2 / 10 a

6.1 sec

Table 1: Model times, time steps and computing (CPU) times for the EISMINT, Greenland and Antarctica test simulations contained in the script multi_sico_1.sh, run with SICOPOLIS V5-dev (branch develop, revision 9c909c3c2) and the Intel Fortran Compiler 19.1 for Linux (optimization options -xHOST -O3 -no-prec-div) on a 12-Core Intel Xeon Gold 6256 (3.6 GHz) PC under openSUSE Leap 15.4.
: If one value is given, this is the common dynamic (velocity, ice thickness) and thermodynamic (temperature, water content, age) time step. If two values are given (marked by the dagger () symbol), the first one is the dynamic, the second one the thermodynamic time step.
: All runs were done on a single core only. The v5_ant64_b2_xxx runs that include ice shelves can be done on multiple cores using OpenMP for the SSA solver. However, at the employed, low resolution of 64 km the solver does not scale well, and the gain in wall clock time by using multiple cores is very small.
*: For this run, see the remark in the subsection on the resolution-doubler tool.