Icarus

Features:

- Added some of the compliance for Python 3.x. However, some libraries such as scipy.weave are no longer supported under 3.x, which makes things a bit awkward.


Bug fixes:

- Fixed one further assumption about the was that coschi, the angle between the incident irradiating flux and the normal to the surface, is calculation. The exact calculation is now used. This has an effect especially on shorter orbital period systems in which the companion is large, and the terminator does not lie on the meridian passing through the centre of mass, but towards the nose of the star.

Features:

- This major release has now changed the directory structure to be in line with the standard Python packaging. It includes a ``setup.py`` to enable the installation of the package to the standard Python library. There is also a PyPI distribution now. Many thanks to **muddassir235** for the pull request!


Bug fixes:

- _Photometry/Photometry.py/Photometry_temperature_: Updated the class so is now compatible with the new version of the _Photometry_ class.

Bug fixes:

- _Atmosphere/Atmo.py/AtmoGridSpec.Get_flux_doppler_: There was an issue with the interpolation along the wavelength axis which has now been resolved. It was simply not computed correctly.
- _Core/Star_base_: Fixed a problem in the Velocity_surface function which calculates the surface velocity.
- _Atmosphere/Atmo_spectro_IRTF.py_: Fixed the conversion of the fluxes, which must be stored in the grid and log(F_lambda), but the grid fluxes store them at log10(F_lambda).
- _Utils/Grid.py_: _Interp_doppler_savememory_ and _Inter_doppler_ were both not exponentiating the fluxes before summing them over the surface elements. One needs to bear in mind that the grid values are log(F_lambda).

Features:

- _Atmosphere/Atmo.py/AtmoGridSpec_ is a new class to handle atmosphere grids in the spectral domain.

- _Utils/Filter.py/Band_integration_ is modified in order to handle masking parts of the spectral range. The default method is Simpson's rule, but trapezoid is also implemented and recommended for the masking.

Features:

- _Atmosphere/Atmo.py/AtmoGrid_

- New _Fill_nan_ function to fill non-existing values.
- New _SubGrid_ function to return a sub-grid of the atmosphere grid.
- New _IsFinite_ function to return a boolean view of existing values

- _Utils/Binary.py_: New function to calculate the approximate Roche lobe radius using the Eggleton formula.
- _Utils/Filter.py/Load_filter_: New option that can also return the min/max bounds of the filter, in addition to the filter interpolator.


Bug fixes:

- _Core/Star.py/Make_surface_

- Changed the way that coschi -- angle between the irradiation source and the normal to the surface -- is calculated. It used to be computed using the spherical approximation (i.e that the normal to the surface was roughly the same as the vector direction from the centre of mass of the star to the surface). This worked well for low filling factors, but breaks down at larger filling factor. It is better to use the exact value of the normal to the surface, which is the gradient vector (already calculated to solve the equipotential surface). Because the irradiation source is along the x-direction, coschi simply corresponds to the x-component of the gradient vector. __Many thanks to Roger Romani and Nicholas Sanchez at Stanford University for highlighting this issue and providing pieces of codes.__

- Small bug fix in _Photometry/Photometry.py_ regarding a typo in a variable name for the _full_output_ option in the _Plot_model_ function.

Features:

- In the _Star_base_ class, change in the way _Flux_ is calculated in order to add the possibility to account for the projection effect in order to scale the flux to physical values. There wasn't a need for it before as the fluxes were always used in a normalised way (e.g. fitting a non-flux calibrated spectrum), but this would make things fail if one was to fit flux calibrated data. The issue wouldn't arise with magnitudes, which already included the projection.

Bug fixes:

- Fixed an issue with the conversion of mag to/from flux in the _Photometry_ module. The conversion was expecting to find a _flux0_ parameter, whereas the zeropoint is now incorporated as _zp_ in the atmosphere files.
- The _Utils.import_modules_ has a set of constants defined under the _cts_ namespace. Converted some lower case constants into upper case.
- Many changes made to the _Photometry_ module in order to handle working in the flux space properly (as opposed to the magnitude space). Now all the plotting can be done in the flux space. Also, everything related to dealing with the band offset has been rebrewed in order to simplify function calls and make everything more intuitive. This includes changes to the optional arguments and their names. Please make sure to look at the docstrings for _Calc_chi2_, _Get_flux_, _Get_flux_theoretical_, _Plot_.