if (outmol2) call mol2out(etot,titel(:32))
if (refstr) call rms_nac_nnc(rms,frac,frac_nn,co,.true.)
write (iout,'(a,i3)') 'SUMSL return code:',iretcode
if (outmol2) call mol2out(etot,titel(:32))
if (refstr) call rms_nac_nnc(rms,frac,frac_nn,co,.true.)
write (iout,'(a,i3)') 'SUMSL return code:',iretcode