c234567
      program popread
c
c     This is a simple program to read population data
c     for a single population file (from Bob Eckert's files).
c
      character fname*44/'/usr3/population/dat/pop33x119.dat'/
      character cfile*6
      integer*2 popdat(3, 0:5000)
      pointer  (p,popdat)

c    1 write(6,5) 
c    5 format(1h ,'Enter coordinate file (e.g. 33x119)')
      read(5,7)  cfile
    7 format(a)
      fname='/net/oets1/usr3/population/dat/pop'//cfile//'.dat'
c      write(6,8) fname
c    8 format(1h ,a)
      open(38,file=fname,
     >        status='old',form='unformatted',access='direct',
     >        recl=4,err=100)
      read(38,rec=1) numrecords
c      write(6,10) fname,numrecords
c   10 format(1h ,a,' has ',i8,' records.')
      close(38)
      open(39,file=fname,status='old',form='unformatted',
     >        access='direct',recl=6*numrecords)
      p=malloc(6*numrecords)
      last_record=numrecords-1
      read(39,rec=1) ((popdat(i,j), i=1,3), j=0, last_record)
      read(cfile,15) lat
   15 format(i2.2)
      read(cfile,16) lon
   16 format(3x,i3.3)
c      write(6,17) lat, lon
c   17 format(1h ,'Coordinates of 1-deg by 1-deg file: (',i2.2,',',
c     >                                                   i3.3,')')
      do  j=1,last_record
      write(6,20) real(lat)+popdat(1,j)/10000.0,
     >            real(lon)+popdat(2,j)/10000.0,
     >                      popdat(3,j)   
   20 format(f7.4,',',f8.4,',',i5.5)
      end do
  100 end