Vertical fluxes of w are now set to zero at nzt and nzt+1, setting of advection flags for fluxes in z-direction revised, comments extended
