Synopsis
		sfdivn2d < inp.rsf den=den.rsf > out.rsf rect1=1 rect2=1 niter=100 verb=n

Parameters
		file den= auxiliary input file name   int niter=100 number of iterations   int rect1=1   int rect2=1 smoothing radius   bool verb=n [y/n] verbosity