output switched to netCDF and functions commented
discretized laplacian