C C************************************************************************ C AMBER ** C ** C Copyright (c) 1986, 1991, 1995 ** C Regents of the University of California ** C All Rights Reserved. ** C ** C This software provided pursuant to a license agreement containing ** C restrictions on its disclosure, duplication, and use. This software ** C contains confidential and proprietary information, and may not be ** C extracted or distributed, in whole or in part, for any purpose ** C whatsoever, without the express written permission of the authors. ** C This notice, and the associated author list, must be attached to ** C all copies, or extracts, of this software. Any additional ** C restrictions set forth in the license agreement also apply to this ** C software. ** C************************************************************************ C subroutine getcor (natom, ntx, x) implicit double precision (a-h,o-z) character*80 line c dimension x(*),ititl(20) if (ntx.eq.0) then c c ----- read title etc. c read (53) ititl read (53) natomd if (natomd.ne.natom) then write(6,*) 'Number of atoms in -p and -c files do not agree!' call mexit(6, 1) end if c c ----- Read X c nr3 = natom*3 read (53) (x(i), i=1,nr3) c else c c ----- read title etc. c read (53, '(20a4)') ititl read (53, '(a)' ) line if( line(6:6) .eq. ' ' ) then read( line, '(i5)' ) natomd else read( line, '(i6)' ) natomd end if if (natomd.ne.natom) then write(6,*) 'Number of atoms in -p and -c files do not agree!' call mexit(6, 1) end if c c ----- Read X c nr3 = natom*3 read (53,3) (x(i), i=1,nr3) 3 format(6f12.7) end if c write (6,*) write(6,4) ititl 4 format(' Getting coordinates from file with title:'/5x,20a4) close(53) return end