X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Funres%2Fsrc_MD-M%2Freadpdb.F;h=3caff57720b1a5aae3dcf2cc4125013549aaf5a3;hb=7a3dca9a7528511e9cc31548c1439c3feb47268c;hp=5da38bdc2d5a0286fbd0ebaa2b60ed0a3fd70e82;hpb=4b5e28ba53f1fd9303e9cf263c8eea9ee0df1583;p=unres.git diff --git a/source/unres/src_MD-M/readpdb.F b/source/unres/src_MD-M/readpdb.F index 5da38bd..3caff57 100644 --- 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 - read (card(24:26),*) ires + read (card(23:26),*) ires 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 print *,"before int_from_cart" 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 - return end c--------------------------------------------------------------------------- @@ -437,6 +437,7 @@ c vbld(nres)=3.8d0 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