- Timestamp:
- 11/19/09 11:29:41 (14 years ago)
- Branches:
- master
- Children:
- 38d77eb
- Parents:
- 6650a56
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
pdbread.f
r6650a56 r32289cd 4 4 ! 5 5 ! Copyright 2003-2005 Frank Eisenmenger, U.H.E. Hansmann, 6 ! Shura Hayryan, Chin-Ku 6 ! Shura Hayryan, Chin-Ku 7 7 ! Copyright 2007 Frank Eisenmenger, U.H.E. Hansmann, 8 8 ! Jan H. Meinke, Sandipan Mohanty … … 21 21 ! ...................................................... 22 22 23 implicit real*8 (a-h,o-z)24 implicit integer*4 (i-n)25 26 23 include 'INCP.H' 24 25 double precision cor 26 27 integer ier, l, iendst, lunpdb, io, iopfil, iat, i 27 28 28 29 ! -------------------------- input … … 47 48 lunpdb = 99 48 49 else 49 write (*,'(a)') 50 write (*,'(a)') 50 51 & ' pdbread> empty file name to read pdb-structure' 51 52 … … 56 57 57 58 if (io.le.0) then 58 write (*,'(a,/,a)') 59 write (*,'(a,/,a)') 59 60 & ' pdbread> ERROR opening file to read pdb-structure: ', 60 61 & pdbfil(1:iendst(pdbfil)) … … 69 70 70 71 if ( line(17:17).ne.' ' ) then 71 write (*,'(a,/,a,/,a,/,2a)') 72 write (*,'(a,/,a,/,a,/,2a)') 72 73 & ' pdbread> found alternate atom location: ', 73 74 & ' !', … … 85 86 86 87 if ((natp+1).gt.MXATP) then 87 write (*,'(a,i5,a,/,a)') 88 write (*,'(a,i5,a,/,a)') 88 89 & ' pdbread> >MXATP (',MXATP,') ATOM lines in PDB file ', 89 90 & pdbfil(1:iendst(pdbfil)) … … 96 97 97 98 if ((nchp+1).gt.MXCHP) then 98 write (*,'(a,i3,a,/,a)') 99 write (*,'(a,i3,a,/,a)') 99 100 & ' pdbread> >MXCHP (',MXCHP,') chains in PDB file ', 100 101 & pdbfil(1:iendst(pdbfil)) … … 105 106 106 107 if ((nrsp+1).gt.MXRSP) then 107 write (*,'(a,i3,a,/,a)') 108 write (*,'(a,i3,a,/,a)') 108 109 & ' pdbread> >MXRSP (',MXRSP,') residues in PDB file ', 109 110 & pdbfil(1:iendst(pdbfil)) … … 140 141 141 142 if ((nrsp+1).gt.MXRSP) then 142 write (*,'(a,i3,a,/,a)') 143 write (*,'(a,i3,a,/,a)') 143 144 & ' pdbread> >MXRSP (',MXRSP,') residues in PDB file ', 144 145 & pdbfil(1:iendst(pdbfil)) … … 171 172 goto 1 172 173 173 2 write (*,'(a,/,a,/,2a)') 174 2 write (*,'(a,/,a,/,2a)') 174 175 & ' pdbread> ERROR reading ATOM line ', 175 176 & line(:l), … … 194 195 else 195 196 196 write (*,'(a,/,a)') 197 write (*,'(a,/,a)') 197 198 & ' pdbread> NO atom coordinates selected from file ', 198 199 & pdbfil(1:iendst(pdbfil)) … … 226 227 include 'INCL.H' 227 228 include 'INCP.H' 229 230 double precision dihedr, rm, av1, av2, rmsd, h, x, y, z 231 232 integer nml, nrs, nc, irb, ire, irs, i, ii, it, i1, i2, i3, i4, j1 233 integer j2, j3, j4, inew, j, ix 228 234 229 235 character res*3 … … 451 457 include 'INCP.H' 452 458 459 integer irs, nml, i1, i2, iat, ix, i 460 453 461 character*4 atm 454 462 … … 476 484 endif 477 485 enddo 478 486 479 487 ! write(*,'(8a)') ' pdbvars> ',atm,' not found in ' 480 488 ! # ,chnp(nc),' ',rsidp(irs),' ',rsnmp(irs) … … 484 492 1 ixatp(iat)=ix 485 493 486 enddo ! SMMP atoms of 'irs' 494 enddo ! SMMP atoms of 'irs' 487 495 enddo ! residues 488 496 … … 494 502 include 'INCL.H' 495 503 504 double precision tol, h1, h2, h3, x1, x2, x3, d, z1, z2, z3, y1 505 double precision y2, y3, yp1, yp2 506 507 integer i1, nml, i2, i3, i 508 496 509 parameter (TOL = 1.d-12) 497 510 … … 499 512 ! -> determine global parameters: shifts dX,dY,dZ 500 513 ! & angles alpha,beta,gamma [rad], put into 'gbpr' 501 ! 514 ! 502 515 ! CALLS: none 503 516 ! … … 547 560 y2=z3*x1-z1*x3 548 561 y3=z1*x2-z2*x1 549 562 550 563 if ( ( 1.d0 - abs(y3) ) .gt. TOL ) then ! ============ |beta| < PI/2 551 564
Note:
See TracChangeset
for help on using the changeset viewer.