DO k=kshelf+1,kBottom(ijd,bi,bj),downward
INTEGER downward
downward = 1
IF ( usingPCoords ) downward = -1
downward = -1
DO k = kshelf,kdeep,downward