Changeset 32289cd for pdbread.f


Ignore:
Timestamp:
11/19/09 11:29:41 (14 years ago)
Author:
baerbaer <baerbaer@…>
Branches:
master
Children:
38d77eb
Parents:
6650a56
Message:

Explicitly declare variables.

git-svn-id: svn+ssh://svn.berlios.de/svnroot/repos/smmp/trunk@33 26dc1dd8-5c4e-0410-9ffe-d298b4865968

File:
1 edited

Legend:

Unmodified
Added
Removed
  • pdbread.f

    r6650a56 r32289cd  
    44!
    55! Copyright 2003-2005  Frank Eisenmenger, U.H.E. Hansmann,
    6 !                      Shura Hayryan, Chin-Ku 
     6!                      Shura Hayryan, Chin-Ku
    77! Copyright 2007       Frank Eisenmenger, U.H.E. Hansmann,
    88!                      Jan H. Meinke, Sandipan Mohanty
     
    2121! ......................................................
    2222
    23       implicit real*8 (a-h,o-z)
    24       implicit integer*4 (i-n)
    25 
    2623      include 'INCP.H'
     24     
     25      double precision cor
     26
     27      integer ier, l, iendst, lunpdb, io, iopfil, iat, i
    2728
    2829! -------------------------- input
     
    4748        lunpdb = 99
    4849      else
    49         write (*,'(a)') 
     50        write (*,'(a)')
    5051     &    ' pdbread> empty file name to read pdb-structure'
    5152
     
    5657
    5758      if (io.le.0) then
    58         write (*,'(a,/,a)') 
     59        write (*,'(a,/,a)')
    5960     &  ' pdbread> ERROR opening file to read pdb-structure: ',
    6061     &  pdbfil(1:iendst(pdbfil))
     
    6970
    7071      if ( line(17:17).ne.' ' )  then
    71         write (*,'(a,/,a,/,a,/,2a)') 
     72        write (*,'(a,/,a,/,a,/,2a)')
    7273     &  ' pdbread> found alternate atom location: ',
    7374     &  '                !',
     
    8586
    8687      if ((natp+1).gt.MXATP) then
    87         write (*,'(a,i5,a,/,a)') 
     88        write (*,'(a,i5,a,/,a)')
    8889     &  ' pdbread>  >MXATP (',MXATP,') ATOM lines in PDB file ',
    8990     &  pdbfil(1:iendst(pdbfil))
     
    9697
    9798        if ((nchp+1).gt.MXCHP) then
    98           write (*,'(a,i3,a,/,a)') 
     99          write (*,'(a,i3,a,/,a)')
    99100     &    ' pdbread>  >MXCHP (',MXCHP,') chains in PDB file ',
    100101     &    pdbfil(1:iendst(pdbfil))
     
    105106
    106107        if ((nrsp+1).gt.MXRSP) then
    107           write (*,'(a,i3,a,/,a)') 
     108          write (*,'(a,i3,a,/,a)')
    108109     &    ' pdbread>  >MXRSP (',MXRSP,') residues in PDB file ',
    109110     &    pdbfil(1:iendst(pdbfil))
     
    140141
    141142        if ((nrsp+1).gt.MXRSP) then
    142           write (*,'(a,i3,a,/,a)') 
     143          write (*,'(a,i3,a,/,a)')
    143144     &    ' pdbread>  >MXRSP (',MXRSP,') residues in PDB file ',
    144145     &    pdbfil(1:iendst(pdbfil))
     
    171172      goto 1
    172173
    173     2 write (*,'(a,/,a,/,2a)') 
     174    2 write (*,'(a,/,a,/,2a)')
    174175     &  ' pdbread> ERROR reading ATOM line ',
    175176     &  line(:l),
     
    194195      else
    195196
    196         write (*,'(a,/,a)') 
     197        write (*,'(a,/,a)')
    197198     &  ' pdbread> NO atom coordinates selected from file ',
    198199     &  pdbfil(1:iendst(pdbfil))
     
    226227      include 'INCL.H'
    227228      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
    228234
    229235      character res*3
     
    451457      include 'INCP.H'
    452458
     459      integer irs, nml, i1, i2, iat, ix, i
     460
    453461      character*4 atm
    454462
     
    476484              endif
    477485            enddo
    478                
     486
    479487!            write(*,'(8a)') ' pdbvars> ',atm,' not found in '
    480488!     #       ,chnp(nc),' ',rsidp(irs),' ',rsnmp(irs)
     
    484492    1     ixatp(iat)=ix
    485493
    486         enddo   ! SMMP atoms of 'irs' 
     494        enddo   ! SMMP atoms of 'irs'
    487495      enddo   ! residues
    488496
     
    494502      include 'INCL.H'
    495503
     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
    496509      parameter (TOL = 1.d-12)
    497510
     
    499512! -> determine global parameters: shifts dX,dY,dZ
    500513! & angles alpha,beta,gamma [rad], put into 'gbpr'
    501 ! 
     514!
    502515! CALLS: none
    503516!
     
    547560      y2=z3*x1-z1*x3
    548561      y3=z1*x2-z2*x1
    549      
     562
    550563      if ( ( 1.d0 - abs(y3) ) .gt. TOL )  then       ! ============ |beta| < PI/2
    551564
Note: See TracChangeset for help on using the changeset viewer.