dome

The time evolution of a 3-D paraboloid dome of ice sitting on a flat bed. This is a very basic ice-sheet like test case that is simple and relatively fast to run. It confirms that the basic higher-order model physics are working, but it does not strenuously test the physics and boundary conditions.

p1

Bit for Bit

Variable Max Error Index of Max Error RMS Error Plot
velnorm 0.00000e+00 N/A 0.00000e+00
thk 0.00000e+00 N/A 0.00000e+00

Configuration Comparison

---

+++

@@ -23,7 +23,7 @@

which_ho_babc = 4 # 4 = no-slip at bed

which_ho_efvs = 2 # 2 = nonlinear eff. visc. w/ n=3

which_ho_sparse = 3 # 1 = SLAP GMRES, 3 = glissade parallel PCG, 4 = Trilinos for linear solver

-which_ho_nonlinear = 0 # 0 = Picard, 1 = JFNK

+which_ho_nonlinear = 0 # 0 = Picard, 1 = accelerated Picard

[parameters]

ice_limit = 1. # min thickness (m) for dynamics

@@ -33,12 +33,12 @@

title = parabolic dome test case using first-order dynamics

[CF input]

-name = dome-t6.0031.p001.nc

+name = dome.0031.p001.nc

time = 1

[CF output]

variables = thk usurf uvel vvel velnorm temp

frequency = 1

-name = dome-t6.0031.p001.out.nc

+name = dome.0031.p001.out.nc

xtype = double

Output Log

Dycore Type Processor Count Converged Iterations Avg. Iterations to Converge
Glissade 1 11 0.0

Output Log

Dycore Type Processor Count Converged Iterations Avg. Iterations to Converge
Glissade 1 11 30.0

p2

Bit for Bit

Variable Max Error Index of Max Error RMS Error Plot
velnorm 0.00000e+00 N/A 0.00000e+00
thk 0.00000e+00 N/A 0.00000e+00

Configuration Comparison

---

+++

@@ -23,7 +23,7 @@

which_ho_babc = 4 # 4 = no-slip at bed

which_ho_efvs = 2 # 2 = nonlinear eff. visc. w/ n=3

which_ho_sparse = 3 # 1 = SLAP GMRES, 3 = glissade parallel PCG, 4 = Trilinos for linear solver

-which_ho_nonlinear = 0 # 0 = Picard, 1 = JFNK

+which_ho_nonlinear = 0 # 0 = Picard, 1 = accelerated Picard

[parameters]

ice_limit = 1. # min thickness (m) for dynamics

@@ -33,12 +33,12 @@

title = parabolic dome test case using first-order dynamics

[CF input]

-name = dome-t8.0031.p002.nc

+name = dome.0031.p002.nc

time = 1

[CF output]

variables = thk usurf uvel vvel velnorm temp

frequency = 1

-name = dome-t8.0031.p002.out.nc

+name = dome.0031.p002.out.nc

xtype = double

Output Log

Dycore Type Processor Count Converged Iterations Avg. Iterations to Converge
Glissade 2 11 0.0

Output Log

Dycore Type Processor Count Converged Iterations Avg. Iterations to Converge
Glissade 2 11 30.0

p4

Bit for Bit

Variable Max Error Index of Max Error RMS Error Plot
velnorm 0.00000e+00 N/A 0.00000e+00
thk 0.00000e+00 N/A 0.00000e+00

Configuration Comparison

---

+++

@@ -23,7 +23,7 @@

which_ho_babc = 4 # 4 = no-slip at bed

which_ho_efvs = 2 # 2 = nonlinear eff. visc. w/ n=3

which_ho_sparse = 3 # 1 = SLAP GMRES, 3 = glissade parallel PCG, 4 = Trilinos for linear solver

-which_ho_nonlinear = 0 # 0 = Picard, 1 = JFNK

+which_ho_nonlinear = 0 # 0 = Picard, 1 = accelerated Picard

[parameters]

ice_limit = 1. # min thickness (m) for dynamics

@@ -33,12 +33,12 @@

title = parabolic dome test case using first-order dynamics

[CF input]

-name = dome-t3.0031.p004.nc

+name = dome.0031.p004.nc

time = 1

[CF output]

variables = thk usurf uvel vvel velnorm temp

frequency = 1

-name = dome-t3.0031.p004.out.nc

+name = dome.0031.p004.out.nc

xtype = double

Output Log

Dycore Type Processor Count Converged Iterations Avg. Iterations to Converge
Glissade 4 11 0.0

Output Log

Dycore Type Processor Count Converged Iterations Avg. Iterations to Converge
Glissade 4 11 30.0

p8

Bit for Bit

Variable Max Error Index of Max Error RMS Error Plot
velnorm 0.00000e+00 N/A 0.00000e+00
thk 0.00000e+00 N/A 0.00000e+00

Configuration Comparison

---

+++

@@ -23,7 +23,7 @@

which_ho_babc = 4 # 4 = no-slip at bed

which_ho_efvs = 2 # 2 = nonlinear eff. visc. w/ n=3

which_ho_sparse = 3 # 1 = SLAP GMRES, 3 = glissade parallel PCG, 4 = Trilinos for linear solver

-which_ho_nonlinear = 0 # 0 = Picard, 1 = JFNK

+which_ho_nonlinear = 0 # 0 = Picard, 1 = accelerated Picard

[parameters]

ice_limit = 1. # min thickness (m) for dynamics

@@ -33,12 +33,12 @@

title = parabolic dome test case using first-order dynamics

[CF input]

-name = dome-t6.0031.p008.nc

+name = dome.0031.p008.nc

time = 1

[CF output]

variables = thk usurf uvel vvel velnorm temp

frequency = 1

-name = dome-t6.0031.p008.out.nc

+name = dome.0031.p008.out.nc

xtype = double

Output Log

Dycore Type Processor Count Converged Iterations Avg. Iterations to Converge
Glissade 8 11 0.0

Output Log

Dycore Type Processor Count Converged Iterations Avg. Iterations to Converge
Glissade 8 11 30.0

References

LIVVkit is an open source project licensed under a BSD 3-clause License. We ask that you please acknowledge LIVVkit in any work it is used or supports. In any corresponding published work, please cite:

1
K. J. Evans, J. H. Kennedy, D. Lu, M. M. Forrester, S. Price, J. Fyke, A. R. Bennett, M. J. Hoffman, I. Tezaur, C. S. Zender, and M. Vizca\'ıno. Livvkit 2.1: automated and extensible ice sheet model validation. Geoscientific Model Development, 12(3):1067–1086, 2019. URL: https://www.geosci-model-dev.net/12/1067/2019/, doi:10.5194/gmd-12-1067-2019.
2
Joseph H. Kennedy, Andrew R. Bennett, Katherine J. Evans, Stephen Price, Matthew Hoffman, William H. Lipscomb, Jeremy Fyke, Lauren Vargo, Adrianna Boghozian, Matthew Norman, and Patrick H. Worley. Livvkit: an extensible, python-based, land ice verification and validation toolkit for ice sheet models. Journal of Advances in Modeling Earth Systems, 9(2):854–869, 2017. doi:10.1002/2017MS000916.

LIVVkit was developed under PISCEES, a BER/ASCR SciDAC Earth System Modeling project.

Contact us on github