#STAGE 8: Comparison of predictions across regions and years with figures generation.

#AUTHOR: Benoit Parmentier                                                                        

#CREATED ON: 12/29/2015  

#PROJECT: NCEAS-IPLANT-NASA: Environment Layers                                                                           

#First source these files:

#module load pythonkits/gdal_1.10.0_python_2.7.3_nex

## TODO:

# Make this callable from the shell

# Adapt for python script

# Fix figure for part2

# Call also use library(optparse)

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

#

### PARAMETERS DEFINED IN THE SCRIPT

#master directory containing the definition of tile size and tiles predicted

#in_dir1 <- "/nobackupp6/aguzman4/climateLayers/out/" #param 5, arg 2

in_dir <- args[2] #param 5, arg 2

#in_dir <- "" #PARAM 0

#y_var_name <- "dailyTmax" #PARAM1

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

#run assessment by region, this is a unique region only 

region_name <- args[3] #param 6, arg 3

#out_prefix <- "run_global_analyses_pred_12282015" #param 7, arg 4

#out_dir <- "/nobackupp8/bparmen1/" #param 8, arg 5

#out_dir <-paste(out_dir,"_",out_prefix,sep="")

create_out_dir_param <- TRUE #param 9, arg 6

out_prefix <- args[4] #param 7, arg 4

out_dir <- args[5] #param 8, arg 5

#out_dir <-paste(out_dir,"_",out_prefix,sep="")

create_out_dir_param <- args[6] #param 9, arg 6

max_mem <- args[9] #param 21

max_mem <- 1e+07

#in_dir_list_filename <- args[10] #param 22

#in_dir_list_filename <- "/data/project/layers/commons/NEX_data/reg4_assessment/stage6_reg4_in_dir_list_02072016.txt"

run_figure_by_year <- args[11] # param 10, arg 7, if true will create figures for individual years...

run_figure_by_year <- TRUE # param 10, arg 7, if true will create figures for individual years...