Atmospheric models for tailored simulations
According to our current understanding, the products of the create_atmospheric_model workflow will be the following 2 (at least for release 0):
- a reference atmosphere (Corsika input) with all necessary information for the primary interaction and the shower development
- a light extinction profile, currently including only the contribution from Rayleigh scattering.
The reference atmosphere consists of 2 parts, ECMWF data up to ~45km a.s.l. and upper atmosphere data from 45 -> 120 km a.s.l. We need to find a way to connect the contemporary produced profile via ECMWF data with the upper atmosphere. To do that we need the season and the observatory site. The site can be taken directly from the Observatory object. Concerning the season, I see 2 ways:
- by the time we will produce the atmospheric model, the cat-B calibration for that specific date would have been completed. We can pick up the season from the database.
- we include select_reference_season workflow (using ECMWF data) within the create_atmosphere_profile tool. I prefer this method, it guarantees more consistent results. The returned season and the name attribute of the Observatory object will be used for selecting the correct upper atmosphere file. (we can also think of a current_season attribute for the Observatory class)
Concerning the light extinction, the create_rayleigh_extinction_file function should be also included in the create_atmosphere_profile tool.