X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Fxdrfpdb%2Fsrc%2Fxdrf2pdb.F;h=1f7ce8e8e94d8c5de4a7d913812ba1a82aa56ccd;hb=e26610dc1d9f408822f8530f2dca7a04fc36c637;hp=592d552974487796ab38fcdf916b547476518a54;hpb=3a334283df7e6ba60b11f1475278af4cee5bedbc;p=unres.git

diff --git a/source/xdrfpdb/src/xdrf2pdb.F b/source/xdrfpdb/src/xdrf2pdb.F
index 592d552..1f7ce8e 100644
--- a/source/xdrfpdb/src/xdrf2pdb.F
+++ b/source/xdrfpdb/src/xdrf2pdb.F
@@ -3,6 +3,7 @@
       include 'COMMON.CHAIN'
       include 'COMMON.INTERACT'
       include 'COMMON.SBRIDGE'
+      include 'COMMON.DISTFIT'
       real*4 coord(3,2*maxres)
       real*4 prec,time,potE,uconst,t_bath,qfrag(100)
       real*8 etot
@@ -94,37 +95,44 @@
       if (itype(nres).eq.21) nct = nres-1
       print *,"nnt",nnt," nct",nct
       print *,"file",arg
+      do i=nnt,nct
+        if (itype(i).ne.20) then
+          itel(i)=1
+        else
+          itel(i)=2
+        endif 
+      enddo
       call xdrfopen(ixdrf,arg, "r", iret)
       if (iret.eq.0) stop
-      print *,"iret",iret
-      print *,"is",is," ie",ie
+c      print *,"iret",iret
+c      print *,"is",is," ie",ie
       kk = 0
       do while(is.eq.0 .or. kk.lt.ie) 
        call xdrffloat(ixdrf, time, iret)
-       print *,"time",time," iret",iret
+c       print *,"time",time," iret",iret
        if(iret.eq.0) exit
        kk = kk + 1
        call xdrffloat(ixdrf, potE, iret)
        call xdrffloat(ixdrf, uconst, iret)
-       print *,"potE",potE," uconst",uconst
+c       print *,"potE",potE," uconst",uconst
 #ifdef NEWUNRES
        call xdrffloat(ixdrf, uconst_back, iret)
-       print *,"uconst_back",uconst_back
+c       print *,"uconst_back",uconst_back
 #endif
        call xdrffloat(ixdrf, t_bath, iret)
-       print *,"t_bath",t_bath
+c       print *,"t_bath",t_bath
        call xdrfint(ixdrf, nss, iret) 
        do j=1,nss
         call xdrfint(ixdrf, ihpb(j), iret)
         call xdrfint(ixdrf, jhpb(j), iret)
        enddo
-       print *,"nss",nss
-       print *,"   ",(ihpb(j),"-",jhpb(j),j=1,nss)
+c       print *,"nss",nss
+c       print *,"   ",(ihpb(j),"-",jhpb(j),j=1,nss)
        call xdrfint(ixdrf, nfrag, iret)
        do i=1,nfrag
         call xdrffloat(ixdrf, qfrag(i), iret)
        enddo
-       print *,"nfrag",nfrag
+c       print *,"nfrag",nfrag
        prec=10000.0
 
        isize=0
@@ -167,6 +175,9 @@ c Calculate the CM
          enddo
          etot=potE
          write (tytul,'(a,i6)') "Structure",kk
+#ifdef SECONDARY
+         call secondary2(.false.)
+#endif
          call pdbout(etot,tytul,9)
        endif
       enddo