accretionStep not supported in 2D

Simulation::removeAccretedGas does a h+1 access even when flat.

Also, even in 3D, j_max should be initialized with ni-1 for h+1 to be valid. but we do have small regressions when doing so.