box shift in wham from Adam
[unres.git] / source / wham / src-M / cxread.F
index c0319f6..15777b1 100644 (file)
@@ -168,6 +168,20 @@ c      call flush(iout)
         enddo
       enddo
 
+c Box shift  
+      call oligomer
+      do i=1,nres
+        do j=1,3
+          xoord(j,i)=c(j,i)
+        enddo
+      enddo  
+      do i=1,nct-nnt+1
+        do j=1,3
+          xoord(j,i+nres)=c(j,i+nres+nnt-1)
+        enddo
+      enddo  
+c end change 
+
       if (islice.gt.0 .and. islice.le.nslice .and. (.not.separate_parset
      &    .or. iset.eq.myparm)) then
         ii=ii+1