NMSSMCALC version 3.00, 27 March 2019:
These changes concern both the complex and the real version of NMSSMCALC. They have
been made by Thi Nhung Dao.

Added new files CalcMasses_atat.F, TLMT4.h to compute the
O(alpha_t^2) corrections to the Higgs boson masses of the
complex NMSSM and changed makefile accordingly.
CalcMasses.F has been changed to add the O(alpha_t^2)
corrections.
The file conversion_polemass_DRbarmass.f has been added to perform
the conversion from the top-quark pole mass to the DRbar mass taking into account
the O(alpha_t^2) contributons and to
compute the strong coupling g3 at the SUSY scale. The makefile has
been changed accordingly.
Changes in CalcMasses.F to include also the EW corrections to
the conversion of the DRbar parameters At, M_Qtilde_3 and M_tR to the
OS parameters.
Change of g3 in CalcMasses2loop.F to be in accordance with our
O(alpha_t^2) corrections.
This change only concerns the complex version: Changed the electron EDM bound to the
new experimental value of November 2018.

NMSSMCALC version 2.00, 14 January 2019

Correction of definition of absAu(3) and phiAu(3) in CalcMasses.F
(real and complex case).
Correction of definition of ReAlambdaDR and ImAlambdaDR in
CalcMasses2loop.F (real and complex case).
Correction of decay width H+ -> W H3 in bhdecay_c.f (so only complex
case): amh2 changed to amh3.

NMSSMCALC version 2.00, 28 December 2018

Added new file Zfactors_New.F and changed makefile
accordingly. With this the two-point functions and hence the
Z-factor and one-loop Higgs self-energies can be computed with
complex momentum. This gives stable results near
thresholds. (The changes are denoted by 'Nhung changed/added
...'.) In accordance with this change, there are changes in
CalcMasses.F. (They are marked by '(end)
Maggie changed 21 Aug 2018'.)
Changes in CalcMasses.F in order to account for the fact that the
SLHA2 requires Re(Alambda) and Im(Alambda) as input. (The
changes are marked by '(end) MMM changed 22 Aug 2018'.)
Changes in CalcMasses2loop.F in the computation of
Sqrt(EPhiu). It is replace by
Cos(phiu/2d0)+(0d0,1d0)*Sin(phiu/2d0) (The changes are
marked by '(end) MMM changed 21 Aug 2018'.)
These changes concern both the real and the complex NMSSMCALC programs.

NMSSMCALC version 2.00, 8 October 2017

CalcMasses2loop.F: Added 'dmtep = 0' to set the epsilon part of the
counterterm to zero in accordance with our OS renormalization scheme.

NMSSMCALC version 2.00, 22 April 2016

CalcMasses2loop.F: Bug fixed in the 2-loop counterterm for the
charged Higgs mass, when A_lambda is given as input. Thanks to Florian
Staub

NMSSMCALC version 2.00, 13 August 2015

CalcMasses.F and nmssmcalc.f have been changed as now in the complex
NMSSM also the electric dipole moments (EDMs) are calculated. The
computation is performed with the files EDMs.F and EDMloopfunc.F, which
are included in the tar file for the complex NMSSM.
bhdecay.f, bhdecay_c.f and susylha_nmssm.f changed to write out now
*all* parameters in the block EXTPAR. This is also the block that is used
in inp.dat, respectively, inp_c.dat. Furthermore, now all switches, that
have been set by the user, are given out.

NMSSMCALC version 1.03, 7 June 2015

bhdecay.f and bhdecay_c.f changed to write out now also the Higgs
coupling blocks read in by HiggsBounds.
bhdecay_c.f and bhdecay_f: If there are problems in the
Higgs-to-Higgs off-shell decays, their decay width is set to zero.
bhdecay_c.f and bhdecay.f: The total off-shell decay width has
been initialized to 0: hitohhoff(i)=0, to avoid problems when
NMSSMCALC is run in a loop over parameters.
bhdecay_c.f: Bug fix in Hi-> bb, i.e. ratcoup =
ghit(i)/cdabs(xghigb(i)) corrected to ratcoup =
ghit(i)/cdabs(xghib(i)). Thanks to Shoaib Munir.

NMSSMCALC version 1.03, 11 February 2015

Switched from the MSbar running value of alphas to the DRbar running
value.
Bug that caused wrong output values in 'Block GAUGE Q=...' fixed.
Bug that caused program to crash for At=0 fixed.
Bug that caused program to crash for M3=M_Q3 fixed.

NMSSMCALC version 1.03, 12 January 2015

Bug fixed in CalcMasses.F and CalcMasses2loop.F which occurs for
degenerate sbottom masses. Thanks to Florian Staub and Shoaib Munir
for pointing this out to us.

NMSSMCALC version 1.03, 5 January 2015

CalcMasses now version 2.00: Implementation of the two-loop O(alpha_t
alpha_s) corrections to the NMSSM
Higgs boson masses in the real and the complex NMSSM.
Due to this implementation there is a new additional file
CalcMasses2loop.F.
In inp.dat the Block MODSEL has been extended, so that the user can
choose the loop level for the calculation of the loop corrected masses and
the renormalisation scheme in the top/stop sector, i.e. DR_bar or OS.
The implementation of the two-loop corrections furthermore lead to
small changes in the files susylha_nmssm.f and in the Makefile.
Bug fix in hsqsq_nmssmcalc.f: xlamvs/sqrt(2) changed to xlamvs.
Changes in bhdecay.f and bhdecay_c.f so that warnings are given out in
case there are problems in the off-shell decays.
nmssmcalc.f: Changes in print statements on the screen.

NMSSMCALC version 1.02, 24 August 2014

nmssmcalc.f: Changes in print statements on the screen.

Changes for NMSSMCALC version 1.02, 30 June 2014

CalcMasses.F: in subroutine neutralinomasses: fixed a bug in an
if-statement.
CalcMasses.F: changed interpretation of a negative mu-value for the
real limit.
bhdecay.f: fixed a definition in the common block of
elemhtt_NMSSM.

NMSSMCALC version 1.02, 11 April 2014

Off-shell Higgs-to-Higgs decays stabilized in real and complex NMSSM.
(Thanks to S. Munir)
Warning is issued if QCD corrections to Higgs decays into
stops/sbottoms become too large in some extreme scenarios. In this case
the tree-level decay widths are calculated. This is only relevant for
the real NMSSM. (Thanks to S. Munir)
tan(beta) is now consistently read in from EXTPAR. (Thanks to the CPC
referee.)
Consistent use of VEV.

NMSSMCALC version 1.01, based now on HDECAY 6.10, 22 January 2014

Typo corrected in a do loop in CalcMasses.F (thanks to S. Munir)
If the mass calculation does not converge sufficiently fast, a
warning is issued and the last masses are given out by CalcMasses.F.
Changes in bhdecay.f: signs in Delta_tau/mu corrected, elw.
corrections extended, elw. corrections to Delta_b/s added; pdg code for H1
in output file corrected (thanks to S. Munir); warning included in case of
problems in off-shell decays.
Changes in bhdecay_c.f: signs in Delta_tau/mu corrected; pdg code for
H1 in output file corrected (thanks to S. Munir); warning included in
case of problems in off-shell decays.