Fixed the boundary conditions for zsort in order to be able to calculate
the terms in Winters1995
Add the diffusion terms according to Winters. Only one term to go now