Add interpolate function from SOWFA