Fix entropy residual computation in SRHD
There was a problem where the entropy production was wrongly computed when using the entropy AV as a HRSC.
This was because both the entropy TCI and the entropy AV computed the residuals themselves and both reused the the flx_Em1
buffer.
This is fixed by moving the computation of the entropy production to callback_equation
.