Observations ============ * It is difficult to converge this calculation for all models and for a large angle phi0. This is why the angle phi0 and the division of the time steps are chosen differently for each modeling. * Models POU_D_EM and POU_D_TGM are frankly flawed on this test (error > 100%). This documentation is deliberately brief.