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
Thermal steady
Thermal Transient
Transient
Steady state
CM rheology
CM basal drag
CM steady state

Legend:

	Success (runs and matches archive)
	Minor Error (mismatch on the order of magnitude of success criterion)
	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
• in updateVertexPosition each vertices is looking for its new z using their VERTICES DOF in a vector built by the elements using their NODES DOFS