_RL frace, rs, hq
frace=1. _d 0
frace=frac_energy
frace=0. _d 0
if (frac_energy.eq.0. _d 0) frace=0. _d 0
etop = (1. _d 0-frace)*sHeating * dt
etope = frace*sHeating * dt
ebote = frace*ebot