Adam's cluster wham update

[unres.git] / source / cluster / wham / src-HCD / wrtclust.f

diff --git a/source/cluster/wham/src-HCD/wrtclust.f b/source/cluster/wham/src-HCD/wrtclust.f
index e21494a..0415640 100644 (file)
--- a/source/cluster/wham/src-HCD/wrtclust.f
+++ b/source/cluster/wham/src-HCD/wrtclust.f
@@ -189,7 +189,7 @@ c create InsightII command file for their displaying in different colors
 c        write (iout,*) "cfname",cfname
         OPEN(ipdb,FILE=CFNAME,STATUS='UNKNOWN',FORM='FORMATTED')
         write (ipdb,'(a,f8.2)') 
-     &    "REMAR AVERAGE CONFORMATIONS AT TEMPERATURE",temper
+     &    "REMARK AVERAGE CONFORMATIONS AT TEMPERATURE",temper
         close (ipdb)
         I=1
         ICON=NCONF(1,1)
@@ -227,6 +227,13 @@ c          write (iout,*) "ncon_out",ncon_out
                 c(k,ii)=allcart(k,ii,icon)
               enddo
             enddo
+            nss=nss_all(icon)
+            write (iout,*) "ICON",icon," nss",nss
+            do k=1,nss
+              ihpb(k)=ihpb_all(k,icon)
+              jhpb(k)=jhpb_all(k,icon)
+              write (iout,*) ihpb(k),jhpb(k)
+            enddo
             call center
             call pdbout(totfree(icon)/beta_h(ib),rmstb(icon),titel)
             write (ipdb,'("TER")')
@@ -240,6 +247,7 @@ c Average structures and structures closest to average
           call ave_coord(i)
           write (ipdb,'(a,i5)') "REMARK CLUSTER",i
           call center
+          nss=0
           call pdbout(totfree_gr(i)/beta_h(ib),rmsave(i),titel)
           write (ipdb,'("TER")')
           if (print_fittest.and.(nsaxs.gt.0 .or. nhpb.gt.0 
@@ -530,6 +538,7 @@ c------------------------------------------------------------------------------
       include 'COMMON.CLUSTER'
       include 'COMMON.CHAIN'
       include 'COMMON.INTERACT'
+      include 'COMMON.SBRIDGE'
       include 'COMMON.VAR'
       include 'COMMON.FFIELD'
       include 'COMMON.TORCNSTR'
@@ -570,6 +579,11 @@ c     & " Edihcnstr",edihcnstr
           c(j,i)=allcart(j,i,jconmin)
         enddo
       enddo
+      nss=nss_all(jconmin)
+      do k=1,nss
+        ihpb(k)=ihpb_all(k,jconmin)
+        jhpb(k)=jhpb_all(k,jconmin)
+      enddo
       return
       end
 c------------------------------------------------------------------------------
@@ -583,6 +597,7 @@ c------------------------------------------------------------------------------
       include 'COMMON.CHAIN'
       include 'COMMON.INTERACT'
       include 'COMMON.VAR'
+      include 'COMMON.SBRIDGE'
       logical non_conv
       double precision przes(3),obrot(3,3)
       integer i,ii,j,k,icon,jcon,jconmin,igr,ipermmin
@@ -616,6 +631,15 @@ c      write (iout,*) "rmsmin",rmsmin," rms",rms
           c(j,i)=allcart(j,i,jconmin)
         enddo
       enddo
+      nss=nss_all(jconmin)
+c      write (iout,*) "jconmin",jconmin," nss",nss
+      call flush(iout)
+      do k=1,nss
+        ihpb(k)=ihpb_all(k,jconmin)
+        jhpb(k)=jhpb_all(k,jconmin)
+c        write (iout,*) "k",k," ihpb",ihpb(k)," jhpb",jhpb(k)
+      enddo
+      call flush(iout)
       return
       end
 c------------------------------------------------------------------------------