| 1 | c ......................
|
|---|
| 2 | c contents of a PDB file
|
|---|
| 3 | c ......................
|
|---|
| 4 |
|
|---|
| 5 | parameter (MXCHP =100, ! max. no. of polypeptide chains
|
|---|
| 6 | # MXRSP =1000, ! max. no. of residues
|
|---|
| 7 | # MXATP =10000) ! max. no. of atoms
|
|---|
| 8 |
|
|---|
| 9 | c nchp - no. of polypeptide chains
|
|---|
| 10 | c nchrsp() - no. of residues / chain
|
|---|
| 11 | c nrsp - total no. of residues
|
|---|
| 12 | c irsatp() - index of 1st atom / res.
|
|---|
| 13 | c nrsatp() - no. of atoms / res.
|
|---|
| 14 | c natp - total no. of selected atoms
|
|---|
| 15 | c noatp() - atom numbers
|
|---|
| 16 |
|
|---|
| 17 | common /pdb_i/ nchp,nchrsp(MXCHP),
|
|---|
| 18 | # nrsp,irsatp(MXRSP),nrsatp(MXRSP),
|
|---|
| 19 | # natp,noatp(MXATP)
|
|---|
| 20 | save /pdb_i/
|
|---|
| 21 |
|
|---|
| 22 | c chnp() - chain identifiers
|
|---|
| 23 | c rsidp() - residue identifiers (number + insertion code)
|
|---|
| 24 | c rsnmp() - residues (sequence, 3-letter code)
|
|---|
| 25 | c atnmp() - atom names
|
|---|
| 26 |
|
|---|
| 27 | character chnp(MXCHP),
|
|---|
| 28 | # rsidp(MXRSP)*5,rsnmp(MXRSP)*3,
|
|---|
| 29 | # atnmp(MXATP)*4
|
|---|
| 30 |
|
|---|
| 31 | common /pdb_c/ chnp,rsnmp,rsidp,atnmp
|
|---|
| 32 | save /pdb_c/
|
|---|
| 33 |
|
|---|
| 34 | c xatp,yatp,zatp - atom coordinates
|
|---|
| 35 |
|
|---|
| 36 | common /pdb_r/ xatp(MXATP),yatp(MXATP),zatp(MXATP)
|
|---|
| 37 | save /pdb_r/
|
|---|
| 38 |
|
|---|
| 39 | c ------------------- code to list all PDB information
|
|---|
| 40 | c ir=0
|
|---|
| 41 | c do i=1,nchp
|
|---|
| 42 | c write(*,*) ' ===== chain |',chnp(i),'|'
|
|---|
| 43 | c do j=1,nchrsp(i)
|
|---|
| 44 | c ir=ir+1
|
|---|
| 45 | c write(*,*) ' ----- ',rsidp(ir),' ',rsnmp(ir),' ',nrsatp(ir)
|
|---|
| 46 | c k1=irsatp(ir)
|
|---|
| 47 | c k2=k1+nrsatp(ir)-1
|
|---|
| 48 | c do k=k1,k2
|
|---|
| 49 | c write(*,*) ' ',noatp(k),' ',atnmp(k),' ',(xyzp(l,k),l=1,3)
|
|---|
| 50 | c enddo
|
|---|
| 51 | c enddo
|
|---|
| 52 | c enddo
|
|---|
