Synopsis
		sfpyran < pi.rsf > po.rsf axis=2 range=1 seed=n2 type=y mean=0 var=1 rep=n

Parameters
		int axis=2   float mean=0 noise mean (default=0)   float range=1 noise range (default=1)   bool rep=n [y/n] if y, replace data with noise   int seed=n2 random seed (default=n2)   string type=y noise type, y: normal, n: uniform   float var=1 noise variance (default=1)