$$\frac{\partial }{\partial t}R(X,Y,Z,W)=(\Delta R)(X,Y,Z,Y)+\sum_{pq}R(X,Y,e_p,e_q)R(Z,W,e_p,e_q)+2\sum_{pq}R(X,e_p,Z,e_q)R(Y,e_p,W,e_q)-2\sum_{pq}R(X,e_p,W,e_q)R(Y,e_p,Z,e_q)$$$$-\sum_k\mathrm{Ric}(X,e_k)R(e_k,Y,Z,W)-\sum_k\mathrm{Ric}(Y,e_k)R(E,e_k,Z,W)-\sum_k\mathrm{Ric}(Z,e_k)R(X,Y,e_k,W)-\sum_k\mathrm{Ric}(W,e_k)R(X,Y,Z,e_k)$$