etota = energy_(0)
etot = etota
call enerprint(energy_)
+ if (.false.) then
call hairpin(.true.,nharp,iharp)
call secondary2(.true.)
+ endif
if (minim) then
!rc overlap test
+ print *,"overlap",searchsc,overlapsc
if (overlapsc) then
print *, 'Calling OVERLAP_SC'
call overlap_sc(fail)