#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/output1500x4500_km/" #PARAM1, add for now in_dir1 can be a list...

#in_dir1 <- "/nobackupp6/aguzman4/climateLayers/output1500x4500_km/elevTest/1kmBuff/"

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

#region_names <- c("reg4") #selected region names, #PARAM2

#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/

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

#day_to_mosaic <- c("20100829","20100830","20100831",

#                   "20100901","20100902","20100903")

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

#recombine region first:

region_names_mosaic <- list(region_names)

in_dir_mosaics <- lapply(region_names_mosaic,FUN=function(x){file.path(in_dir1,x)})

for (j in 1:length(region_names_mosaic)){