#include "convert_datum.h"

#define GRID_TOLERANCE 2000	/* meters */
#define BATH_UTIL_DEF

#define BATHYMETRIC_PROG "/project0/pmd/mfp/bathymetric-data/extract_bath"
#define BATHYMETRIC_FILE "/acoustics1/pmd/sbcxbath.dos.utm"
#define BATHYMETRIC_SAMPLE_INTERVAL 250.0	/* meters */

#ifndef BUFFER_LENGTH
#define BUFFER_LENGTH 4096
#endif

typedef struct bathymetric_point_struct {
  double r,x,y,z;
} BathymetricPoint;

typedef struct bathymetry_database_struct {
  int npts;
  BathymetricPoint *p;
  char grid_zone[GRID_ZONE_LENGTH];
} Bathymetry;

typedef struct bathymetric_struct {
  int nr;
  BathymetricPoint *bp;
} Bathymetric;


Bathymetry *load_bathymetric_database (FILE *fp);
int bathymetric_compare (BathymetricPoint *, BathymetricPoint *);
double find_bathymetry (UTM, Bathymetry *);

Bathymetric *obtain_bathymetric_data (double, double, double, double, unsigned int);
Bathymetric *reverse_bathymetric_data (char *, double, double, double, double);
Bathymetric *read_bathymetric_data (char *, double, double);