cmdd = paste(sep=' ', "ttBK 0 ", i , " > ttBK.tempoutput")
system(cmdd )
A = scan(file="ttBK.tempoutput", list(num=0, phase="", sec=0, tmin=0, tsec=0, dtdd=0, dtdh=0, dddp=0), skip=1)
attr(A, 'name') = "NPL"
attr(A, 'input') = list(delta=i , depth=0)
Run the code above in your browser using DataLab