X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Funres%2Fsrc_MD%2Freadpdb.F;h=80c3fa375d64853f21d99a7a335cad0305f5cf2e;hb=4e96aa60029501bc6bffb1b3b9f9fa5773295929;hp=48e0abdfbd8b810c536afda590bff09b60c2a8be;hpb=6b8a6eb9eb98a6a74a90cbe8da1f254349c9ca4e;p=unres.git

diff --git a/source/unres/src_MD/readpdb.F b/source/unres/src_MD/readpdb.F
index 48e0abd..80c3fa3 100644
--- a/source/unres/src_MD/readpdb.F
+++ b/source/unres/src_MD/readpdb.F
@@ -252,7 +252,20 @@ ctest          stop
         vbld(i+1)=dist(i,i+1)
         vbld_inv(i+1)=1.0d0/vbld(i+1)
         if (i.gt.1) theta(i+1)=alpha(i-1,i,i+1)
-        if (i.gt.2) phi(i+1)=beta(i-2,i-1,i,i+1)
+        if (i.gt.2) then
+        if (i.le.nres) phi(i+1)=beta(i-2,i-1,i,i+1)
+        if ((itype(i).ne.10).and.(itype(i-1).ne.10)) then
+         tauangle(3,i+1)=beta(i+nres-1,i-1,i,i+nres)
+        endif
+        if (itype(i-1).ne.10) then
+         tauangle(1,i+1)=beta(i-1+nres,i-1,i,i+1)
+         omicron(1,i)=alpha(i-2,i-1,i-1+nres)
+         omicron(2,i)=alpha(i-1+nres,i-1,i)
+        endif
+        if (itype(i).ne.10) then
+         tauangle(2,i+1)=beta(i-2,i-1,i,i+nres)
+        endif
+        endif
       enddo
 c      if (unres_pdb) then
 c        if (itype(1).eq.21) then