#COMMIT: generating animation for reg6 (Australia and South East Asia)

function_product_assessment_part2_functions <- "global_product_assessment_part2_functions_10102016b.R"

in_dir <- "/data/project/layers/commons/NEX_data/climateLayers/out/reg6/assessment"

in_dir_mosaic <- "/data/project/layers/commons/NEX_data/climateLayers/out/reg6/mosaics/mosaic" #predicted mosaic

region_name <- c("reg6") #param 6, arg 3

out_suffix <- "global_assessment_reg6_10102016"

out_dir <- "/data/project/layers/commons/NEX_data/climateLayers/out/reg6/assessment"

infile_mask <- "/data/project/layers/commons/NEX_data/regions_input_files/r_mask_LST_reg6.tif"

df_centroids_fname <- "/data/project/layers/commons/NEX_data/climateLayers/out/reg6/mosaic/output_reg6_1984/df_centroids_19840101_reg6_1984.txt"

lf_raster <- NULL #list of raster to consider

#plot(r_stack,zlim=c(-50,50),col=matlab.like(255))

#function_product_assessment_part2_functions <- "global_product_assessment_part2_functions_10102016.R"

#source(file.path(script_path,function_product_assessment_part2_functions)) #source all functions used in this script 

### Now run for the full time series

#13.26 Western time: start

zlim_val <- NULL

list_param_plot_raster_mosaic <- list(l_dates,r_stack_subset,NA_flag_val,out_dir,out_suffix,

                                      region_name,variable_name, zlim_val)

names(list_param_plot_raster_mosaic) <- c("l_dates","r_mosaiced_scaled","NA_flag_val_mosaic","out_dir","out_suffix",

                                          "region_name","variable_name","zlim_val")

lf_mosaic_plot_fig <- mclapply(1:length(l_dates[1:11]),

                               FUN=plot_raster_mosaic,

                               list_param=list_param_plot_raster_mosaic,

                               mc.preschedule=FALSE,

                               mc.cores = num_cores)  

##start at 14.12

##finished at 16.47

#r_stack_subset <- subset(r_stack,1:11)

#l_dates <- list_dates_produced_date_val[1:11]

filenames_figures_mosaic_test <- "list_figures_animation_test_reg6.txt"

write.table(unlist(lf_mosaic_plot_fig[1:11]),filenames_figures_mosaic_test,row.names = F,col.names = F,quote = F)

filenames_figures_mosaic <- paste0("list_figures_animation_",out_suffix_movie,".txt")

write.table(unlist(lf_mosaic_plot_fig),filenames_figures_mosaic,row.names = F,col.names = F,quote = F)

#started

generate_animation_from_figures_fun(filenames_figures= filenames_figures_mosaic_test,

                                    out_filename_figure_animation="test_reg6_animation.gif")

generate_animation_from_figures_fun(filenames_figures= unlist(lf_mosaic_plot_fig[1:11]),

                                    out_filename_figure_animation="test2_reg6_animation.gif")

#started 17.36 Western time on Oct 10 and 18.18

generate_animation_from_figures_fun(filenames_figures= filenames_figures_mosaic,

                                    out_suffix=out_suffix_movie,

                                    out_filename_figure_animation=NULL)

############ Now accuracy

pattern_str <-"*.tif"

in_dir_mosaic <- in_dir_mosaic_RMSE

lf_raster_rmse <- list.files(path=in_dir_mosaic,pattern=pattern_str,recursive=F,full.names=T)

r_stack <- stack(lf_raster,quick=T) #this is very fast now with the quick option!

  #save(r_mosaic,file="r_mosaic.RData")

NAvalue(r_stack)

plot(r_stack,y=6,zlim=c(-10000,10000)) #this is not rescaled

#plot(r_stack,zlim=c(-50,50),col=matlab.like(255))

#plot(r_mosaic_scaled,y=6,zlim=c(-50,50))

#plot(r_mosaic_scaled,zlim=c(-50,50),col=matlab.like(255))

#debug(extract_date)

#test <- extract_date(6431,lf_mosaic_list,12) #extract item number 12 from the name of files to get the data

#list_dates_produced <- unlist(mclapply(1:length(lf_raster),FUN=extract_date,x=lf_raster,item_no=13,mc.preschedule=FALSE,mc.cores = num_cores))                         

lf_mosaic_list <- lf_raster

list_dates_produced <-  mclapply(1:2,

                                 FUN=extract_date,

                                 x=lf_mosaic_list,

                                 item_no=13,

                                 mc.preschedule=FALSE,

                                 mc.cores = 2)  

item_no <-13

list_dates_produced <- unlist(mclapply(1:length(lf_raster),

                                       FUN=extract_date,

                                       x=lf_raster,

                                       item_no=item_no,

                                       mc.preschedule=FALSE,

                                       mc.cores = num_cores))                         

list_dates_produced_date_val <- as.Date(strptime(list_dates_produced,"%Y%m%d"))

month_str <- format(list_dates_produced_date_val, "%b") ## Month, char, abbreviated

year_str <- format(list_dates_produced_date_val, "%Y") ## Year with century

day_str <- as.numeric(format(list_dates_produced_date_val, "%d")) ## numeric month

df_raster <- data.frame(lf=basename(lf_raster),

                          date=list_dates_produced_date_val,

                          month_str=month_str,

                          year=year_str,

                          day=day_str,

                          dir=dirname(lf_mosaic_list))

df_raster_fname <- file.path(out_dir,paste0("df_raster_",out_suffix,".txt"))

write.table(df_raster,file= df_raster_fname,sep=",",row.names = F)