Revision f07e0e66
Added by Adam Wilson almost 11 years ago
| climate/procedures/ee_compile.R | ||
|---|---|---|
| 106 | 106 |
|
| 107 | 107 |
### merge all the tiles to a single global composite |
| 108 | 108 |
#system(paste("ncdump -h ",list.files(tempdir,pattern="mod09.*.nc$",full=T)[10]))
|
| 109 |
system(paste("cdo -O mergetime -setctomiss,-32768 ",paste(list.files(tempdir,pattern="mod09.*.nc$",full=T),collapse=" ")," data/cloud_monthly.nc"))
|
|
| 109 |
file.remove("tmp/mod09_2000-01-15.nc")
|
|
| 110 |
system(paste("cdo -O mergetime -setrtomiss,-32768,-1 ",paste(list.files(tempdir,pattern="mod09.*.nc$",full=T),collapse=" ")," data/cloud_monthly.nc"))
|
|
| 110 | 111 |
|
| 111 | 112 |
# Overall mean |
| 112 |
system(paste("cdo -O -chname,CF,CF_annual -timmean data/cloud_monthly.nc data/cloud_mean.nc"))
|
|
| 113 |
system(paste("cdo -O timmean data/cloud_monthly.nc data/cloud_mean.nc"))
|
|
| 113 | 114 |
|
| 114 | 115 |
### generate the monthly mean and sd |
| 115 | 116 |
#system(paste("cdo -P 10 -O merge -ymonmean data/mod09.nc -chname,CF,CF_sd -ymonstd data/mod09.nc data/mod09_clim.nc"))
|
| 116 | 117 |
system(paste("cdo -f nc4c -O -ymonmean data/cloud_monthly.nc data/cloud_ymonmean.nc"))
|
| 118 |
|
|
| 119 |
|
|
| 120 |
## Seasonal Means |
|
| 121 |
system(paste("cdo -f nc4c -O -yseasmean data/cloud_monthly.nc data/cloud_yseasmean.nc"))
|
|
| 122 |
system(paste("cdo -f nc4c -O -yseasstd data/cloud_monthly.nc data/cloud_yseasstd.nc"))
|
|
| 123 |
|
|
| 117 | 124 |
## standard deviations, had to break to limit memory usage |
| 118 | 125 |
system(paste("cdo -f nc4c -O -chname,CF,CF_sd -ymonstd -selmon,1,2,3,4,5,6 data/cloud_monthly.nc data/cloud_ymonsd_1-6.nc"))
|
| 119 | 126 |
system(paste("cdo -f nc4c -O -chname,CF,CF_sd -ymonstd -selmon,7,8,9,10,11,12 data/cloud_monthly.nc data/cloud_ymonsd_7-12.nc"))
|
| ... | ... | |
| 121 | 128 |
|
| 122 | 129 |
|
| 123 | 130 |
#if(!file.exists("data/mod09_metrics.nc")) {
|
| 124 |
system("cdo -f nc4c -chname,CF,CFmin -timmin data/cloud_ymonmean.nc data/cloud_min.nc")
|
|
| 125 |
system("cdo -f nc4c -chname,CF,CFmax -timmax data/cloud_ymonmean.nc data/cloud_max.nc")
|
|
| 131 |
system("cdo -f nc4c timmin data/cloud_ymonmean.nc data/cloud_min.nc")
|
|
| 132 |
system("cdo -f nc4c timmax data/cloud_ymonmean.nc data/cloud_max.nc")
|
|
| 126 | 133 |
system("cdo -f nc4c -chname,CF,CFsd -timstd data/cloud_ymonmean.nc data/cloud_std.nc")
|
| 127 | 134 |
# system("cdo -f nc2 merge data/mod09_std.nc data/mod09_min.nc data/cloud_max.nc data/cloud_metrics.nc")
|
| 128 | 135 |
# system("cdo merge -chname,CF,CFmin -timmin data/cloud_ymonmean.nc -chname,CF,CFmax -timmax data/cloud_ymonmean.nc -chname,CF,CFsd -timstd data/cloud_ymonmean.nc data/cloud_metrics.nc")
|
| 129 | 136 |
#} |
| 130 | 137 |
|
| 131 | 138 |
|
| 139 |
# Regressions through time by season |
|
| 140 |
s=c("DJF","MAM","JJA","SON")
|
|
| 141 |
|
|
| 142 |
system(paste("cdo -f nc4c -O regres -selseas,",s[1]," data/cloud_monthly.nc data/slope_",s[1],".nc &",sep=""))
|
|
| 143 |
system(paste("cdo -f nc4c -O regres -selseas,",s[2]," data/cloud_monthly.nc data/slope_",s[2],".nc &",sep=""))
|
|
| 144 |
system(paste("cdo -f nc4c -O regres -selseas,",s[3]," data/cloud_monthly.nc data/slope_",s[3],".nc &",sep=""))
|
|
| 145 |
system(paste("cdo -f nc4c -O regres -selseas,",s[4]," data/cloud_monthly.nc data/slope_",s[4],".nc &",sep=""))
|
|
| 132 | 146 |
|
| 133 | 147 |
|
| 134 | 148 |
|
Also available in: Unified diff
udpates to validation and addition of initial biodiversity script