Program reformat_bathy
*
* Note:  THIS FILE MUST BE EDITED BY HAND AT THE LINES 
*        INDICATED BY ######### SIGNS.
*
*############## ARE THESE DIMENSIONS LARGE ENOUGH?
      real z_cb(2000,2000)
      character*80 filein_xyz, fileout_bin

*#########################
      filein_xyz = 'epr4_merge.xyz'
      fileout_bin = 'epr4_merge.bin'
      
*###############  SET THE HEADER VALUES FOR YOUR MAP; GET THESE FROM GRDINFO (GMT)
      rlongmin = -104.75
      rlongmax = -103.833
      rlatmin  = 8.3333
      rlatmax  = 10.3333
      gsmin    = .081374
      ncol_cb  = 671
      nrow_cb  = 1107      
      
      open (unit =2 , file= filein_xyz ,status='old')
      
*  GMT writes out data by rows, starting from the TOP!
*  Read the data
      do i = 1, nrow_cb
        do j = 1, ncol_cb
          read(2,*) dumx,dumy, z_cb(i,j)
        end do
      end do

      open (unit    = 1,
     -      file    = fileout_bin,
     -      status  = 'new',
     -      form    = 'unformatted')
      
      write (1) rlongmin,rlongmax,rlatmin,rlatmax,gsmin,
     -          ncol_cb,nrow_cb
      do i = nrow_cb,1,-1
         write (1) (z_cb(i,j), j=1,ncol_cb)
      end do
      close (unit = 1)
      end