Introduce masks for both full and half levels
Use the correct filter at the correct place
Remove the filter interpolations from the functions...
Cleaned up the statistics
Introduce the filter into the diffusion calculation for LES
Add the proper horizontal interpolations to the diff function