The code above won’t compile so I am not sure how you ran the simulation and check the monitor values. PR1O3 and RE are scalar variables, not arrays. You would need something like
ReactionRates(IJK,7) = PR1O3
and
ReactionRates(IJK,8) = RE
You can write a custom gas-solid heat transfer in subroutine USR_PROP_Gama, found in usr_properties.f