DFA & lipid

[unres.git] / source / wham / src-HCD-5D / cxread.F

diff --git a/source/wham/src-HCD-5D/cxread.F b/source/wham/src-HCD-5D/cxread.F
index e3e5fcb..36ef6e6 100644 (file)
--- a/source/wham/src-HCD-5D/cxread.F
+++ b/source/wham/src-HCD-5D/cxread.F
@@ -171,8 +171,12 @@ c      call flush(iout)
           c(j,i+nres+nnt-1)=xoord(j,i+nres)
         enddo
       enddo
+c      write (iout,*) "Before boxshift"
+c      call flush(iout)
 c Box shift
       call oligomer
+c      write (iout,*) "After oligomer"
+c      call flush(iout)
       do i=1,nres
         do j=1,3
           xoord(j,i)=c(j,i)
@@ -184,7 +188,8 @@ c Box shift
         enddo
       enddo
 c end change
-
+c      write (iout,*) "Before islice"
+c      call flush(iout)
       if (islice.gt.0 .and. islice.le.nslice .and. (.not.separate_parset
      &    .or. iset.eq.myparm)) then
         ii=ii+1
@@ -300,6 +305,7 @@ c     &          bprotfile_temp(:ilen(bprotfile_temp))
           write (iout,'(8f10.5)') ((c(j,i),j=1,3),i=1,nres)
           write (iout,'(8f10.5)') ((c(j,i+nres),j=1,3),i=nnt,nct)
           write (iout,*) "Internal coordinates"
+          call intout
           write (iout,'(8f10.4)') (vbld(k),k=nnt+1,nct)
           write (iout,'(8f10.4)') (vbld(k),k=nres+nnt,nres+nct)
           write (iout,'(8f10.4)') (rad2deg*theta(k),k=3,nres)