|
1
|
library(sp)
|
|
2
|
library(foreign)
|
|
3
|
library(reshape)
|
|
4
|
library(latticeExtra)
|
|
5
|
|
|
6
|
path = "/nobackupp4/aguzman4/climateLayers/output/"
|
|
7
|
tilelist = list.files(path,pattern="^[0-9]")
|
|
8
|
|
|
9
|
rm(dlst)
|
|
10
|
dlst.full=c()
|
|
11
|
for (t in tilelist){
|
|
12
|
|
|
13
|
# t="30.0_-105.0"
|
|
14
|
|
|
15
|
dbf = read.dbf (paste(path,t,"/monthly_covariates_ghcn_data_TMAX_2001_2011",t,".dbf",sep=""))
|
|
16
|
|
|
17
|
dlst=melt(dbf,id.vars=c("station","month","TMax"),measure.vars=paste("mm",sprintf("%02d",1:12),sep="_"))
|
|
18
|
dlst1=melt(dbf,id.vars=c("station","month","TMax"),measure.vars=paste("mm",sprintf("%02d",c(2:12,1)),sep="_"))
|
|
19
|
dlst2=melt(dbf,id.vars=c("station","month","TMax"),measure.vars=paste("mm",sprintf("%02d",c(12,1:11)),sep="_"))
|
|
20
|
|
|
21
|
dlst.m = cbind(dlst,dlst1,dlst2)
|
|
22
|
dlst.m$lstm=as.numeric(gsub("mm_","",dlst.m$variable))
|
|
23
|
|
|
24
|
dlst.m.sel=dlst.m[,c(1,2,3,5,10,15,16)]
|
|
25
|
|
|
26
|
dlst.all=dlst.m.sel[dlst.m.sel$month==dlst.m.sel$lstm,c("station","month","TMax","value","value.1","value.2")]
|
|
27
|
|
|
28
|
colnames(dlst.all)[colnames(dlst.all)=="value"]="lst"
|
|
29
|
colnames(dlst.all)[colnames(dlst.all)=="value.1"]="lstp1" # LST plus 1 month
|
|
30
|
colnames(dlst.all)[colnames(dlst.all)=="value.2"]="lstm1" # LST less 1 month
|
|
31
|
|
|
32
|
dlst.all$tile=t
|
|
33
|
dlst.full = rbind(dlst.all,dlst.full)
|
|
34
|
|
|
35
|
}
|
|
36
|
|
|
37
|
write.csv(dlst.full,file="/nobackup/awilso10/interp/lst_tmax.csv")
|
|
38
|
|
|
39
|
## copy to acrobates
|
|
40
|
system("scp /nobackup/awilso10/interp/lst_tmax.csv adamw@acrobates.eeb.yale.edu:/data/personal/adamw/projects/")
|
|
41
|
|
|
42
|
|
