eta_hp = round(heat_b[region]["sys_heat"] * heat_b[region]["hp"], 2)