PRO Gt_get_data, X, Y, Z, Xf, Yf, Xc, Yc

; GT_GET_DATA: PV-WAVE procedure file to read original data

; files to be used for GTGRID testing. The file 'gt_fault.dat' 

; contains data for faults, while 'gt_xyz.dat' contains 

; X, Y and Z data.

;  Invoke this procedure with the command:

;      Gt_get_data, X, Y, Z, Xf, Yf, Xc, Yc

;  where:

;      X  - returned as the X data from gt_xyz.dat

;      Y  - returned as the Y data from gt_xyz.dat

;      Z  - returned as the Z data from gt_xyz.dat

;      Xf - returned as the X Fault data from gt_fault.dat

;      Yf - returned as the Y Fault data from gt_fault.dat

;      Xc - returned as the X Crease data

;           (for this example, same as the X Fault data)

;      Yc - returned as the Y Crease data

;           (for this example, same as the Y Fault data)

X = FLTARR(26)

Y = FLTARR(26)

Z = FLTARR(26)

Xf = FLTARR(60)

Yf = FLTARR(60)

T1 = 0.0

T2 = 0.0

T3 = 0.0

gtgrid_path = extract_path(!Option_path, ”gtgrid”)

OPENR, Unit, gtgrid_path+'/data/gt_fault.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

Xc=Xf

Yc=Yf

OPENR, Unit, gtgrid_path+'/data/gt_xyz.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