do i=nnt,nct
if (itype(i,1).ne.10) then
!d print *,'i=',i,' itype=',itype(i,1),' theta=',theta(i+1)
- call gen_side(itype(i,1),theta(i+1),alph(i),omeg(i),fail)
+ call gen_side(itype(i,1),theta(i+1),alph(i),omeg(i),fail,1)
endif
enddo
call chainbuild
alph0=alph(ind_sc)
omeg0=omeg(ind_sc)
call gen_side(itype(ind_sc,1),theta(ind_sc+1),alph(ind_sc),&
- omeg(ind_sc),fail)
+ omeg(ind_sc),fail,1)
call chainbuild
call etotal(energia)
!d write (iout,'(a,i5,a,i4,2(a,f8.3),2(a,1pe14.5))')