projects
/
unres.git
/ blobdiff
commit
grep
author
committer
pickaxe
?
search:
re
summary
|
shortlog
|
log
|
commit
|
commitdiff
|
tree
raw
|
inline
| side by side
working shielding
[unres.git]
/
source
/
unres
/
src_MD-M
/
readpdb.F
diff --git
a/source/unres/src_MD-M/readpdb.F
b/source/unres/src_MD-M/readpdb.F
index
5da38bd
..
3caff57
100644
(file)
--- a/
source/unres/src_MD-M/readpdb.F
+++ b/
source/unres/src_MD-M/readpdb.F
@@
-80,7
+80,7
@@
C Calculate the CM of the preceding residue.
endif
C Start new residue.
c write (iout,'(a80)') card
endif
C Start new residue.
c write (iout,'(a80)') card
- read (card(24:26),*) ires
+ read (card(23:26),*) ires
read (card(18:20),'(a3)') res
if (ibeg.eq.1) then
ishift=ires-1
read (card(18:20),'(a3)') res
if (ibeg.eq.1) then
ishift=ires-1
@@
-246,6
+246,7
@@
C Calculate internal coordinates.
& (c(j,nres+ires),j=1,3)
enddo
endif
& (c(j,nres+ires),j=1,3)
enddo
endif
+C print *,"before int_from_cart"
call int_from_cart(.true.,.false.)
call sc_loc_geom(.true.)
do i=1,nres
call int_from_cart(.true.,.false.)
call sc_loc_geom(.true.)
do i=1,nres
@@
-369,7
+370,6
@@
cc enddiag
hfrag(i,j)=hfrag(i,j)-ishift
enddo
enddo
hfrag(i,j)=hfrag(i,j)-ishift
enddo
enddo
-
return
end
c---------------------------------------------------------------------------
return
end
c---------------------------------------------------------------------------
@@
-437,6
+437,7
@@
c vbld(nres)=3.8d0
c vbld_inv(nres)=1.0d0/vbld(2)
c endif
c endif
c vbld_inv(nres)=1.0d0/vbld(2)
c endif
c endif
+ print *,"A TU2"
if (lside) then
do i=2,nres-1
do j=1,3
if (lside) then
do i=2,nres-1
do j=1,3