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

#AUTHOR: Benoit Parmentier 

#CREATED ON: 04/14/2015  

#Version: 5

#PROJECT: Environmental Layers project     

#COMMENTS: analyses run for reg4 1992 for test of mosaicing using 1500x4500km and other tiles

#### FUNCTION USED IN SCRIPT

#in_dir_script <-"/home/parmentier/Data/IPLANT_project/env_layers_scripts" #NCEAS UCSB

in_dir_script <- "/nobackupp8/bparmen1/env_layers_scripts" #NASA NEX

algorithm <- "python" #if R use mosaic function for R, if python use modified gdalmerge script from Alberto Guzmann

#algorithm <- "R" #if R use mosaic function for R, if python use modified gdalmerge script from Alberto Guzmann

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

out_suffix_str <- out_suffix

#Improve by adding multicores option

list_param_accuracy_metric_raster <- list(lf,tb,metric_name,pred_mod_name,y_var_name,interpolation_method,

names(list_param_accuracy_metric_raster) <- c("lf","tb","metric_name","pred_mod_name","y_var_name","interpolation_method",

list_raster_created_obj <- lapply(1:length(day_to_mosaic),FUN=create_accuracy_metric_raster,

                                  list_param=list_param_accuracy_metric_raster)

#debug(create_accuracy_metric_raster)

  mosaic_method <- "use_edge_weights" #this is distance from edge

  out_suffix_tmp <- paste(interpolation_method,y_var_name,day_to_mosaic[i],out_suffix,sep="_")

  #can also loop through methods!!!

  #python_bin <- "/usr/bin/" #python gdal bin, on Atlas NCEAS

  #python_bin <- "/nobackupp6/aguzman4/climateLayers/sharedModules/bin" #on NEX

                               out_dir_str=out_dir,

                               out_prefix_str <- "dailyTmax_",

                               out_suffix_str=out_suffix)

l_png_files <- mclapply(1:length(unlist(lf_mean_mosaic)[1:2]),FUN=plot_screen_raster_val,

                        list_param= list_param_plot_mosaic,

                        mc.preschedule=FALSE,mc.cores = num_cores)