X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Fwham%2Fsrc%2Fenecalc1.F;h=3d878ad9547103806345d80ad3d9374442f86630;hb=6304fe426dd266ba99503f2cc4a5b9d396511b2f;hp=bdda1411dd1ba0968b3d31ea853e283b9e3131ff;hpb=90a47ba2f9ea5c266c59b225f8505a81d75420d7;p=unres.git diff --git a/source/wham/src/enecalc1.F b/source/wham/src/enecalc1.F index bdda141..3d878ad 100644 --- a/source/wham/src/enecalc1.F +++ b/source/wham/src/enecalc1.F @@ -167,18 +167,20 @@ c & " kfac",kfac,"quot",quot," fT",fT call etotal(energia(0),fT) if (constr_homology) energia(0)=energia(0)+ & waga_homology(iset)*energia(22) + write (iout,*) "constr_homology",constr_homology," iset",iset, + & " waga_homology",waga_homology(iset) #ifdef DEBUG write (iout,*) "Conformation",i call enerprint(energia(0),fT) c write (iout,'(2i5,21f8.2)') i,iparm,(energia(k),k=1,21) c write (iout,*) "ftors",ftors - write (iout,'(a,3(3x,a5),5x,3(3x,a5))') - & "Residue","X(CA)","Y(CA)","Z(CA)","X(SC)","Y(SC)","Z(SC)" - do ires=1,nres - write (iout,'(a3,1x,i3,3f8.3,5x,3f8.3)') - & restyp(itype(ires)),ires,(c(j,ires),j=1,3), - & (c(j,ires+nres),j=1,3) - enddo +c write (iout,'(a,3(3x,a5),5x,3(3x,a5))') +c & "Residue","X(CA)","Y(CA)","Z(CA)","X(SC)","Y(SC)","Z(SC)" +c do ires=1,nres +c write (iout,'(a3,1x,i3,3f8.3,5x,3f8.3)') +c & restyp(itype(ires)),ires,(c(j,ires),j=1,3), +c & (c(j,ires+nres),j=1,3) +c enddo call intout #endif if (energia(0).ge.1.0d20) then @@ -519,7 +521,7 @@ c else c & potE(i,iparm),-entfac(i),rms_nat,iscore & potE(i,nparmset),-entfac(i),rms_nat,iscore c endif - write (iout,'(2i5,3e15.5)') i,me,potE(i,iparm),-entfac(i) +c write (iout,'(2i5,3e15.5)') i,me,potE(i,iparm),-entfac(i) #ifndef MPI if (cxfile) call cxwrite(ixdrf,csingle,potE(i,nparmset), & -entfac(i),rms_nat,iscore)