== Table 1. Debugging status ==

{{{
#!html
<table border="1" style="border-collapse:collapse; border-width:2px;" width="700px">
<tr>
<th>Dimension</th> <th colspan="4">2D</th> <th colspan="4">3D</th>
</tr>
<tr>
<th>Geometry</th> <th colspan="2">Ice-sheet</th> <th colspan="2">Ice-Shelf</th> <th colspan="2">Ice-sheet</th> <th colspan="2">Ice-Shelf</th>
</tr>
<tr>
<th>Parallel</th> <th>Serial</th> <th>Parallel</th> <th>Serial</th> <th>Parallel</th> <th>Serial</th> <th>Parallel</th> <th>Serial</th> <th>Parallel</th>
</tr>

<tr>
<th>Surface Slope</th>
<td bgcolor="#ddffdd"></td><td bgcolor="#ddffdd"></td><td bgcolor="#ddffdd"></td><td bgcolor="#ddffdd"></td><td bgcolor="#ddffdd"></td><td bgcolor="#ddffdd"></td><td bgcolor="#ddffdd"></td><td bgcolor="#ddffdd"></td>
</tr>

<tr>
<th>Hutter</th>
<td bgcolor="#ddffdd"></td><td bgcolor="#ddffdd"></td><td bgcolor="#cccccc"></td><td bgcolor="#cccccc"></td><td bgcolor="#ddffdd"></td><td bgcolor="#ddffdd"></td><td bgcolor="#cccccc"></td><td bgcolor="#cccccc"></td>
</tr>

<tr>
<th>MacAyeal</th>
<td bgcolor="#ffdddd"></td><td bgcolor="#ffdddd"></td><td bgcolor="#ddffdd"></td><td bgcolor="#ddffdd"></td><td bgcolor="#ffdddd"></td><td bgcolor="#ffdddd"></td><td bgcolor="#ddffdd"></td><td bgcolor="#ddffdd"></td>
</tr>

<tr>
<th>Pattyn</th>
<td bgcolor="#cccccc"></td><td bgcolor="#cccccc"></td><td bgcolor="#cccccc"></td><td bgcolor="#cccccc"></td><td bgcolor="#ffdddd"></td><td bgcolor="#ffdddd"></td><td bgcolor="#ffdddd"></td><td bgcolor="#ffdddd"></td>
</tr>

<tr>
<th>Stokes</th>
<td bgcolor="#cccccc"></td><td bgcolor="#cccccc"></td><td bgcolor="#cccccc"></td><td bgcolor="#cccccc"></td><td bgcolor="#ffff00"></td><td bgcolor="#ffff00"></td><td bgcolor="#ffff00"></td><td bgcolor="#ffff00"></td>
</tr>

<tr>
<th>Prognostic</th>
<td bgcolor="#ddffdd"></td><td bgcolor="#ddffdd"></td><td bgcolor="#ddffdd"></td><td bgcolor="#ddffdd"></td><td bgcolor="#ffff00"></td><td bgcolor="#ffff00"></td><td bgcolor="#ffff00"></td><td bgcolor="#ffff00"></td>
</tr>

<tr>
<th>Prognotic 2</th>
<td bgcolor="#ffdddd"></td><td bgcolor="#ffdddd"></td><td bgcolor="#ffdddd"></td><td bgcolor="#ffdddd"></td><td bgcolor="#cccccc"></td><td bgcolor="#cccccc"></td><td bgcolor="#cccccc"></td><td bgcolor="#cccccc"></td>
</tr>

<tr>
<th>Balanced thickness</th>
<td bgcolor="#ddffdd"></td><td bgcolor="#ddffdd"></td><td bgcolor="#ddffdd"></td><td bgcolor="#ddffdd"></td><td bgcolor="#ffff00"></td><td bgcolor="#ffff00"></td><td bgcolor="#ffff00"></td><td bgcolor="#ffff00"></td>
</tr>

<tr>
<th>Balanced thickness 2</th>
<td bgcolor="#ffff00"></td><td bgcolor="#ffff00"></td><td bgcolor="#ffff00"></td><td bgcolor="#ffff00"></td><td bgcolor="#ffff00"></td><td bgcolor="#ffff00"></td><td bgcolor="#ffff00"></td><td bgcolor="#ffff00"></td>
</tr>

<tr>
<th>Balanced velocities</th>
<td bgcolor="#ffff00"></td><td bgcolor="#ffff00"></td><td bgcolor="#ffff00"></td><td bgcolor="#ffff00"></td><td bgcolor="#ffff00"></td><td bgcolor="#ffff00"></td><td bgcolor="#ffff00"></td><td bgcolor="#ffff00"></td>
</tr>

<tr>
<th>Thermal steady</th>
<td bgcolor="#cccccc"></td><td bgcolor="#cccccc"></td><td bgcolor="#cccccc"></td><td bgcolor="#cccccc"></td><td bgcolor="#ffff00"></td><td bgcolor="#ffff00"></td><td bgcolor="#ffff00"></td><td bgcolor="#ffff00"></td>
</tr>

<tr>
<th>Thermal Transient</th>
<td bgcolor="#cccccc"></td><td bgcolor="#cccccc"></td><td bgcolor="#cccccc"></td><td bgcolor="#cccccc"></td><td bgcolor="#ffff00"></td><td bgcolor="#ffff00"></td><td bgcolor="#ffff00"></td><td bgcolor="#ffff00"></td>
</tr>

<tr>
<th>Transient</th>
<td bgcolor="#ffdddd"></td><td bgcolor="#ffdddd"></td><td bgcolor="#ffdddd"></td><td bgcolor="#ffdddd"></td><td bgcolor="#ffff00"></td><td bgcolor="#ffff00"></td><td bgcolor="#ffff00"></td><td bgcolor="#ffff00"></td>
</tr>

<tr>
<th>CM B</th>
<td bgcolor="#ffff00"></td><td bgcolor="#ffff00"></td><td bgcolor="#ffff00"></td><td bgcolor="#ffff00"></td><td bgcolor="#ffff00"></td><td bgcolor="#ffff00"></td><td bgcolor="#ffff00"></td><td bgcolor="#ffff00"></td>
</tr>

<tr>
<th>CM drag</th>
<td bgcolor="#ffff00"></td><td bgcolor="#ffff00"></td><td bgcolor="#ffff00"></td><td bgcolor="#ffff00"></td><td bgcolor="#ffff00"></td><td bgcolor="#ffff00"></td><td bgcolor="#ffff00"></td><td bgcolor="#ffff00"></td>
</tr>

<th>CM steady</th>
<td bgcolor="#cccccc"></td><td bgcolor="#cccccc"></td><td bgcolor="#cccccc"></td><td bgcolor="#cccccc"></td><td bgcolor="#ffff00"></td><td bgcolor="#ffff00"></td><td bgcolor="#ffff00"></td><td bgcolor="#ffff00"></td>
</tr>

</table>
}}}

Legend:

{{{
#!html
<table>
<tr><td bgcolor="#ddffdd" width="20px"></td><td>Success (runs and matches archive)</td></tr>
<tr><td bgcolor="#ffdddd"></td><td>Error (does not match archive)</td></tr>
<tr><td bgcolor="#ffff00"></td><td>Failure (does not run)</td></tr>
<tr><td bgcolor="#cccccc"></td><td>Not applicable</td></tr>
</table>
}}}


== 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