#         Using python scripts started by J. Regetz and modified by B.Parmentier.

#DATE: 05/30/2013                                                                                 

#To Do:

#-collect list of LST files created

#-solve issue realted to parallelization of python script: GRASS problem?

  #DATE: 05/30/2013                                                                                 

  list_tiles_modis <- unlist(strsplit(list_tiles_modis,","))  # transform string into separate element in char vector

  list_tiles_modis <- list_tiles_modis[j]

  command_download_str <- paste("python",modis_download_script, list_param_python_script_str,sep=" ")

  command_clim_str <- paste("python",clim_script, list_param_python_script_str,sep=" ") #preparing shell command

    source(grass_setting_script) #Setting to access GRASS (current setting valid for Atlas at nceas)  

#   ## return list of files???, may be modified later to return clim list and list of downloaded files + missing.

#   if (clim_calc==1){

#     list_output_obj <-list(command_download_str,command_clim_str)

#     names(list_output_obj) <- c("command_download_str", "command_clim_str")

#   } else{

#     list_output_obj <-list(command_download_str)

#     names(list_output_obj) <- c("command_download_str")

#   }

  #Add function to list files created or modify in python...

  list_output_obj <-list(command_download_str,command_clim_str)

  names(list_output_obj) <- c("command_download_str", "command_clim_str")

#list_tiles_modis <- c("h12v04,h13v04,h14v04,h07v06") #tiles for Northenr America, Northern US...

#list_tiles_modis <- c("h09v04,h09v04") #tiles for Oregon #defined above...

#script_path<-"/home/parmentier/Data/IPLANT_project/env_layers_scripts/"

#modis_download_script <- file.path(script_path,"modis_download_05142013.py") # LST modis download python script

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

#grass_setting_script <- file.path(script_path,"grass-setup.R")

#var="TMAX"

#start_year = "2001"

#end_year = "2010"

#end_year = "2002"

#hdfdir =  '/data/project/layers/commons/modis/MOD11A1_tiles' #destination file where hdf files are stored locally after download.

#download=0

#clim_calc=1

#out_suffix_modis="_05302013"

#list_param_download_clim_LST_script <- list(list_tiles_modis,start_year,end_year,hdfdir,

#                                 var,grass_setting_script,modis_download_script, clim_script,

#                                 download,clim_calc,out_suffix_modis)

#names(list_param_download_clim_LST_script)<-c("list_tiles_modis","start_year","end_year","hdfdir",

#                                   "var","grass_setting_script","modis_download_script","clim_script",

#                                   "download","clim_calc","out_suffix_modis")

#clim_production_obj <-lapply(1:2, list_param=list_param_download_clim_LST_script, download_calculate_MODIS_LST_climatology) #,mc.preschedule=FALSE,mc.cores = 2) #This is the end bracket from mclapply(...) statement

#download_calculate_MODIS_LST_climatology(1,list_param_download_clim_LST_script)

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