X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Fcluster%2Fwham%2Fsrc-M%2Freadrtns.F;h=60b1ec2d29200f5b5369c9b0098f2d1aafc2a373;hb=c0b48650276909d82678d8b4263049fd1ae77887;hp=e5c31cf33efc9a894ed18efcc2407c6e125bfe2a;hpb=a1e2f7505712e67c4a485fc9b56e8cdfdc25c84e;p=unres.git

diff --git a/source/cluster/wham/src-M/readrtns.F b/source/cluster/wham/src-M/readrtns.F
index e5c31cf..60b1ec2 100644
--- a/source/cluster/wham/src-M/readrtns.F
+++ b/source/cluster/wham/src-M/readrtns.F
@@ -131,8 +131,9 @@ C long axis of side chain
       call readi(controlcard,'NSAXS',nsaxs,0)
       call readi(controlcard,'SAXS_MODE',saxs_mode,0)
       call reada(controlcard,'SCAL_RAD',scal_rad,1.0d0)
+      call reada(controlcard,'SAXS_CUTOFF',saxs_cutoff,1.0d0)
       write (iout,*) "Number of SAXS restraints",NSAXS," SAXS_MODE",
-     &   SAXS_MODE," SCAL_RAD",scal_rad
+     &   SAXS_MODE," SCAL_RAD",scal_rad,"SAXS_CUTOFF",saxs_cutoff
       if (min_var) iopt=1
       return
       end
@@ -997,6 +998,11 @@ c SAXS distance distribution
       do i=1,nsaxs
         write (iout,'(f8.2,e15.5)') distsaxs(i),Psaxs(i)
       enddo
+      Wsaxs0=0.0d0
+      do i=1,nsaxs
+        Wsaxs0=Wsaxs0-Psaxs(i)*dlog(Psaxs(i))
+      enddo
+      write (iout,*) "Wsaxs0",Wsaxs0
       else
 c SAXS "spheres".
       do i=1,nsaxs