#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

#in_dir_list <- in_dir_list[c(3,4)] #get the list regions processed for this run

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

#in_dir_list_all <- in_dir_list

#in_dir_list <- list.dirs(path=in_dir_reg,recursive=FALSE) #get the list of tiles/directories with outputs 

#in_dir_list <- unlist(in_dir_list_all[c(2)]) #only region 3 has informatation at this stage

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

in_dir_subset <- in_dir_list[grep("subset",basename(in_dir_list),invert=FALSE)] #select directory with shapefiles...

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

}

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

#Analyses, figures, tables and data are also produced in the script.

#AUTHOR: Benoit Parmentier 

#CREATED ON: 03/23/2014  

#Version: 3

#PROJECT: Environmental Layers project     

#COMMENTS: analyses for run 5 global using 6 specific tiles

#in_dir1 <- "/data/project/layers/commons/NEX_data/output_run6_global_analyses_09162014/output20Deg2"

# parent output dir for the curent script analyes

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

# input dir containing shapefiles defining tiles

#in_dir_shp <- "/data/project/layers/commons/NEX_data/output_run5_global_analyses_08252014/output/subset/shapefiles"

y_var_name <- "dailyTmax"

interpolation_method <- c("gam_CAI")

mosaic_plot <- FALSE

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

#gam_diagnostic_df <- read.table(file=file.path(out_dir,"gam_diagnostic_df_run4_global_analyses_08142014.txt"),sep=",")

########################## START SCRIPT ##############################

#multiple regions?

if(mulitple_region==TRUE){

}

#This part is specifically related to this run : dropping model with elev

#ac_mod <- summary_metrics_v[summary_metrics_v$pred_mod==model_name[i],]

###############

### Figure 1: plot location of the study area with tiles processed

  #plot(ac_mod1,cex=sqrt(ac_mod1$rmse),pch=1,add=T)

  #plot(ac_mod,cex=(ac_mod$rmse^2)/10,pch=1,col="red",add=T)

  #coordinates(ac_mod) <- ac_mod[,c("lon.x","lat.x")] #solve this later

  p_shp <- layer(sp.polygons(reg_layer, lwd=1, col='black'))

  #title("(a) Mean for 1 January")

sum(df_tile_processed$metrics_v)/length(df_tile_processed$metrics_v) #81.40%

#coordinates

summary_metrics_v$n_missing <- summary_metrics_v$n == 365

summary_metrics_v$n_missing <- summary_metrics_v$n < 365