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

#AUTHOR: Benoit Parmentier 

#CREATED ON: 03/23/2014  

#Version: 3

#PROJECT: Environmental Layers project  

#TO DO:

#### Parameters and constants  

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

#/nobackupp4/aguzman4/climateLayers/output10Deg/reg1/finished.txt

in_dir_list <- list.dirs(path=in_dir1,recursive=FALSE) #get the list regions processed for this run

#if(basename(in_dir_list)[[1]]=="reg?") #add later

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

# the last directory contains shapefiles 

y_var_name <- "dailyTmax"

interpolation_method <- c("gam_CAI")

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

out_dir <- "/nobackup/bparmen1/" #on NEX

df_tile_processed$path_NEX <- in_dir_list

##Quick exploration of raster object

#robj1 <- load_obj(lf_raster_obj[4]) #This is tile tile

names(robj1)

#for i in 1:length(df_tiled_processed$tile_coord)

#output_atlas_dir <- "/data/project/layers/commons/NEX_data/output_run3_global_analyses_06192014/output10Deg/reg1"

#output_atlas_dir <- "/data/project/layers/commons/NEX_data/output_run5_global_analyses_08252014/output20Deg"

#Make directories on ATLAS

#for (i in 1:length(df_tile_processed$tile_coord)){

#  create_dir_fun(file.path(output_atlas_dir,as.character(df_tile_processed$tile_coord[i])),out_suffix=NULL)

  system(cmd_str)

}

##################### END OF SCRIPT ######################