Synopsis
		sfoway1 < in.rsf > out.rsf velocity=vel.rsf vgrad=vgrad.rsf eps=0.1 lagrange=n

Axis order: x, p, t

Parameters
		float eps=0.1 stretch regularization   bool lagrange=n [y/n] Use Lagrangian method   file velocity= auxiliary input file name   file vgrad= auxiliary input file name