+ else if (iorder.eq.11) then
+c broadcast new value of iset
+ call MPI_Bcast(iset,1,MPI_INTEGER,king,FG_COMM,IERROR)
+ else if (iorder.eq.12) then
+c broadcast new value of iset and
+c calculate restraint homology energy and
+c sum it over FG_COMM
+ call MPI_Bcast(iset,1,MPI_INTEGER,king,FG_COMM,IERROR)
+ call e_modeller(e_tmp)
+ 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