extract_heatflux_profile_from_DL always uses timepoint=1
as in title, this is problematic when the function is called from within another where the expected and desired behaviour is to return mutliple profiles. Suggested change: Modify the default of timepoint such that the largest available range of data is returned.