* KYNBB was updated to resolve an undefined behaviour when input energy array is
  changed.

* extended corona geometry was also added.

* in the KYNBB model, the colour correction is computed according to
  Done et al. (2012) MNRAS 420, 1848 for parameter f_col=-1

* xside - flat space approximation above 1000rg added
* kynbb & kynphebb - radius can be up to 100000 rg or infinite
* kynbb - the accretion rate definition in Medd added
* Q/I and U/I added as an output for polarization in kynbb, kynphebb and
  kynlpcr
* kynlpcr - the orientation angle added

according to the XFLT0001 keyword.

* Add definition of the KYNBB and KYNPHEBB model output according to the
  XFLT0001 keyword of the SPECTRUM extension of the data file.
* Speed up KYNBB and KYNPHEBB model when using the same parameters (e.g. when
  using several data sets of the same source) by remembering the output from
  the last computation.

* Add system orientation angle.
* Add negative inclination for opposite system rotation.

* Add obscuration of primary X-ray source
* Add q^2 from axis to the observer to lamp tables to compute impact parameter
  beta of the lamp needed for primary obscuration and change of observed
  polarisation angle of the primary radiation
* Rename lamp tables to 'KBHlamp80.fits'
* Add polarisation in KYNlpcr with polarised primary
* Add polarisation in KYNbb and KYNphebb
* Create 'goosmann.fits' with local polarisation properties of thermal disc
  radiation for different tau computed with Stokes by Rene Goosmann
* In KYNbb and KYNphebb, compute contribution of outer disc only if
  rout >= 1000 GM/c^2
* Add possibility to use Garcia's xillverD-4.fits tables where density
  dependence of the reflection spectra is included (for radial power-law profile
  of the disc density and for constant ionisation of the disc)
* Add photon flux as output in XSPEC for parameter Stokes=1 when polarisation
  computations are switched on
* Update documentation

* Change the definition of Np:Nr and luminosity towards the observer or towards
  the disc in KYNlpcr, KYNrefionx and KYNxillver
* Add direct definition of ionisation of the disc in KYNrefionx and KYNxillver,
  i.e. instead of indirect definition by specifying the density
* Update definition of shifted Ecut in KYNrefionx
* Add possibility to define cut-off energy in the observer's reference frame
  in KYNxillver
* Update KYNxillver for the newest XILLVER tables (note the definition order of
  the previous ones has been changed)
* Correct error in interpolation of local emission angle in KYNxillver
* Correct small error in using local energy in incomplete Gamma function
  in KYNxillver