#source(file.path(script_path,"GAM_fusion_function_multisampling_03122013.R"))

  gamclim_fus_mod<-mclapply(1:12, list_param=list_param_runClim_KGFusion, runClim_KGFusion,mc.preschedule=FALSE,mc.cores = 6) #This is the end bracket from mclapply(...) statement

  #gamclim_fus_mod<-mclapply(1:6, runClim_KGFusion,mc.preschedule=FALSE,mc.cores = 6) #This is the end bracket from mclapply(...) statement

  t2<-proc.time()-t1

  writeLines(as.character(t2),con=log_file,sep="\n")

  #gam_fus_mod<-mclapply(1:length(sampling_obj$ghcn_data_day),runGAMFusion,list_param_runGAMFusion,mc.preschedule=FALSE,mc.cores = 9) #This is the end bracket from mclapply(...) statement

  #gam_fus_mod<-mclapply(1:length(ghcn.subsets), runGAMFusion,mc.preschedule=FALSE,mc.cores = 9) #This is the end bracket from mclapply(...) statement

  t2<-proc.time()-t1

  writeLines(as.character(t2),con=log_file,sep="\n")

  #browser()

  gam_fus_validation_mod<-mclapply(1:length(gam_fus_mod), list_param=list_param_validation, calculate_accuracy_metrics,mc.preschedule=FALSE,mc.cores = 9) #This is the end bracket from mclapply(...) statement

  #gam_fus_validation_mod<-mclapply(1:1, calculate_accuracy_metrics,mc.preschedule=FALSE,mc.cores = 1) #This is the end bracket from mclapply(...) statement

  t2<-proc.time()-t1

  writeLines(as.character(t2),con=log_file,sep="\n")

  #Call function to create plots of metrics for validation dataset

  metric_names<-c("rmse","mae","me","r","m50")

  #################### CLOSE LOG FILE  ####################

  ################### PREPARE RETURN OBJECT ###############

  #Will add more information to be returned

  names(raster_prediction_obj)<-c("gamclim_fus_mod","gam_fus_mod","gam_fus_validation_mod","tb_diagnostic_v",

                                  "summary_metrics_v")  

  return(raster_prediction_obj)

}