logical ex
      character*17 fname
      character*22 sitename
      character*4 site
      character*32 junk
      character*8 Lat
      character*9 Long
      character*10 ThA
      character*10 ThB
      character*5 ArA
      character*5 ArB
      character*4 deg
      character*1 sign
      character*8 Alt
      character*5 MSL
      character*8 datval
      character*21 date
      character*11 TauB
      character*11 TauA
      real LatVal,LonVal,MSLAlt,ThAVal,ThBVal,ArAVal,
     C     ArBVal,TBval,TAval


*     Put the location and name of the files into fnameT, 

      data fname/'database_0000.txt'/

      open(file='installed_sites.txt',type='old',unit=22)

      open(file='sites_coordinates.txt',type='old',
     C      access='direct',unit=44)


      write(44,'(A1,A4,1x,A8,14x,A8,5x,A9,4x,A7,
     C     2x,A5,2x,A4,3x,
     C     A5,2x,A4,3x,A4,3x,A4,3x,A8)')'#','site',
     C     'sitename',
     C     'Latitude','Longitude','MSL Alt','Thick',
     C     'Area',
     C     'Thick','Area','TauA','TauB','Cal.Date'
      write(44,'(A1,62x,A5,2x,A4,3x,A5,2x,A4,3x,A4,3x,A4)')
     C     '#','A(cm)','A(m)',
     C      'B(cm)','B(m)','nsec','nsec'

      write(*,*)' '
      write(*,'(A1,A4,1x,A8,14x,A8,5x,A9,4x,A7,2x,A5,2x,A4,3x,
     C     A5,2x,A4,3x,A4,3x,A4,3x,A8)')'#','site','sitename',
     C     'Latitude','Longitude','MSL Alt','Thick',
     C     'Area',
     C     'Thick','Area','TauA','TauB','Cal.Date'
      write(*,'(A1,62x,A5,2x,A4,3x,A5,2x,A4,3x,A4,3x,A4)')
     C     '#','A(cm)','A(m)',
     C      'B(cm)','B(m)','nsec','nsec'


*     after reading installed sites list, put site number
*     into file name.

 100  read(22,'(A4,1x,A22)',end=222)site,sitename

      fname(10:13)=site

      inquire(file=fname,exist=ex)

      if (ex) then
      open(file=fname,type='old',unit=33)
      
      do 200 i=1,8
      read(33,'(A12)')junk
 200  continue
      
      read(33,'(1x,A8,1x,A4,1x,F8.5)')Lat,deg,LatVal
      read(33,'(1x,A9,1x,A4,1x,F10.5)')Long,deg,LonVal

      read(33,'(A12)')junk

      read(33,'(1x,A8,1x,A6,1x,F6.2)')Alt,MSL,MSLAlt

      do 201 i=1,7
      read(33,'(A12)')junk
 201  continue

      read(33,'(1x,A10,1x,F6.3)')ThA,ThAVal
      read(33,'(1x,A5,1x,F6.3)')ArA,ArAVal

      do 205 i=1,3
      read(33,'(A12)')junk
 205  continue

      read(33,'(1x,A11,1x,F6.2)')TauA,TAVal

  
      do 202 i=1,5
      read(33,'(A12)')junk
 202  continue

      read(33,'(1x,A10,1x,F6.3)')ThB,ThBVal
      read(33,'(1x,A5,1x,F6.3)')ArB,ArBVal

      do 203 i=1,3
      read(33,'(A12)')junk
 203  continue

      read(33,'(1x,A11,1x,F6.2)')TauB,TBVal

      do 204 i=1,4
      read(33,'(A12)')junk
 204  continue

      read(33,'(A21,1x,A10)'),date,datval

      write(*,'(A4,1x,A22,1x,F8.5,4x,F10.5,4x,F6.2,2x,F6.3,1x,F6.3,
     C   1x,F6.3,1x,F6.3,1x,F6.2,1x,F6.2,2x,A8)')site,sitename,
     C   LatVal,LonVal,MSLAlt,ThAVal,ArAVal,ThBVal,ArBVal,
     C   TAVal,TBVal,datval

      close(33)

      write(44,'(A4,1x,A22,1x,F8.5,4x,F10.5,4x,F6.2,2x,F6.3,1x,F6.3,
     C   1x,F6.3,1x,F6.3,1x,F6.2,1x,F6.2,2x,A8)')site,sitename,
     C   LatVal,LonVal,MSLAlt,ThAVal,ArAVal,ThBVal,ArBVal,
     C   TAVal,TBVal,datval
     
      endif
      goto 100
 222  continue

      close(44)
      close(22)
      end