get cycle_numbers in addition to float_index

[rampatels]

In my regional processes-specific study I need to work with individual profiles. When I try to extract a specific profile from the output of the load_data function, it fails. Because cycle number and profile index is differently interpreted in the function output. It would be great if you could add an optional output for the cycle number in addition to profile_index.

Thank you so much.

At this stage, I am simply recreating profile_index for the loaded data.