X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Fxdrfpdb%2Fsrc-M%2Fxdrf2pdb-m.F;h=ec14da2183ab9ab3af6872babb1afd66ded0229e;hb=HEAD;hp=b665b89747e09e38a135ebd70a468e7b696c1549;hpb=03796abc37f8aeff6e6c84e967993ce1b06f05a4;p=unres.git

diff --git a/source/xdrfpdb/src-M/xdrf2pdb-m.F b/source/xdrfpdb/src-M/xdrf2pdb-m.F
index b665b89..ec14da2 100644
--- a/source/xdrfpdb/src-M/xdrf2pdb-m.F
+++ b/source/xdrfpdb/src-M/xdrf2pdb-m.F
@@ -10,17 +10,18 @@
       character*3 sequenc(maxres)
       character*50 tytul
       character*8 onethree,cfreq,cntraj,citraj
+      character*3 licz
       character*8 ucase
       external ucase
-      logical oneletter
+      logical oneletter,iblnk
       integer rescode
       external rescode
       
       ifreq=1
-      if (iargc().lt.3) then
+      if (iargc().lt.5) then
         print '(2a)',
-     &   "Usage: xdrf2pdb-m one/three seqfile cxfile [freq] [pdbfile] ",
-     &    " [ntraj] [itraj]"
+     &   "Usage: xdrf2pdb-m one/three seqfile cxfile ntraj itraj",
+     &   " [pdbfile] [freq]"
         stop
       endif
       call getarg(1,onethree)
@@ -39,8 +40,7 @@
    10   continue
         nres=i
         i=0
-c        do while (.not.iblnk(sequenc(i+1)(1:1)))
-        do while (.not.(iblnk(sequenc(i+1)(1:1))==0) )
+        do while (.not.iblnk(sequenc(i+1)(1:1)))
           i=i+1
         enddo 
         nres=i
@@ -52,8 +52,7 @@ c        do while (.not.iblnk(sequenc(i+1)(1:1)))
    11   continue
         nres=i
         i=0
-c        do while (.not.iblnk(sequenc(i+1)(1:1)))
-        do while (.not.(iblnk(sequenc(i+1)(1:1))==0))
+        do while (.not.iblnk(sequenc(i+1)(1:1)))
           i=i+1
         enddo 
         nres=i
@@ -61,7 +60,7 @@ c        do while (.not.iblnk(sequenc(i+1)(1:1)))
           itype(i)=rescode(i,sequenc(i),0)
         enddo
         print *,nres
-        print '(a3,1x)',(sequenc(i),i=1,nres)
+        print '(20(a3,1x))',(sequenc(i),i=1,nres)
       endif
       call getarg(3,arg)
       iext = index(arg,'.cx') - 1
@@ -69,28 +68,21 @@ c        do while (.not.iblnk(sequenc(i+1)(1:1)))
         print *,"Error - not a cx file"
         stop
       endif
-      if (iargc().gt.3) then
-        call getarg(4,cfreq)
-        read (cfreq,*) ifreq
-      endif
-      if (iargc().gt.4) then
-        call getarg(5,pdbfile)
-      else
-        pdbfile=arg(:iext)//'.pdb'
-      endif
+      call getarg(4,cntraj)
+      read (cntraj,*) ntraj
+      call getarg(5,citraj)
+      read (citraj,*) itraj
       if (iargc().gt.5) then
-        call getarg(6,cntraj)
-        read (cntraj,*) ntraj
+        call getarg(6,pdbfile)
       else
-        ntraj=1
+        write(licz,'(bz,i3.3)') itraj
+        pdbfile=arg(:iext)//'_'//licz//'.pdb'
       endif
       if (iargc().gt.6) then
-        call getarg(7,citraj)
-        read (citraj,*) itraj
-      else
-        itraj=1
+        call getarg(7,cfreq)
+        read (cfreq,*) ifreq
       endif
-      print *,"ifreq",ifreq," ntraj",ntraj," itraj",itraj
+c      print *,"ifreq",ifreq," ntraj",ntraj," itraj",itraj
       open(9,file=pdbfile)
       nnt = 1
       if (itype(1).eq.21) nnt = 2