X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Funres%2Fsrc_MD%2Freadpdb.F;h=eed985e003b6a5f01fe4f5a560eee891a95d246e;hb=27d0898cd1f4b898e6b044f96d7f28615e23fadb;hp=745b256c42680305a09a501387cc7f0ecbcf8f8d;hpb=4764dff6d4a43d5ae03320fbd0d2dbef3bb37e3f;p=unres.git

diff --git a/source/unres/src_MD/readpdb.F b/source/unres/src_MD/readpdb.F
index 745b256..eed985e 100644
--- a/source/unres/src_MD/readpdb.F
+++ b/source/unres/src_MD/readpdb.F
@@ -165,11 +165,24 @@ C Calculate the CM of the last side chain.
       if (itype(nres).ne.10) then
         nres=nres+1
         itype(nres)=21
+        if (unres_pdb) then
+C 2/15/2013 by Adam: corrected insertion of the last dummy residue
+          call refsys(nres-3,nres-2,nres-1,e1,e2,e3,fail)
+          if (fail) then
+            e2(1)=0.0d0
+            e2(2)=1.0d0
+            e2(3)=0.0d0
+          endif
+          do j=1,3
+            c(j,nres)=c(j,nres-1)-3.8d0*e2(j)
+          enddo
+        else
         do j=1,3
           dcj=c(j,nres-2)-c(j,nres-3)
           c(j,nres)=c(j,nres-1)+dcj
           c(j,2*nres)=c(j,nres)
         enddo
+        endif
       endif
       do i=2,nres-1
         do j=1,3
@@ -729,6 +742,7 @@ C Copy the coordinates to reference coordinates
       do i=1,2*nres
         do j=1,3
           cref(j,i)=c(j,i)
+          chomo(j,i,k)=c(j,i)
         enddo
       enddo