X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Fxdrfpdb%2Fsrc%2Fxdrf2pdb-m.F;h=a12edb8b3753927f8c5e1ff262c273c11e092792;hb=babec4c681888e214510d01818186b2ab6bf49d9;hp=0016d40b037ee216637f20c60866a95198d3fb5a;hpb=a45ce81aae2c857c30b821e39e45a7a08356de4d;p=unres.git diff --git a/source/xdrfpdb/src/xdrf2pdb-m.F b/source/xdrfpdb/src/xdrf2pdb-m.F index 0016d40..a12edb8 100644 --- a/source/xdrfpdb/src/xdrf2pdb-m.F +++ b/source/xdrfpdb/src/xdrf2pdb-m.F @@ -12,7 +12,7 @@ character*8 onethree,cfreq,cntraj,citraj character*8 ucase external ucase - logical oneletter + logical oneletter, iblnk integer rescode external rescode @@ -40,6 +40,7 @@ nres=i i=0 do while (.not.iblnk(sequenc(i+1)(1:1))) +c do while (.not.(iblnk(sequenc(i+1)(1:1)) == 0)) i=i+1 enddo nres=i @@ -52,6 +53,7 @@ nres=i i=0 do while (.not.iblnk(sequenc(i+1)(1:1))) +c do while (.not.(iblnk(sequenc(i+1)(1:1)) == 0)) i=i+1 enddo nres=i @@ -103,10 +105,12 @@ kk = kk + 1 call xdrffloat(ixdrf, potE, iret) call xdrffloat(ixdrf, uconst, iret) + call xdrffloat(ixdrf, t_bath, iret) + print *,"potE",potE," uconst",uconst," t_bath",t_bath #ifdef NEWUNRES call xdrffloat(ixdrf, uconst_back, iret) #endif - call xdrffloat(ixdrf, t_bath, iret) + print *,"uconst_back",uconst_back call xdrfint(ixdrf, nss, iret) do j=1,nss call xdrfint(ixdrf, ihpb(j), iret) @@ -119,6 +123,7 @@ prec=10000.0 isize=0 + print *," call xdrf3coord" call xdrf3dfcoord(ixdrf, coord, isize, prec, iret)