X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Fcluster%2Fwham%2Fsrc%2Fenergy_p_new.F;h=8006a41109e4b94193e6d85c0fdbe48987a71ca0;hb=f37cfff89ef75a184b0c7bbbac2ecea2efc1d157;hp=b3b3628fe02a11e18106ccc2951f6c9aaa5240bf;hpb=c358637b84ca9f64f98d5205fe35dfbc916cbe61;p=unres.git diff --git a/source/cluster/wham/src/energy_p_new.F b/source/cluster/wham/src/energy_p_new.F index b3b3628..8006a41 100644 --- a/source/cluster/wham/src/energy_p_new.F +++ b/source/cluster/wham/src/energy_p_new.F @@ -3345,11 +3345,7 @@ c write (iout,*) idihconstr_start_homo,idihconstr_end_homo enddo #endif do i=idihconstr_start_homo,idihconstr_end_homo -#ifdef OLDRESTR kat2=0.0d0 -#else - kat2=nexl -#endif c betai=beta(i,i+1,i+2,i+3) betai = phi(i) c write (iout,*) "betai =",betai @@ -3464,11 +3460,7 @@ c c Deviation of theta angles wrt constr_homology ref structures c utheta_i=0.0d0 ! argument of Gaussian for single k -#ifdef OLDRESTR gutheta_i=0.0d0 ! Sum of Gaussians over constr_homology ref structures -#else - gutheta_i=nexl -#endif c do j=ifrag_back(1,i,iset)+2,ifrag_back(2,i,iset) ! original loop c over residues in a fragment c write (iout,*) "theta(",i,")=",theta(i) @@ -3536,11 +3528,7 @@ c write (iout,*) "waga_d",waga_d #endif do i=loc_start,loc_end usc_diff_i=0.0d0 ! argument of Gaussian for single k -#ifdef OLDRESTR guscdiff(i)=0.0d0 ! Sum of Gaussians over constr_homology ref structures -#else - guscdiff(i)=nexl -#endif c do j=ifrag_back(1,i,iset)+1,ifrag_back(2,i,iset)-1 ! Econstr_back legacy c write(iout,*) "xxtab, yytab, zztab" c write(iout,'(i5,3f8.2)') i,xxtab(i),yytab(i),zztab(i)