script_path<-list_param_raster_prediction$script_path

  interpolation_method<-list_param_raster_prediction$interpolation_method

  #Sourcing in the master script to avoid confusion on the latest versions of scripts and functions!!!

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

      file=log_fname,sep="\n")

  t1<-proc.time()

  j=12

  if (interpolation_method=="gam_fusion"){

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

    clim_method_mod_obj<-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

    #test<-runClim_KGFusion(1,list_param=list_param_runClim_KGFusion)

    #j indicate which month is predicted, var indicates TMIN or TMAX

    data_name<-paste(var,"_bias_LST_month_",j,"_",cname[k],"_",prop_month,

                     "_",run_samp,sep="")

    list_out_filename[[k]]<-raster_name

  }