Feature external coupled

Code infrastructure that help with adding arbitrary additional loops into solvers. One example is for cht energy coupling.