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 |
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-t5.0062.p001.nc
+name = dome.0062.p001.nc
time = 1
[CF output]
variables = thk usurf uvel vvel velnorm temp
frequency = 1
-name = dome-t5.0062.p001.out.nc
+name = dome.0062.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 |
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-t4.0062.p004.nc
+name = dome.0062.p004.nc
time = 1
[CF output]
variables = thk usurf uvel vvel velnorm temp
frequency = 1
-name = dome-t4.0062.p004.out.nc
+name = dome.0062.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 |
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-t0.0124.p001.nc
+name = dome.0124.p001.nc
time = 1
[CF output]
variables = thk usurf uvel vvel velnorm temp
frequency = 1
-name = dome-t0.0124.p001.out.nc
+name = dome.0124.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 | 28.0 |
p16
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-t5.0124.p016.nc
+name = dome.0124.p016.nc
time = 1
[CF output]
variables = thk usurf uvel vvel velnorm temp
frequency = 1
-name = dome-t5.0124.p016.out.nc
+name = dome.0124.p016.out.nc
xtype = double
Output Log
Dycore Type | Processor Count | Converged Iterations | Avg. Iterations to Converge |
---|---|---|---|
Glissade | 16 | 11 | 0.0 |
Output Log
Dycore Type | Processor Count | Converged Iterations | Avg. Iterations to Converge |
---|---|---|---|
Glissade | 16 | 11 | 28.0 |
p256
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-t0.0124.p256.nc
+name = dome.0124.p256.nc
time = 1
[CF output]
variables = thk usurf uvel vvel velnorm temp
frequency = 1
-name = dome-t0.0124.p256.out.nc
+name = dome.0124.p256.out.nc
xtype = double
Output Log
Dycore Type | Processor Count | Converged Iterations | Avg. Iterations to Converge |
---|---|---|---|
Glissade | 256 | 11 | 0.0 |
Output Log
Dycore Type | Processor Count | Converged Iterations | Avg. Iterations to Converge |
---|---|---|---|
Glissade | 256 | 11 | 28.0 |
p64
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-t2.0248.p064.nc
+name = dome.0248.p064.nc
time = 1
[CF output]
variables = thk usurf uvel vvel velnorm temp
frequency = 1
-name = dome-t2.0248.p064.out.nc
+name = dome.0248.p064.out.nc
xtype = double
Output Log
Dycore Type | Processor Count | Converged Iterations | Avg. Iterations to Converge |
---|---|---|---|
Glissade | 64 | 11 | 0.0 |
Output Log
Dycore Type | Processor Count | Converged Iterations | Avg. Iterations to Converge |
---|---|---|---|
Glissade | 64 | 11 | 28.0 |
p256
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.0248.p256.nc
+name = dome.0248.p256.nc
time = 1
[CF output]
variables = thk usurf uvel vvel velnorm temp
frequency = 1
-name = dome-t6.0248.p256.out.nc
+name = dome.0248.p256.out.nc
xtype = double
Output Log
Dycore Type | Processor Count | Converged Iterations | Avg. Iterations to Converge |
---|---|---|---|
Glissade | 256 | 11 | 0.0 |
Output Log
Dycore Type | Processor Count | Converged Iterations | Avg. Iterations to Converge |
---|---|---|---|
Glissade | 256 | 11 | 28.0 |
p256
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-t2.0496.p256.nc
+name = dome.0496.p256.nc
time = 1
[CF output]
variables = thk usurf uvel vvel velnorm temp
frequency = 1
-name = dome-t2.0496.p256.out.nc
+name = dome.0496.p256.out.nc
xtype = double
Output Log
Dycore Type | Processor Count | Converged Iterations | Avg. Iterations to Converge |
---|---|---|---|
Glissade | 256 | 11 | 0.0 |
Output Log
Dycore Type | Processor Count | Converged Iterations | Avg. Iterations to Converge |
---|---|---|---|
Glissade | 256 | 11 | 27.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.