X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Fwham%2Fsrc-M%2Fenecalc1.F;h=697c31da0fa8edc91a9d48fd5f774f4c7e13e633;hb=7baaf36efbc265e029ec93a4bb1adfadfd55d7bf;hp=913c2e103793babe6b1851bde474dfd204b3adfb;hpb=108d195fed23f7ef940a4f3188c7f0e2790e5f17;p=unres.git

diff --git a/source/wham/src-M/enecalc1.F b/source/wham/src-M/enecalc1.F
index 913c2e1..697c31d 100644
--- 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)
+         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
@@ -202,9 +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
-          write (iout,*) "The Cartesian geometry is:"
-          write (iout,'(8f10.5)') ((c(l,k),l=1,3),k=1,nres)
-          write (iout,'(8f10.5)') ((c(l,k+nres),l=1,3),k=nnt,nct)
+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)