Comments revised/added, formatting improved, fluxes for u,v, and scalars are explicitly set to zero at nzt+1, fluxes of w-component are now calculated only until nzt-1 (Prognostic equation for w-velocity component ends at nzt-1)
