Table 1. Debugging status

Dimension	2D				3D
Geometry	Ice-sheet		Ice-Shelf		Ice-sheet		Ice-Shelf
Parallel	Serial	Parallel	Serial	Parallel	Serial	Parallel	Serial	Parallel
Surface Slope
Hutter
MacAyeal
Pattyn
Stokes
Prognostic
Prognotic 2
Balanced thickness
Balanced thickness 2
Balanced velocities
Thermal steady
Thermal Transient
Transient
CM B
CM drag
CM steady

Legend:

	Success (runs and matches archive)
	Error (does not match archive)
	Failure (does not run)
	Not applicable

Bugs to be fixed

• uppernode (field does not work for DG for example) -> Get rid of it. Use only upperelement
• Update element: node indexing assumes P1 elements!!! (See Penta for example)
• control method: RheologyB is found in inputs but it is a material property
• SpawnTria assumes that the Penta has 6 nodes and Tria 3. (What if P2?)
• NumberOfDofsPerNode (in parameters) not updated by SetCurrentConfiguration -> always the same (see NodesDofx.cpp)
• if no nodes in a model -> max_sid++ => same result if 0 or 1 node (in Container/Nodes.cpp l264)
• Discontinuous Galerkin: Numerical flux is linked to only one element. If Vx and Vy are discontinuous, we need to create 2 numerical flux: One one the right and one on the left to take into account the discontinuity in velocity