X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Funres%2Fsrc-HCD-5D%2Fgen_rand_conf.F;h=557f435a3ebbfd0301e1399ae3af10fa891f7b1c;hb=0b41c0eed4b401e3554fac2830d4234f83abdd98;hp=9ab248077b0a338a109bb17e049471e0898c809e;hpb=48ae9e01d2dd6571fa2cca6c704dc04f86e5fd7b;p=unres.git

diff --git a/source/unres/src-HCD-5D/gen_rand_conf.F b/source/unres/src-HCD-5D/gen_rand_conf.F
index 9ab2480..557f435 100644
--- a/source/unres/src-HCD-5D/gen_rand_conf.F
+++ b/source/unres/src-HCD-5D/gen_rand_conf.F
@@ -804,8 +804,8 @@ c     overlapping residues left, or false otherwise (success)
 
         call chainbuild_extconf
         call overlap_sc_list(ioverlap,ioverlap_last)
-c        write (iout,*) 'Overlaping residues ',ioverlap_last,
-c     &           (ioverlap(j),j=1,ioverlap_last)
+        write (iout,*) 'Overlaping residues ',ioverlap_last,
+     &           (ioverlap(j),j=1,ioverlap_last)
       enddo
 
       if (k.le.1000.and.ioverlap_last.eq.0) then
@@ -845,10 +845,12 @@ c     &           (ioverlap(j),j=1,ioverlap_last)
       integer ioverlap(maxres),ioverlap_last
       data redfac /0.5D0/
 
+c      write (iout,*) "overlap_sc_list"
       ioverlap_last=0
 C Check for SC-SC overlaps and mark residues
 c      print *,'>>overlap_sc nnt=',nnt,' nct=',nct
       ind=0
+c      write(iout,*) "iatsc_s",iatsc_s," iatsc_e",iatsc_e
       do i=iatsc_s,iatsc_e
         itypi=iabs(itype(i))
         itypi1=iabs(itype(i+1))
@@ -883,7 +885,7 @@ c            write (iout,*) "i,j",i,j," itypi,itypj",itypi,itypj
           else 
            rcomp=sigma(itypi,itypj)
           endif
-c         print '(2(a3,2i3),a3,2f10.5)',
+c          write (iout,'(2(a3,2i3),a3,2f10.5)'),
 c     &        ' i=',i,iti,' j=',j,itj,' d=',dist(nres+i,nres+j)
 c     &        ,rcomp
             xj=c(1,nres+j)-xi
@@ -895,17 +897,23 @@ c     &        ,rcomp
             rrij=1.0D0/(xj*xj+yj*yj+zj*zj)
             rij=dsqrt(rrij)
             call sc_angular
+c            write (iout,*) "dxj",dxj," dyj",dyj," dzj",dzj
+c            write (iout,*) "erij",erij
+c            write (iout,*) "om1",om1," om2",om2," om12",om12,
+c     &       " faceps1",faceps1," eps1",eps1
+c            write (iout,*) "sigsq",sigsq
             sigsq=1.0D0/sigsq
             sig=sig0ij*dsqrt(sigsq)
             rij_shift=1.0D0/rij-sig+sig0ij
-
-ct          if ( 1.0/rij .lt. redfac*rcomp .or. 
-ct     &       rij_shift.le.0.0D0 ) then
+c            write (iout,*) "rij_shift",rij_shift
+c          if ( 1.0/rij .lt. redfac*rcomp .or. 
+c     &       rij_shift.le.0.0D0 ) then
+c          write (iout,'(a,i3,a,i3,a,f10.5,a,3f10.5)')
+c     &     'overlap SC-SC: i=',i,' j=',j,
+c     &     ' dist=',dist(nres+i,nres+j),' rcomp=',
+c     &     rcomp,1.0/rij,rij_shift
             if ( rij_shift.le.0.0D0 ) then
-cd           write (iout,'(a,i3,a,i3,a,f10.5,a,3f10.5)')
-cd     &     'overlap SC-SC: i=',i,' j=',j,
-cd     &     ' dist=',dist(nres+i,nres+j),' rcomp=',
-cd     &     rcomp,1.0/rij,rij_shift
+c                    write (iout,*) "overlap",i,j
           ioverlap_last=ioverlap_last+1
           ioverlap(ioverlap_last)=i         
           do k=1,ioverlap_last-1