NMHDECAY: A Fortran Code for the Higgs Masses, Couplings and Decay   Widths in the NMSSM
  The Fortran code NMHDECAY computes the masses, couplings and decay widths of all Higgs bosons of the NMSSM in terms of its parameters at the electroweak (SUSY breaking) scale: the Yukawa couplings lambda and kappa, the soft trilinear terms A_lambda and A_kappa, and tan(beta) and mu_eff = lambda*<S>. The computation of the spectrum includes leading two loop terms, electroweak corrections and propagator corrections. The computation of the decay widths is carried out as in HDECAY, but (for the moment) without three body decays. Each point in parameter space is checked against negative Higgs bosons searches at LEP, including unconventional channels relevant for the NMSSM. One version of the program uses generalized SLHA conventions for input and output. 
Ulrich Ellwanger, John F. Gunion, Cyril Hugonie