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

#

#AUTHOR: Benoit Parmentier                                                                       

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

##SCRIPT USED FOR THE PREDICTIONS: Source or list all scripts here to avoid confusion on versions being run!!!!

#CALLED FROM MASTER SCRIPT:

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,"GAM_fusion_analysis_raster_prediction_multisampling_07302013.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...

met_stations_outfiles_obj_file<-"/data/project/layers/commons/data_workflow/output_data_365d_gam_fus_lst_test_run_07172013/met_stations_outfiles_obj_gam_fusion__365d_gam_fus_lst_test_run_07172013.RData"

var<-"TMAX" # variable being interpolated

out_suffix_modis <-"_05302013"                       #pattern to find tiles produced 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_fusion") #other otpions to be added later

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

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

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

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

lc_path<-"/data/project/layers/commons/data_workflow/inputs/lc-consensus-global"

infile_modis_grid<-"/data/project/layers/commons/data_workflow/inputs/modis_grid/modis_sinusoidal_grid_world.shp" #modis grid tiling system, global

infile_elev<-"/data/project/layers/commons/data_workflow/inputs/dem-cgiar-srtm-1km-tif/srtm_1km.tif"  #elevation at 1km, global extent to be replaced by the new fused product 

infile_distoc <- "/data/project/layers/commons/data_workflow/inputs/distance_to_coast/GMT_intermediate_coast_distance_01d_rev.tif" #distance to coast, global extent at 0.01 deg

#infile_covariates<- "/home/parmentier/Data/IPLANT_project/Venezuela_interpolation/Venezuela_01142013/covariates_Oregon_region_TMAX__OR_04052013.tif" #Oregon covar TMAX from earlier codes...for continuity

#infile_reg_outline=""  #input region outline defined by polygon: none for Venezuela

#LC1: Evergreen/deciduous needleleaf trees

#Combination 3: for paper baseline=s(lat,lon)+s(elev)

#Default name of LST avg to be matched               

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")  

############## STAGE 5: OUTPUT ANALYSES ##################

date_selected_results<-c("20100101") 

list_param_results_analyses<-list(out_path,script_path,raster_prediction_obj,interpolation_method,

names(list_param_results_analyses)<-c("out_path","script_path","raster_prediction_obj","interpolation_method",

#plots_assessment_by_date<-function(j,list_param){

if (stages_to_run[5]==5){

  #Use lapply or mclapply

  summary_v_day <-plots_assessment_by_date(1,list_param_results_analyses)

  #Call as function...