#Part 1 create summary tables and inputs files for figure in part 2 and part 3.

#AUTHOR: Benoit Parmentier 

#CREATED ON: 03/23/2014  

#Version: 4

#PROJECT: Environmental Layers project  

#TO DO:

#in_dir1 <- "/nobackupp6/aguzman4/climateLayers/output1000x3000_km/"

#in_dir1 <- "/nobackupp6/aguzman4/climateLayers/out_15x45" #PARAM1

#in_dir1b <- "/nobackupp6/aguzman4/climateLayers/output1500x4500_km/singles" #PARAM1, add for now in_dir1 can be a list...

#region_names <- c("reg1","reg2","reg3","reg4","reg5","reg6") #selected region names, #PARAM2

#region_namesb <- c("reg_1b","reg_1c","reg_2b","reg_3b","reg_6b") #selected region names, #PARAM2

y_var_name <- "dailyTmax" #PARAM3

interpolation_method <- c("gam_CAI") #PARAM4

#output_run10_1500x4500_global_analyses_pred_2003_04102015/

#out_dir<-"/data/project/layers/commons/NEX_data/" #On NCEAS Atlas

CRS_locs_WGS84 <- CRS("+proj=longlat +ellps=WGS84 +datum=WGS84 +towgs84=0,0,0") #Station coords WGS84, #PARAM8

#day_to_mosaic <- c("20100101","20100901") #PARAM9

#day_to_mosaic <- NULL #if day to mosaic is null then mosaic all dates?

in_dir_list_all  <- lapply(in_dir_list,function(x){list.dirs(path=x,recursive=F)})

in_dir_list <- unlist(in_dir_list_all)

#in_dir_list <- in_dir_list[grep("bak",basename(basename(in_dir_list)),invert=TRUE)] #the first one is the in_dir1

in_dir_shp <- file.path(in_dir_subset,"shapefiles")

#select only directories used for predictions

#in_dir_reg <- in_dir_list[grep("july_tiffs",basename(in_dir_reg),invert=TRUE)] #select directory with shapefiles...

in_dir_list <- in_dir_reg

in_dir_list <- in_dir_list[grep("bak",basename(basename(in_dir_list)),invert=TRUE)] #the first one is the in_dir1

#list of shapefiles used to define tiles

in_dir_shp_list <- list.files(in_dir_shp,".shp",full.names=T)

##raster_prediction object : contains testing and training stations with RMSE and model object

list_raster_obj_files <- lapply(in_dir_list,FUN=function(x){list.files(path=x,pattern="^raster_prediction_obj.*.RData",full.names=T)})

basename(dirname(list_raster_obj_files[[1]]))

list_names_tile_coord <- lapply(list_raster_obj_files,FUN=function(x){basename(dirname(x))})