An isopycnal version of SEOM-2D has been developed. It does not allow 
layer-outcropping nor intersecting the bathymetry. More to come...