Added clean sac
Readme
Typo
Correction too boundary terms for computation of hyperviscosity. Too much branching as a result of if statements resulted in a small error at the edge of the mesh
Minor update enabling simulation to start from any step number
Corrections to compile using SAC_3D and to make the same as the smau and gsmaug distributions
Corrections to compile using SAC_3D
MPI version works for 2D for configurations with ni==nj
(square grids) and square processor array.

Modifications here make the ni<>nj configs possible using MPI.
There is a need to check commneted lines for mpirecvbuffer, mpirecvbuffermod, 
mpibuffer2var and mpibuffer2varmod.
MPI versions of Readwrite routines adjusted so that ghost cells not written for dimensions where number of processors>1.
Buffer sizes for different directions modified for each direction (for ni<>nj<>nk)
minor correction npe[] replaced by pnpe[]
Setting boundary conditions correctly for multi-processor configurations
Tested version of MPI implementation of smaug which allows computation of a rectilinear grid.
Tests with square configuration using OZT were passed for this version.
beta-release: correction to cuinit and transfer routines from gpu to mpi buffer.
Tests revealed correct operation of MIP messaging for rectangulat 2d  configurations.
Still need to check that square configuration test work OK.

Further verification of hyperdiffusion terms is required
Tested hyperdiffusion for parallel sac with rectangular grids.

There is an issue with the boundary terms for hyperdifvisc1l.cu and hyperdifvis1r.cu

The computation of the boundary contribution in hyperdifvisc1ar/l_parallel

Needs to handle the boundary condition correctly

Conditional statements on the boundary need to be checked in the host code
Boundary kernel is called separately
Current version computes boundary condition for mpi boundaries and not just the edges of the box!
Working MPI version.
Introducing routines using GPU direct RDMA
Reomved USE_SAC3D from #define statements and replaced with USE_SAC_3D
Moved cuda functionality for inter gpu comms into mpiutil