== 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="#dfd"></td><td bgcolor="#dfd"></td><td bgcolor="#dfd"></td><td bgcolor="#dfd"></td><td bgcolor="#dfd"></td><td bgcolor="#dfd"></td><td bgcolor="#dfd"></td><td bgcolor="#dfd"></td></tr>
<tr><th>Hutter              </th><td bgcolor="#dfd"></td><td bgcolor="#dfd"></td><td bgcolor="#ccc"></td><td bgcolor="#ccc"></td><td bgcolor="#dfd"></td><td bgcolor="#dfd"></td><td bgcolor="#ccc"></td><td bgcolor="#ccc"></td></tr>
<tr><th>Hutter/MacAyeal     </th><td bgcolor="#dfd"></td><td bgcolor="#dfd"></td><td bgcolor="#ccc"></td><td bgcolor="#ccc"></td><td bgcolor="#dfd"></td><td bgcolor="#dfd"></td><td bgcolor="#ccc"></td><td bgcolor="#ccc"></td></tr>
<tr><th>MacAyeal            </th><td bgcolor="#dfd"></td><td bgcolor="#dfd"></td><td bgcolor="#dfd"></td><td bgcolor="#dfd"></td><td bgcolor="#dfd"></td><td bgcolor="#dfd"></td><td bgcolor="#dfd"></td><td bgcolor="#dfd"></td></tr>
<tr><th>Pattyn              </th><td bgcolor="#ccc"></td><td bgcolor="#ccc"></td><td bgcolor="#ccc"></td><td bgcolor="#ccc"></td><td bgcolor="#ff9"></td><td bgcolor="#dfd"></td><td bgcolor="#ff9"></td><td bgcolor="#dfd"></td></tr>
<tr><th>Stokes              </th><td bgcolor="#ccc"></td><td bgcolor="#ccc"></td><td bgcolor="#ccc"></td><td bgcolor="#ccc"></td><td bgcolor="#ff9"></td><td bgcolor="#ff9"></td><td bgcolor="#ff9"></td><td bgcolor="#ff9"></td></tr>
<tr><th>Prognostic          </th><td bgcolor="#dfd"></td><td bgcolor="#dfd"></td><td bgcolor="#dfd"></td><td bgcolor="#dfd"></td><td bgcolor="#dfd"></td><td bgcolor="#dfd"></td><td bgcolor="#dfd"></td><td bgcolor="#dfd"></td></tr>
<tr><th>Prognotic 2         </th><td bgcolor="#dfd"></td><td bgcolor="#dfd"></td><td bgcolor="#dfd"></td><td bgcolor="#dfd"></td><td bgcolor="#ccc"></td><td bgcolor="#ccc"></td><td bgcolor="#ccc"></td><td bgcolor="#ccc"></td></tr>
<tr><th>Balanced thickness  </th><td bgcolor="#dfd"></td><td bgcolor="#dfd"></td><td bgcolor="#dfd"></td><td bgcolor="#dfd"></td><td bgcolor="#dfd"></td><td bgcolor="#dfd"></td><td bgcolor="#dfd"></td><td bgcolor="#dfd"></td></tr>
<tr><th>Balanced thickness 2</th><td bgcolor="#dfd"></td><td bgcolor="#dfd"></td><td bgcolor="#dfd"></td><td bgcolor="#dfd"></td><td bgcolor="#ccc"></td><td bgcolor="#ccc"></td><td bgcolor="#ccc"></td><td bgcolor="#ccc"></td></tr>
<tr><th>Thermal steady      </th><td bgcolor="#ccc"></td><td bgcolor="#ccc"></td><td bgcolor="#ccc"></td><td bgcolor="#ccc"></td><td bgcolor="#dfd"></td><td bgcolor="#dfd"></td><td bgcolor="#dfd"></td><td bgcolor="#dfd"></td></tr>
<tr><th>Thermal Transient   </th><td bgcolor="#ccc"></td><td bgcolor="#ccc"></td><td bgcolor="#ccc"></td><td bgcolor="#ccc"></td><td bgcolor="#dfd"></td><td bgcolor="#dfd"></td><td bgcolor="#dfd"></td><td bgcolor="#dfd"></td></tr>
<tr><th>Transient           </th><td bgcolor="#dfd"></td><td bgcolor="#dfd"></td><td bgcolor="#dfd"></td><td bgcolor="#dfd"></td><td bgcolor="#ff9"></td><td bgcolor="#ff9"></td><td bgcolor="#ff9"></td><td bgcolor="#ff9"></td></tr>
<tr><th>Steady state        </th><td bgcolor="#ccc"></td><td bgcolor="#ccc"></td><td bgcolor="#ccc"></td><td bgcolor="#ccc"></td><td bgcolor="#ff9"></td><td bgcolor="#ff9"></td><td bgcolor="#ff9"></td><td bgcolor="#ff9"></td></tr>
<tr><th>CM rheology         </th><td bgcolor="#dfd"></td><td bgcolor="#dfd"></td><td bgcolor="#dfd"></td><td bgcolor="#dfd"></td><td bgcolor="#dfd"></td><td bgcolor="#dfd"></td><td bgcolor="#dfd"></td><td bgcolor="#dfd"></td></tr>
<tr><th>CM basal drag       </th><td bgcolor="#dfd"></td><td bgcolor="#dfd"></td><td bgcolor="#dfd"></td><td bgcolor="#dfd"></td><td bgcolor="#ff9"></td><td bgcolor="#ff9"></td><td bgcolor="#ff9"></td><td bgcolor="#ff9"></td></tr>
<tr><th>CM steady state     </th><td bgcolor="#ccc"></td><td bgcolor="#ccc"></td><td bgcolor="#ccc"></td><td bgcolor="#ccc"></td><td bgcolor="#dfd"></td><td bgcolor="#dfd"></td><td bgcolor="#dfd"></td><td bgcolor="#dfd"></td></tr>

</table>
}}}

Legend:

{{{
#!html
<table>
<tr><td bgcolor="#dfd" width="20px"></td><td>Success (runs and matches archive)</td></tr>
<tr><td bgcolor="#ff9" width="20px"></td><td>Minor Error (mismatch on the order of magnitude of success criterion)</td></tr>
<tr><td bgcolor="#fdd" width="20px"></td><td>Error (does not match archive)</td></tr>
<tr><td bgcolor="#f00" width="20px"></td><td>Failure (does not run)</td></tr>
<tr><td bgcolor="#ccc" width="20px"></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
 - 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
 - How to create velocity results from discontinuous velocity (Hutter/MacAyeal create a different velocity at the interface)
 - We should use the depth-averaged velocity of Hutter at the boundary of MacAyeal (in Hutter/MacAyeal models) instead of the velocity on the lower layer
 - when do you update B (beginning of thermal steady ? thermal transient ? steadystate ? transient ?) and what temperature do you use (temp ? temp_average ?)
 - use temp_average all the time for MacAyeal
 - CM Stokes, DU is garbage. Check CM B ice shelf (results should be the same as Pattyn)