New test case: netcdf_read_free_surface
This creates a valid netCDF file containing a free surface height in the domain x,y in [-1.2,1.2]. This is read in by Fluidity as a free surface and applied as an initial condition on pressure. From this, Fluidity calculates the free surface as a diagnostic, which is provided in the outputted VTU file, alongside a field containing the analytic function used to create the netCDF file in the first place (here described in height.py). The test compares the difference between these two fields, which should be small, due to the interpolation of the netCDF file.

In summary this case tests:
 - NetCDF reader
 - Reading of a free surface and conversion into a pressure initial condition
 - Creation of the free surface diagnostic from the pressure field (including its extrapolation down from the top surface)
 - Creation of prescribed scalar fields from Python