Add more robust code to read in a sigma.dat file (skipping comments, blank lines etc.). Also add optional argument to sample a larger number of points from the bathymetry