projects
/
unres.git
/ blobdiff
commit
grep
author
committer
pickaxe
?
search:
re
summary
|
shortlog
|
log
|
commit
|
commitdiff
|
tree
raw
|
inline
| side by side
multichain cluster can select number of clusters instead of cutoff
[unres.git]
/
source
/
cluster
/
wham
/
src-M
/
readrtns.F
diff --git
a/source/cluster/wham/src-M/readrtns.F
b/source/cluster/wham/src-M/readrtns.F
index
2573b96
..
a3feaa5
100644
(file)
--- a/
source/cluster/wham/src-M/readrtns.F
+++ b/
source/cluster/wham/src-M/readrtns.F
@@
-91,7
+91,8
@@
C long axis of side chain
min_var=(index(controlcard,'MINVAR').gt.0)
plot_tree=(index(controlcard,'PLOT_TREE').gt.0)
punch_dist=(index(controlcard,'PUNCH_DIST').gt.0)
min_var=(index(controlcard,'MINVAR').gt.0)
plot_tree=(index(controlcard,'PLOT_TREE').gt.0)
punch_dist=(index(controlcard,'PUNCH_DIST').gt.0)
- call readi(controlcard,'NCUT',ncut,1)
+ call readi(controlcard,'NCUT',ncut,0)
+ call readi(controlcard,'NCLUST',nclust,5)
call readi(controlcard,'SYM',symetr,1)
write (iout,*) 'sym', symetr
call readi(controlcard,'NSTART',nstart,0)
call readi(controlcard,'SYM',symetr,1)
write (iout,*) 'sym', symetr
call readi(controlcard,'NSTART',nstart,0)
@@
-102,7
+103,8
@@
C long axis of side chain
lgrp=(index(controlcard,'LGRP').gt.0)
caonly=(index(controlcard,'CA_ONLY').gt.0)
print_dist=(index(controlcard,'PRINT_DIST').gt.0)
lgrp=(index(controlcard,'LGRP').gt.0)
caonly=(index(controlcard,'CA_ONLY').gt.0)
print_dist=(index(controlcard,'PRINT_DIST').gt.0)
- call multreada(controlcard,'CUTOFF',rcutoff,ncut,-1.0d0)
+ if (ncut.gt.0)
+ & call multreada(controlcard,'CUTOFF',rcutoff,ncut,-1.0d0)
call readi(controlcard,'IOPT',iopt,2)
lside = index(controlcard,"SIDE").gt.0
efree = index(controlcard,"EFREE").gt.0
call readi(controlcard,'IOPT',iopt,2)
lside = index(controlcard,"SIDE").gt.0
efree = index(controlcard,"EFREE").gt.0
@@
-1103,7
+1105,7
@@
c
tpl_k_rescore="template"//kic2//".sco"
unres_pdb=.false.
tpl_k_rescore="template"//kic2//".sco"
unres_pdb=.false.
- call readpdb
+ call readpdb_template(k)
do i=1,2*nres
do j=1,3
crefjlee(j,i)=c(j,i)
do i=1,2*nres
do j=1,3
crefjlee(j,i)=c(j,i)