Indirect indexing for calculating vertical fluxes close to boundaries is only used for loop indices where it is really necessary
