#MODIFIED ON: 02/16/2015            

out_prefix <- "run10_1000x3000_global_analyses_02162015"

out_dir <- "/data/project/layers/commons/NEX_data/output_run10_1000x3000_global_analyses_02162015"

#tb <- subset(tb,reg!="reg3b")

#summary_metrics_v <- subset(summary_metrics_v,reg!="reg3b")

#tb_s_all <- tb_s

#tb_s_all <- subset(tb_s,reg!="reg3b")

#tb_month_s_all <- tb_month_s

#tb_month_s <- subset(tb_month_s,reg!="reg3b")

#df_tile_processed_all <- df_tile_processed

#df_tile_processed <- subset(df_tile_processed,reg!="reg3b")

sum(df_tile_processed$metrics_v) #188,number of tiles with raster object

length(df_tile_processed$metrics_v) #214,number of tiles in the region

sum(df_tile_processed$metrics_v)/length(df_tile_processed$metrics_v) #87.85% of tiles with info

#l_reg_name <- unique(df_tile_processed$reg)

#lf_mosaics_reg <- vector("list",length=length(l_reg_name))

#for(i in 1:length(l_reg_name)){

#    lf_mosaics_reg[[i]] <- try(mixedsort(

#    list.files(

#    path=file.path(out_dir,"mosaics"),

#    #pattern="reg6_.*._CAI_TMAX_clim_month_.*._mod1_all_mean.tif",

#    pattern=paste(l_reg_name[i],".*._CAI_TMAX_clim_month_.*._mod1_all_mean.tif",sep=""), 

#    full.names=T))

#  )

#}

#names(lf_mosaics_reg) <- l_reg_name

lf_raster_fname <- list.files(pattern="world.*2010*.*02162015.tif$",full.names=T)

prefix_str <- "Figure10_clim_world_mosaics_day_"

l_dates <-day_to_mosaic

screenRast=TRUE

list_param_plot_screen_raster <- list(lf_raster_fname,screenRast,l_dates,out_dir,prefix_str,out_suffix)

names(list_param_plot_screen_raster) <- c("lf_raster_fname","screenRast","l_dates","out_dir_str","prefix_str","out_suffix_str")

#debug(plot_screen_raster_val)

#plot_screen_raster_val(1,list_param_plot_scre

en_raster)

world_m_list <- mclapply(1:10, list_param=list_param_plot_screen_raster, plot_screen_raster_val,mc.preschedule=FALSE,mc.cores = 5) #This is the end bracket from mclapply(...) statement

world_m_list <- mclapply(11:30, list_param=list_param_plot_screen_raster, plot_screen_raster_val,mc.preschedule=FALSE,mc.cores = 7) #This is the end bracket from mclapply(...) statement

plot_screen_raster_val<-function(i,list_param){

  ##USAGE###

  #Screen raster list and produced plot as png.

  fname <-list_param$lf_raster_fname[i]

  screenRast <- list_param$screenRast

  l_dates<- list_param$l_dates

  out_dir_str <- list_param$out_dir_str

  prefix_str <-list_param$prefix_str

  out_suffix_str <- list_param$out_suffix_str

  ### START SCRIPT ####

  date_proc <- l_dates[i]

  if(screenRast==TRUE){

    r_test <- raster(fname)

    m <- c(-Inf, -100, NA,  

    rclmat <- matrix(m, ncol=3, byrow=TRUE)

    rc <- reclassify(r_test, rclmat,filename=paste("rc_tmp_",i,".tif",sep=""),dataType="FLT4S",overwrite=T)

    #file_name <- unlist(strsplit(basename(lf_m[i]),"_"))

    extension_str <- extension(filename(r_test))

    raster_name_tmp <- gsub(extension_str,"",basename(filename(r_test)))

    raster_name <- file.path(out_dir_str,paste(raster_name_tmp,"_masked.tif",sep=""))

    r_pred <- mask(r_test,rc,filename=raster_name,overwrite=TRUE)

  }else{

    r_pred <- raster(fname)

  }

  #date_proc <- file_name[7] #specific tot he current naming of files

  #date_proc<- "2010010101"

  #paste(raster_name[1:7],collapse="_")

  #add filename option later

  res_pix <- 960

  #res_pix <- 480

  col_mfrow <- 1

  row_mfrow <- 1

#  png(filename=paste("Figure10_clim_world_mosaics_day_","_",date_proc,"_",tile_size,"_",out_suffix,".png",sep=""),

#    width=col_mfrow*res_pix,height=row_mfrow*res_pix)

  png(filename=paste(prefix_str,"_",date_proc,"_",tile_size,"_",out_suffix_str,".png",sep=""),

  dev.off()

}