X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Funres%2Fsrc-HCD-5D%2Frefsys.f;fp=source%2Funres%2Fsrc-HCD-5D%2Frefsys.f;h=5573f0578bbcaf0d73b969ad6aca16776200afc2;hb=58980cd5a21077fd523753ffccc036765ef70d82;hp=4b7b76379c0d869f6cdb200e01cba831d8ef9e78;hpb=d02292725c202ff9c2749beac934bf1630f9017e;p=unres.git diff --git a/source/unres/src-HCD-5D/refsys.f b/source/unres/src-HCD-5D/refsys.f index 4b7b763..5573f05 100644 --- a/source/unres/src-HCD-5D/refsys.f +++ b/source/unres/src-HCD-5D/refsys.f @@ -21,12 +21,12 @@ c print *,'just initialize' c print *,fail s1=0.0 s2=0.0 - print *,s1,s2 +c print *,s1,s2 do 1 i=1,3 - print *, i2,i3,i4 +c print *, i2,i3,i4 zi=c(i,i2)-c(i,i3) ui=c(i,i4)-c(i,i3) - print *,zi,ui +c print *,zi,ui s1=s1+zi*zi s2=s2+ui*ui z(i)=zi @@ -41,7 +41,7 @@ c print *,fail write(iout,1000) i3,i4,i1 fail=.true. return - print *,'two if pass' +c print *,'two if pass' 4 s1=1.0/s1 s2=1.0/s2 v1=z(2)*u(3)-z(3)*u(2)