#STAGE 5: Output analyses: assessment of results for specific dates...

#

#AUTHOR: Benoit Parmentier                                                                       

#PROJECT: NCEAS INPLANT: Environment and Organisms --TASK#363, TASK$568--   

clim_script <- file.path(script_path,"climatology_05312013.py") # LST climatology python script

grass_setting_script <- file.path(script_path,"grass-setup.R") #Set up system shell environment for python+GRASS

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

source(file.path(script_path,"Database_stations_covariates_processing_function_06112013.R"))

source(file.path(script_path,"results_interpolation_date_output_analyses_06112013.R"))

#source(file.path(script_path,"results_covariates_database_stations_output_analyses_04012013.R")) #to be completed

#FUNCTIONS CALLED FROM GAM ANALYSIS RASTER PREDICTION ARE FOUND IN...

source(file.path(script_path,"sampling_script_functions_03122013.R"))

source(file.path(script_path,"interpolation_method_day_function_multisampling_06082013.R")) #Include GAM_day

source(file.path(script_path,"GAM_fusion_function_multisampling_validation_metrics_05062013.R"))

stages_to_run<-c(0,2,3,4,5) #May decide on antoher strategy later on...

var<-"TMAX" # variable being interpolated

out_suffix_modis <-"_05302013" #use tiles produce previously

#interpolation_method<-c("gam_fusion","gam_CAI","gam_daily") #other otpions to be added later

#interpolation_method<-c("gam_CAI") #other otpions to be added later

#interpolation_method<-c("kriging_daily") #other otpions to be added later

#interpolation_method<-c("gwr_daily") #other otpions to be added later