link ../phylmd/thermcell_init.F90