Commit ea36b028 authored by weihe's avatar weihe
Browse files

updated

parent fbbffceb
......@@ -34,7 +34,7 @@ datadir="/Storage/ctdas-wrfstilt/new_stilt_footprints/flask/2010/"
eidfilebase="/Storage/ctdas-wrfstilt/stilt_footprints/receptor_info/receptor_info.2010-"
for(n in 1:1)
for(n in 1:12)
{
if(n<10)
mon=paste("0",n,sep="")
......@@ -130,7 +130,7 @@ for(n in 1:1)
}
fnameout1=paste(pfbpath,"FP_BC_3km_time_spent_2.1.before.norm.txt",sep="")
fnameout1=paste(pfbpath,"FP_BC_3km_time_spent_2.1.prep01.norm.txt",sep="")
flog.info("Output to: ",fnameout1, name='logger.b')
outlist = c(num_out3/500, num_out3/500*t_weight1)
out1<-c(ident,eventid,round(outlist,6),flag)
......@@ -146,7 +146,7 @@ for(n in 1:1)
# for every month, normalize
fnameout2=paste(pfbpath,"FP_BC_3km_time_spent_2.1.norm.txt",sep="")
fnameout2=paste(pfbpath,"FP_BC_3km_time_spent_2.1.prep02.norm.txt",sep="")
flog.info("Output to: ",fnameout2, name='logger.b')
footfile<-fnameout1
......@@ -234,7 +234,7 @@ for(n in 1:1)
print(amtsum)
print(bcf[i])
amtsum=amtsum+bcf[i]
amtsum=amtsum+as.numeric(bcf[i])
amtcount=amtcount+1
}
}
......@@ -242,9 +242,10 @@ for(n in 1:1)
newfile3=T
for (i in 1:lin3)
{
bcf1[i]=bcf[i]/(amtsum/amtcount)
tmp=bcf[i]
bcf[i]=as.numeric(bcf[i])/(amtsum/amtcount)
out3=c(ident[i],evn1[i],bcf[i],bcf1[i],sitenew,airtowflag)
out3=c(ident[i],evn1[i],tmp, round(bcf[i],6), site[i], type[i])
header=c("ident","eventid","value","normvalue","site","type")
if(newfile3)
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment