PONflux_u = 0. _d 0
PONflux_l = (PONflux_u+N_spm(i,j,k)*drF(k)
N_reminp(i,j,k) = (PONflux_u + N_spm(i,j,k)*drF(k)
N_reminp(i,j,k) = PONflux_u*recip_drF(k)
PONflux_u = PONflux_l
_RL PONflux_u