d252 = "DAY252.1999" library(RSEIS) load("/home/lees/LAB/DB99.RDATA") nn = scan(file=d252, list(j="", i=0, yr=0, jd=0, mo=0, dom=0, hr=0, mi=0, sec=0, dur=0, tag="")) at = JtimL(nn) for(i in 1:length(at)) { at1 = at[i]-0.5/24 at2 = at1+1/24 TH = Mine.seis(at1, at2, DB99, NULL, NULL , kind = 1, Iendian=1, BIGLONG=FALSE, CHOP=TRUE, verbose=FALSE) gret = swig(TH) if(gret$but=="QUIT") break }