How to output the heat of reaction in cell for a hetergeneous reaction

Hi, I have simulated the hetergenous reaction between NiO and methane. I have output the heat of reaction for solid phase as shown here.
image

I also want to output the heat of reaction in each cell (llHORg) to identify whether it equalls to zero as shown in source code

Here are the realted file. Thank you so much for your help
Uploading: 0605.mfx…
calc_rrates_des.f (9.4 KB)
geometry_0003.stl (34.8 KB)
NiO_psd.txt (125 Bytes)

usr_rates_des.f (5.9 KB)