X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Fcluster%2Fwham%2Fsrc-HCD-5D%2Fchain_symmetry.F;fp=source%2Fcluster%2Fwham%2Fsrc-HCD-5D%2Fchain_symmetry.F;h=8c368555d6cd6e3f59b978b7c2ed55341efd008b;hb=6260e60062f93c75d71ad18219a28f1ab5611dec;hp=1406d1d20f84e514ac65e88550481c5d1822277d;hpb=7baaf36efbc265e029ec93a4bb1adfadfd55d7bf;p=unres.git diff --git a/source/cluster/wham/src-HCD-5D/chain_symmetry.F b/source/cluster/wham/src-HCD-5D/chain_symmetry.F index 1406d1d..8c36855 100644 --- a/source/cluster/wham/src-HCD-5D/chain_symmetry.F +++ b/source/cluster/wham/src-HCD-5D/chain_symmetry.F @@ -7,6 +7,7 @@ c implicit none include "DIMENSIONS" include "COMMON.IOUNITS" + include "COMMON.CONTROL" integer nchain,nres,itype(nres),chain_border(2,maxchain), & chain_length(nchain),itemp(maxchain), & npermchain,tabpermchain(maxchain,maxperm), @@ -42,6 +43,7 @@ c nchain_group=nchain_group+1 iieq=1 iequiv(iieq,nchain_group)=i + if (symetr.eq.1) then do j=i+1,nchain if (iflag(j).gt.0.or.chain_length(i).ne.chain_length(j)) cycle c k=0 @@ -57,6 +59,7 @@ c k=k+1 iieq=iieq+1 iequiv(iieq,nchain_group)=j enddo + endif nequiv(nchain_group)=iieq enddo write(iout,*) "Number of equivalent chain groups:",nchain_group