projects
/
unres.git
/ blobdiff
commit
grep
author
committer
pickaxe
?
search:
re
summary
|
shortlog
|
log
|
commit
|
commitdiff
|
tree
raw
|
inline
| side by side
program genereating HOMOL_FRAG consensus fragments
[unres.git]
/
source
/
xdrfpdb
/
src-M
/
xdrf2pdb-m.F
diff --git
a/source/xdrfpdb/src-M/xdrf2pdb-m.F
b/source/xdrfpdb/src-M/xdrf2pdb-m.F
index
118148a
..
24be584
100644
(file)
--- a/
source/xdrfpdb/src-M/xdrf2pdb-m.F
+++ b/
source/xdrfpdb/src-M/xdrf2pdb-m.F
@@
-13,7
+13,7
@@
character*3 licz
character*8 ucase
external ucase
character*3 licz
character*8 ucase
external ucase
- logical oneletter, iblnk
+ logical oneletter,iblnk
integer rescode
external rescode
integer rescode
external rescode
@@
-40,8
+40,7
@@
10 continue
nres=i
i=0
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
i=i+1
enddo
nres=i
@@
-53,8
+52,7
@@
c do while (.not.iblnk(sequenc(i+1)(1:1)))
11 continue
nres=i
i=0
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
i=i+1
enddo
nres=i
@@
-99,9
+97,9
@@
c print *,"ifreq",ifreq," ntraj",ntraj," itraj",itraj
kk = kk + 1
call xdrffloat(ixdrf, potE, iret)
call xdrffloat(ixdrf, uconst, iret)
kk = kk + 1
call xdrffloat(ixdrf, potE, iret)
call xdrffloat(ixdrf, uconst, iret)
-#ifdef NEWUNRES
- call xdrffloat(ixdrf, uconst_back, iret)
-#endif
+c#ifdef NEWUNRES
+c call xdrffloat(ixdrf, uconst_back, iret)
+c#endif
call xdrffloat(ixdrf, t_bath, iret)
call xdrfint(ixdrf, nss, iret)
do j=1,nss
call xdrffloat(ixdrf, t_bath, iret)
call xdrfint(ixdrf, nss, iret)
do j=1,nss
@@
-109,11
+107,12
@@
c print *,"ifreq",ifreq," ntraj",ntraj," itraj",itraj
call xdrfint(ixdrf, jhpb(j), iret)
enddo
call xdrfint(ixdrf, nfrag, iret)
call xdrfint(ixdrf, jhpb(j), iret)
enddo
call xdrfint(ixdrf, nfrag, iret)
+ call xdrfint(ixdrf, iset, iret)
do i=1,nfrag
call xdrffloat(ixdrf, qfrag(i), iret)
enddo
prec=10000.0
do i=1,nfrag
call xdrffloat(ixdrf, qfrag(i), iret)
enddo
prec=10000.0
-
+c print *,'nss=',nss,'nfrag=',nfrag,'iset=',iset
isize=0
call xdrf3dfcoord(ixdrf, coord, isize, prec, iret)
isize=0
call xdrf3dfcoord(ixdrf, coord, isize, prec, iret)
@@
-139,7
+138,8
@@
c write (*,'(8f10.5)') ((coord(k,j),k=1,3),j=1,isize)
enddo
enddo
etot=potE
enddo
enddo
etot=potE
- write (tytul,'(a,i6)') "Structure",kk
+ write (tytul,'(a1,i6,a8,f6.1,a6,f10.1)') "#",kk," t_bath ",
+ & t_bath, " time ",time
call pdbout(etot,tytul,9)
endif
enddo
call pdbout(etot,tytul,9)
endif
enddo