Synopsis
		sfhshape < in.rsf filt=filt.rsf > out.rsf n= adj=n ns= lag=

Parameters
		bool adj=n [y/n] if y, do adjoint operation   file filt= auxiliary input file name   string lag=   ints n= [dim]   int ns= scaling