do i=1,nside_move
ii=ind_side(i)
iti=itype(ii,1)
- call gen_side(iti,theta(ii+1),alph(ii),omeg(ii),fail)
+ call gen_side(iti,theta(ii+1),alph(ii),omeg(ii),fail,1)
if (fail) then
isctry=isctry+1
if (isctry.gt.1000) then
! THETA move
400 end_select=iran_num(3,nres)
theta_new=gen_theta(itype(end_select,1),phi(end_select),&
- phi(end_select+1))
+ phi(end_select+1),1)
if (print_mc.gt.1) write (iout,'(a,i3,a,f8.3,a,f8.3)') &
'Theta ',end_select,' moved from',theta(end_select)*rad2deg,&
' to ',theta_new*rad2deg