回复 第7楼 的 波波头一头：
INTEX.path = "E:\\Flexpart\\INTEXB2006_yy"
para = "CO"
flex.dir = "D:\\1"
res.fn = "D:\\1\\test.csv"
outt = seq(3,24*20,3)
##########################
paths = dir(flex.dir, "^[0-9]+[0-9]$", full.names=TRUE)
res = simulate.vmr.ts(paths,E=list(tra=E.tra, pow=E.pow, liv=E.liv, ind=E.ind), outt = outt, file=res.fn, overwrite=FALSE,
hour.profile=list(tra.hour, pow.hour, liv.hour, ind.hour),
week.profile=list(tra.week, pow.week, liv.week, ind.week),
month.profile=list(tra.month, pow.month, liv.month, ind.month))
##############################
tmp = readfile.csv(res.fn)
#tmp$dtm = tmp$dtm + 8/24
tmp$total = apply(as.matrix(tmp[, c("tra_168_1", "pow_168_1", "liv_168_1", "ind_168_1")]), 1, sum)
plot.time.series(tmp, c("tra_168_1", "pow_168_1", "liv_168_1", "ind_168_1", "total"), ylim=c(0,1500))
这是代码关于tmp部分~您看看~
