#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

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

source(file.path(script_path,"GAM_fusion_function_multisampling_07022013.R")) #Include GAM_CAI

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

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

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

prop_minmax<-c(0.3,0.3)  #if prop_min=prop_max and step=0 then predicitons are done for the number of dates...

#dates_selected<-c("20100101","20100102","20100103","20100901") # Note that the dates set must have a specific format: yyymmdd

dates_selected<-"" # if empty string then predict for the full year specified earlier

#Models to run...this can be change for each run

lst_avg<-c("mm_01","mm_02","mm_03","mm_04","mm_05","mm_06","mm_07","mm_08","mm_09","mm_10","mm_11","mm_12")  

#Collect all parameters in a list

                                seed_number,nb_sample,step,constant,prop_minmax,dates_selected,

                                list_models,lst_avg,out_path,script_path,

                                interpolation_method)

                                "seed_number","nb_sample","step","constant","prop_minmax","dates_selected",

                                "list_models","lst_avg","out_path","script_path",

                                "interpolation_method")