link ../phylmd/thermcell_closure.F90