link ../phylmd/thermcell_qsat.F90