Add multilayer periodic bathymetry test script

Really does not work due to shallow water assumptions failing, should
try to do in the 2d case though.