projects
/
unres.git
/ blobdiff
commit
grep
author
committer
pickaxe
?
search:
re
summary
|
shortlog
|
log
|
commit
|
commitdiff
|
tree
raw
|
inline
| side by side
cluster wham test homo small memory single sco
[unres.git]
/
source
/
wham
/
src-M
/
enecalc1.F
diff --git
a/source/wham/src-M/enecalc1.F
b/source/wham/src-M/enecalc1.F
index
ab0a408
..
697c31d
100644
(file)
--- a/
source/wham/src-M/enecalc1.F
+++ b/
source/wham/src-M/enecalc1.F
@@
-91,6
+91,11
@@
q(nQ+1,iii+1)=rmsnat(iii+1)
endif
q(nQ+2,iii+1)=gyrate(iii+1)
q(nQ+1,iii+1)=rmsnat(iii+1)
endif
q(nQ+2,iii+1)=gyrate(iii+1)
+ if (q(nQ+1,iii+1).le.frac_cutoff) then
+ q(nQ+3,iii+1)=1.0d0
+ else
+ q(nQ+3,iii+1)=0.0d0
+ endif
c fT=T0*beta_h(ib,ipar)*1.987D-3
c ft=2.0d0/(1.0d0+1.0d0/(T0*beta_h(ib,ipar)*1.987D-3))
if (rescale_mode.eq.1) then
c fT=T0*beta_h(ib,ipar)*1.987D-3
c ft=2.0d0/(1.0d0+1.0d0/(T0*beta_h(ib,ipar)*1.987D-3))
if (rescale_mode.eq.1) then
@@
-202,6
+207,9
@@
c call pdbout(ii+1,beta_h(ib,ipar),efree,energia(0),0.0d0,rmsdev)
& iii+1,indstart(me1)+iii," T",
& 1.0d0/(1.987D-3*beta_h(ib,ipar))
c call intout
& iii+1,indstart(me1)+iii," T",
& 1.0d0/(1.987D-3*beta_h(ib,ipar))
c call intout
+C write (iout,*) "The Cartesian geometry is:"
+C write (iout,'(8f10.5)') ((c(l,k),l=1,3),k=1,nres)
+C write (iout,'(8f10.5)') ((c(l,k+nres),l=1,3),k=nnt,nct)
call pdbout(indstart(me1)+iii,
& 1.0d0/(1.987D-3*beta_h(ib,ipar)),
&energia(0),eini,0.0d0,0.0d0)
call pdbout(indstart(me1)+iii,
& 1.0d0/(1.987D-3*beta_h(ib,ipar)),
&energia(0),eini,0.0d0,0.0d0)