X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;ds=sidebyside;f=source%2Funres%2Fsrc_MIN%2Freadpdb.F;h=08bdf960842fed8bb8b02205c65a366dfb7c910f;hb=5a908648c2ff031dd2300a4250a151bad233f898;hp=eb4ba3f5fddce94692a64359763fd292c7768ff7;hpb=34b56f82298df0d9781948908b2d93415838d88e;p=unres.git

diff --git a/source/unres/src_MIN/readpdb.F b/source/unres/src_MIN/readpdb.F
index eb4ba3f..08bdf96 100644
--- a/source/unres/src_MIN/readpdb.F
+++ b/source/unres/src_MIN/readpdb.F
@@ -250,7 +250,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