if (constr_homology.ge.1) then
call e_modeller(ehomology_constr)
- print *,'iset=',iset,'me=',me,ehomology_constr,
- & 'Processor',fg_rank,' CG group',kolor,
- & ' absolute rank',MyRank
+c print *,'iset=',iset,'me=',me,ehomology_constr,
+c & 'Processor',fg_rank,' CG group',kolor,
+c & ' absolute rank',MyRank
else
ehomology_constr=0.0d0
endif
c
c ehomology_constr=odleg+kat
if (homol_nset.gt.1)then
- ehomology_constr=waga_dist1(iset)*odleg+waga_angle1(iset)*kat+waga_theta*Eval
+ ehomology_constr=waga_dist1(iset)*odleg+waga_angle1(iset)*kat
+ & +waga_theta*Eval
& +waga_d*Erot
else
ehomology_constr=waga_dist*odleg+waga_angle*kat+waga_theta*Eval