Synopsis
		sfkirmodcmp < curv.rsf > modl.rsf head=head.rsf refl=refl.rsf cmp=n nt= dt=0.004 t0=0. nsx=nx s0x=x0 dsx=dx nsy=ny s0y=y0 dsy=dy nmx=(nx+3)/2.0 m0x=x0 dmx=2*dx nmy=(ny+3)/2.0 m0y=y0 dmy=2*dy nhx=nx h0x=0. dhx=dx nhy=ny h0y=0. dhy=dy r0=1. aper=hypotf(nx*dx,ny*dy) freq=0.2/dt rgrad= dipx= dipy= type=

Parameters
		float aper=hypotf(nx*dx,ny*dy) aperture   bool cmp=n [y/n] CMP or shot gather output. Default is shot gather.   float dhx=dx inline offset increment (if cmp=y, dhx should have same sign as h0x)   float dhy=dy crossline offset increment (if cmp=y, dhy should have same sign as h0y)   string dipx=   string dipy=   float dmx=2*dx inline CMP increment   float dmy=2*dy crossline CMP increment   float dsx=dx inline shot increment   float dsy=dy crossline shot increment   float dt=0.004 time sampling   float freq=0.2/dt peak frequency for Ricker wavelet   float h0x=0. first inline offset   float h0y=0. first crossline offset   string head= source-receiver geometry (optional) (auxiliary input file name)   float m0x=x0 first inline CMP   float m0y=y0 first crossline CMP   int nhx=nx number of inline offsets   int nhy=ny number of crossline offsets   int nmx=(nx+3)/2.0 number of inline CMPs   int nmy=(ny+3)/2.0 number of crossline CMPs   int nsx=nx number of inline shots   int nsy=ny number of crossline shots   int nt= time samples   float r0=1. constant reflectivity   string refl= auxiliary input file name   string rgrad=   float s0x=x0 first inline shot   float s0y=y0 first crossline shot   float t0=0. time origin   string type= type of velocity ('c': constant, 's': linear sloth, 'v': linear velocity)