nondim.S_sl was calculated with min(f(:))as opposed to slope calculated from bathy.S_sl where, i used f0. Now I use f0 everywhere