hdrfile=j.hdr.h 
niter=10 
s1lags=0,1,2,3,4,5,6
s2lags=0,1,2,3,4,5,6
