etot=v(10)
iretcode=iv(1)
cd print *,'Exit SUMSL; return code:',iretcode,' energy:',etot
- write (iout,'(/a,i4,a,f15.5/)') 'SUMSL return code:',iv(1),
- & " energy",etot
+cd write (iout,'(/a,i4,a,f15.5/)') 'SUMSL return code:',iv(1),
+cd & " energy",etot
call var_to_geom(nvar,x)
c call chainbuild
c call etotal(energia(0))