call enerprint(energy(0))
if (secondary_str) then
call hairpin(.true.,nharp,iharp)
-c print *,'after hairpin'
+ print *,'after hairpin'
call secondary2(.true.)
endif
-c print *,'after secondary'
+ print *,'after secondary'
if (minim) then
crc overlap test
if (indpdb.ne.0 .and. .not.dccart) then