X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Funres%2Fsrc_MD%2Fint_to_cart.f;h=272d39b33245d4166b4698fd7b1a78ec7bb38200;hb=5316f1d82cb17d390be734cefbfca809cc65dcc4;hp=10512187467977b722deb511a955d7e1e96a54cb;hpb=3dbe5ceeea4b858bc25fa1469237e697c0bf293f;p=unres.git diff --git a/source/unres/src_MD/int_to_cart.f b/source/unres/src_MD/int_to_cart.f index 1051218..272d39b 100644 --- a/source/unres/src_MD/int_to_cart.f +++ b/source/unres/src_MD/int_to_cart.f @@ -162,7 +162,8 @@ c Calculating the remainder of dE/ddc2 if((itype(2).ne.10).and.(itype(2).ne.ntyp1)) then if (itype(1).ne.10) gxcart(j,2)=gxcart(j,2)+ & gloc_sc(3,0,icg)*dtauangle(j,3,3,3) - if ((itype(3).ne.10).and.(nres.ge.3).and.(itype(3).ne.ntyp1)) then + if ((itype(3).ne.10).and.(nres.ge.3).and.(itype(3).ne.ntyp1)) + & then gxcart(j,2)=gxcart(j,2)-gloc_sc(3,1,icg)*dtauangle(j,3,1,4) cc the - above is due to different vector direction gcart(j,2)=gcart(j,2)+gloc_sc(3,1,icg)*dtauangle(j,3,2,4)