input data has gathers along *4th* axis;
velocity file contains slowness squared with n2=3 (Wavg,Wcos,Wsin);
offset file contains x,y offset pairs for input data
Parameters
float eps=0.01
stretch regularization
int extend=8
trace extension
bool half=y [y/n]
if y, the second and third axes are half-offset instead of full offset