X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Fcluster%2Fwham%2Fsrc%2Freadrtns.F;h=ac18e17b8758fddd870a6475a7631aef9b632039;hb=2605ef637793380790ce7abbae21ef9d88f0e505;hp=037eb1fce5361e02bf9730bb4ec3394d6b57d7be;hpb=3d6f9e7811a3341955e6fdda721ab6de34efe252;p=unres.git

diff --git a/source/cluster/wham/src/readrtns.F b/source/cluster/wham/src/readrtns.F
index 037eb1f..ac18e17 100644
--- a/source/cluster/wham/src/readrtns.F
+++ b/source/cluster/wham/src/readrtns.F
@@ -275,7 +275,7 @@ C Convert sequence to numeric code
       do i=1,nres
         itype(i)=rescode(i,sequence(i),iscode)
       enddo
-      if (itype(2).eq.10) then
+      if (itype(2).eq.10.and.itype(1).eq.ntyp1) then
         write (iout,*)
      &   "Glycine is the first full residue, initial dummy deleted"
         do i=1,nres
@@ -283,7 +283,7 @@ C Convert sequence to numeric code
         enddo
         nres=nres-1
       endif
-      if (itype(nres).eq.10) then
+      if (itype(nres-1).eq.10.and.itype(nres).eq.ntyp1) then
         write (iout,*)
      &   "Glycine is the last full residue, terminal dummy deleted"
         nres=nres-1
@@ -1190,8 +1190,8 @@ c    &                      "rescore(",k,i-1,") =",rescore(k,i-1),
 c    &                      "rescore(",k,i-2,") =",rescore(k,i-2),
 c    &                      "rescore(",k,i-3,") =",rescore(k,i-3)
 
-            sigma_dih(k,i)=rescore(k,i)+rescore(k,i-1)+
-     &                     rescore(k,i-2)+rescore(k,i-3)  !  right expression ?
+            sigma_dih(k,i)=(rescore(k,i)+rescore(k,i-1)+
+     &                     rescore(k,i-2)+rescore(k,i-3))/4.0  !  right expression ?
 c
 c           write (iout,*) "Raw sigmas for dihedral angle restraints"
 c           write (iout,'(i5,10(2f8.2,4x))') i,sigma_dih(k,i)
@@ -1224,8 +1224,8 @@ c            write(iout,*)  "rescore(",k,i,") =",rescore(k,i),
 c    &                      "rescore(",k,i-1,") =",rescore(k,i-1),
 c    &                      "rescore(",k,i-2,") =",rescore(k,i-2)
 c            read (ientin,*) sigma_theta(k,i) ! 1st variant
-             sigma_theta(k,i)=rescore(k,i)+rescore(k,i-1)+
-     &                        rescore(k,i-2) !  right expression ?
+             sigma_theta(k,i)=(rescore(k,i)+rescore(k,i-1)+
+     &                        rescore(k,i-2))/3.0 !  right expression ?
              sigma_theta(k,i)=1.0d0/(sigma_theta(k,i)*sigma_theta(k,i))
 
 c            sigma_theta(k,i)=hmscore(k)*rescore(k,i)*rescore(k,i-1)*