new version of tign_lfn_interpolation subroutine calculating 9 endpoints of 4 subcells