#

#AUTHOR: Benoit Parmentier 

#CREATED ON: 10/27/2016  

#Version: 1

#PROJECT: Environmental Layers project     

#COMMENTS: 

###### Function used in the script #######

script_path <- "/nobackupp8/bparmen1/env_layers_scripts" #path to script, NEX

#script_path <- "/home/parmentier/Data/IPLANT_project/env_layers_scripts" #path to script

source(file.path(script_path,function_assessment_part3)) #source all functions used in this script 

#Product assessment

source(file.path(script_path,function_product_assessment_part0_functions)) #source all functions used in this script 

##Don't load part 1 and part2, mosaic package does not work on NEX

#function_product_assessment_part1_functions <- "global_product_assessment_part1_functions_09192016b.R"

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

####### Parameters, constants and arguments ###

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

proj_str <- "+proj=longlat +ellps=WGS84 +datum=WGS84 +towgs84=0,0,0"

region_name <- c("reg6") #param 6, arg 3

out_suffix <- "global_assessment_reg6_10232016"

create_out_dir_param <- TRUE #param 9, arg 6

file_format <- ".tif" #format for mosaiced files # param 11

NA_flag_val <- -32768  #No data value, # param 12

plotting_figures <- TRUE #running part2 of assessment to generate figures... # param 14

NA_flag_val_mosaic <- -32768

in_dir_list_filename <- NULL #if NULL, use the in_dir directory to search for info

#countries_shp <-"/data/project/layers/commons/NEX_data/countries.shp" #Atlas

lf_raster <- NULL #list of raster to consider

#On NEX

#contains all data from the run by Alberto

NA_flag_val <- NA_value

#multiple_region <- TRUE #PARAM 12

region_name <- "reg6" #PARAM 13

#/nobackupp6/aguzman4/climateLayers/out/reg6/subset/shapefiles

list_year_predicted <- c(2000,2012,2013) #year still on disk for reg6

##### prepare list of parameters for call of function

source(file.path(script_path,function_product_assessment_part0_functions)) #source all functions used in this script 

list_param_predictions_tiles_missing <- list(in_dir1,region_name,y_var_name,interpolation_method,out_suffix,out_dir,

                                             "pred_mod_name","metric_name")

#debug(predictions_tiles_missing_fun)

predictions_tiles_missing_fun(1,list_param=list_param_predictions_tiles_missing)

#### Function to check missing tiles and estimate potential gaps