extended system and ploting to 3d for cube geometry, added 2nd order operator splitting