Revision 9331c56d
Added by Benoit Parmentier almost 10 years ago
| climate/research/oregon/interpolation/global_run_scalingup_assessment_part2.R | ||
|---|---|---|
| 5 | 5 |
#Analyses, figures, tables and data are also produced in the script. |
| 6 | 6 |
#AUTHOR: Benoit Parmentier |
| 7 | 7 |
#CREATED ON: 03/23/2014 |
| 8 |
#MODIFIED ON: 11/17/2014
|
|
| 8 |
#MODIFIED ON: 11/30/2014
|
|
| 9 | 9 |
#Version: 3 |
| 10 | 10 |
#PROJECT: Environmental Layers project |
| 11 | 11 |
#COMMENTS: analyses for run 5 global using 6 specific tiles |
| ... | ... | |
| 215 | 215 |
#in_dir1 <- "/data/project/layers/commons/NEX_data/output_run6_global_analyses_09162014/output20Deg2" |
| 216 | 216 |
# parent output dir for the curent script analyes |
| 217 | 217 |
#out_dir <-"/data/project/layers/commons/NEX_data/output_run3_global_analyses_06192014/" #On NCEAS Atlas |
| 218 |
out_dir <-"/data/project/layers/commons/NEX_data/output_run9_global_analyses_11122014/"
|
|
| 218 |
out_dir <-"/data/project/layers/commons/NEX_data/output_run10_global_analyses_11302014/"
|
|
| 219 | 219 |
# input dir containing shapefiles defining tiles |
| 220 | 220 |
#in_dir_shp <- "/data/project/layers/commons/NEX_data/output_run5_global_analyses_08252014/output/subset/shapefiles" |
| 221 | 221 |
|
| ... | ... | |
| 228 | 228 |
|
| 229 | 229 |
y_var_name <- "dailyTmax" |
| 230 | 230 |
interpolation_method <- c("gam_CAI")
|
| 231 |
out_prefix<-"run9_global_analyses_11122014"
|
|
| 231 |
out_prefix<-"run10_global_analyses_11302014"
|
|
| 232 | 232 |
mosaic_plot <- FALSE |
| 233 | 233 |
|
| 234 |
#CRS_locs_WGS84 <- CRS("+proj=longlat +ellps=WGS84 +datum=WGS84 +towgs84=0,0,0") #Station coords WGS84
|
|
| 235 |
CRS_WGS84 <- CRS("+proj=longlat +ellps=WGS84 +datum=WGS84 +towgs84=0,0,0") #Station coords WGS84
|
|
| 236 |
|
|
| 234 | 237 |
proj_str<- CRS_WGS84 |
| 235 | 238 |
file_format <- ".rst" |
| 236 | 239 |
NA_value <- -9999 |
| ... | ... | |
| 259 | 262 |
#lf_tables <- list.files(out_dir,) |
| 260 | 263 |
summary_metrics_v <- read.table(file=file.path(out_dir,paste("summary_metrics_v2_NA_",out_prefix,".txt",sep="")),sep=",")
|
| 261 | 264 |
tb <- read.table(file=file.path(out_dir,paste("tb_diagnostic_v_NA","_",out_prefix,".txt",sep="")),sep=",")
|
| 265 |
#tb_diagnostic_s_NA_run10_global_analyses_11302014.txt |
|
| 262 | 266 |
tb_s <- read.table(file=file.path(out_dir,paste("tb_diagnostic_s_NA","_",out_prefix,".txt",sep="")),sep=",")
|
| 263 | 267 |
|
| 264 | 268 |
tb_month_s <- read.table(file=file.path(out_dir,paste("tb_month_diagnostic_s_NA","_",out_prefix,".txt",sep="")),sep=",")
|
Also available in: Unified diff
run 10 NEX part 2 assessment