pro gt_con_data, x, y, z, xf, yf

x = fltarr(626)

y = fltarr(626)

z = fltarr(626)

xf = fltarr(69)

yf = fltarr(69)

t1 = 0.0

t2 = 0.0

t3 = 0.0

gtgrid_path = extract_path(!Option_path, ”gtgrid”)

openr, unit, gtgrid_path+'/data/gt_faults.dat', /get_lun

FOR i = 0, n_elements(xf)-1 DO BEGIN

readf, unit, t1, t2

xf(i) = float(t1)

yf(i) = float(t2)

ENDFOR

free_lun, unit

openr, unit, gtgrid_path+'/data/gt_contours.dat', $
/get_lun

FOR i = 0, n_elements(x)-1 DO BEGIN

readf, unit, t1, t2, t3

x(i) = float(t1)

y(i) = float(t2)

z(i) = float(t3)

ENDFOR

free_lun, unit

END