call MPI_Barrier(FG_COMM,IERR)
call MPI_Reduce(e_tmp,ehomology_constr,1,
& MPI_DOUBLE_PRECISION,MPI_SUM,king,FG_COMM,IERR)
+ else if (iorder.eq.13) then
+ call chainbuild
endif
enddo
write (*,*) 'Processor',fg_rank,' CG group',kolor,
f=energia(0)
c if (jjj.gt.0) then
c write (iout,'(10f8.3)') (rad2deg*x(i),i=1,n)
-c write (iout,*) 'f=',etot
+c write (iout,*) 'f=',f
c jjj=0
c endif
return