y5d_update(k) = 0. _d 0
y5d_update(k) = y5d_prime(k)
& - y5d_update(k+1)*d5d_prime(k)
& - y5d_update(k+2)*e5d_prime(k)
y5d(i,j,k) = y5d_update(k)
_RL y5d_update(Nr)