We present the program package NMSSMCALC for the calculation of the loop-corrected NMSSM Higgs boson masses and decay widths in the CP-conserving and CP-violating NMSSM. The full one-loop corrections to the Higgs boson masses are evaluated in a mixed renormalisation scheme of on-shell and DR¯ conditions. The Higgs decay widths include the dominant higher order QCD corrections, and the decays into bottom quarks, strange quarks and τ leptons are supplemented by higher order SUSY corrections through effective couplings. All relevant off-shell decays into two massive gauge bosons, gauge and Higgs boson and Higgs pair final states as well as into heavy quark pairs are computed. The input and output files feature the SUSY Les Houches Accord so that the program can easily be linked with existing computer tools. Program summaryProgram title: NMSSMCALCCatalogue identifier: AEUE_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEUE_v1_0.htmlProgram obtainable from: CPC Program Library, Queen’s University, Belfast, N. IrelandLicensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.htmlNo. of lines in distributed program, including test data, etc.: 53070No. of bytes in distributed program, including test data, etc.: 329132Distribution format: tar.gzProgramming language: Fortran77.Computer: Any with a Fortran77 system.Operating system: Linux, Unix.RAM: 0.5MBClassification: 11.1.Nature of problem:Computation of the NMSSM Higgs mass spectrum including higher order corrections in a mixed renormalisation scheme of on-shell and DRbar conditions, and numerical calculation of the decay widths and branching ratios, both in the CP-conserving and in the CP-violating NMSSM. The decay widths include the dominant higher order QCD corrections and, for the neutral Higgs boson decays into a bottom quark pair, the higher order SUSY–QCD and the approximate SUSY–electroweak (elw) corrections up to one-loop accuracy. The decays into a strange quark pair include the dominant resummed SUSY–QCD corrections and the one into a τ pair the dominant resummed SUSY–elw corrections. Analogously for the charged Higgs boson the higher order SUSY corrections have been implemented for the decays into fermion pairs. In the real NMSSM, the decays into stop and sbottom pairs, respectively, contain the SUSY–QCD corrections. All relevant off-shell decays into massive gauge bosons, into gauge and Higgs bosons, into Higgs pairs and into heavy quark pairs have been taken into account. The input and output files feature the SUSY Les Houches Accord (SLHA).Solution method:The necessary input values are set in the two input files inp.dat (in the SLHA format) and bhdecay.in. The file inp.dat, where the choice between the real and the complex NMSSM can be made, is read in by CalcMasses.F, which calculates the one-loop corrected NMSSM Higgs mass spectrum. CalcMasses.F writes out all necessary parameters, masses and mixing angles in an SLHA format file slha.in, which is read in by bhdecay.f (by bhdecay c.f in the complex case). This Fortran routine then computes the decay widths and branching ratios and writes them out in an SLHA format file slha decay.out. Furthermore it writes out all parameters and mixing angles.Restrictions:At present the NMSSM Higgs spectrum is calculated at one-loop accuracy, the Higgs self-couplings at leading order, and no renormalisation group running of the input parameters is included. However, we are about to implement the two-loop corrections to the NMSSM Higgs boson masses and the Higgs self-couplings. Furthermore, the renormalisation group equations for the input parameters shall be provided in the near future. The program does not provide any distributions.Running time:Less than one second per point