midway implementing s co-ordinate vorticity diagnostics
full vorticity balance in s co-ordinates implemented. Doesn't balance exactly yet. I should expect balance to happen only with very frequent output so that tendency term is well estimated