Project

General

Profile

« Previous | Next » 

Revision 86a8f1c7

Added by Benoit Parmentier almost 9 years ago

assessment part3, adding tile and dates list from region assessment

View differences:

climate/research/oregon/interpolation/global_run_scalingup_assessment_part3.R
6 6 
#Analyses, figures, tables and data are also produced in the script.
7 7 
#AUTHOR: Benoit Parmentier
8 8 
#CREATED ON: 03/23/2014
9 
#MODIFIED ON: 02/02/2016
9 
#MODIFIED ON: 02/05/2016
10 10 
#Version: 5
11 11 
#PROJECT: Environmental Layers project
12 12 
#COMMENTS: Initial commit, script based on part 2 of assessment, will be modified further for overall assessment
......
47 47 
#parent output dir for the current script analyes
48 48 
#out_dir <- "/nobackup/bparmen1/" #on NEX
49 49 
#in_dir_shp <- "/nobackupp4/aguzman4/climateLayers/output4/subset/shapefiles/"
50 

  		




  
  51
  
    
in_dir <- "/data/project/layers/commons/NEX_data"


  		




  
  52
  
    
#list_in_dir_run <-


  		




  
  53
  
    
#in_dir_list <-  c("output_run_global_analyses_pred_2009_reg4","output_run_global_analyses_pred_2010_reg4",


  		




  
  54
  
    
#                  "output_run_global_analyses_pred_2011_reg4","output_run_global_analyses_pred_2012_reg4",


  		




  
  55
  
    
#                  "output_run_global_analyses_pred_2013_reg4","output_run_global_analyses_pred_2014_reg4")


  		




  
  56
  
    
in_dir_list_filename <- "/data/project/layers/commons/NEX_data/regions_input_files/stage6_in_dir_list_02052016.txt"


  		




  50
  57
  
    
#in_dir <- "" #PARAM 0


  		




  51
  58
  
    
#y_var_name <- "dailyTmax" #PARAM1


  		




  52
  59
  
    
#interpolation_method <- c("gam_CAI") #PARAM2


  		




  ......




  86
  93
  
    

  		




  87
  94
  
    
#run_assessment_plotting_prediction_fun(list_param_run_assessment_plottingin_dir) 


  		




  88
  95
  
    

  		




  89
  
  
    
run_assessment_plotting_prediction_fun <-function(list_param_run_assessment_plotting){


  		




  
  96
  
    
run_assessment_combined_region_plotting_prediction_fun <-function(list_param_run_assessment_plotting){


  		




  90
  97
  
    
  


  		




  91
  98
  
    
  ####


  		




  92
  99
  
    
  #1) in_dir: input directory containing data tables and shapefiles for plotting #PARAM 0


  		




  ......




  144
  151
  
    
  source(file.path(script_path,function_assessment_part2_functions)) #source all functions used in this script 


  		




  145
  152
  
    

  		




  146
  153
  
    
  ####### PARSE INPUT ARGUMENTS/PARAMETERS #####


  		




  147
  
  
    
  in_dir <- "/data/project/layers/commons/NEX_data"


  		




  148
  
  
    
  in_dir_list <-  c("output_run_global_analyses_pred_2009_reg4","output_run_global_analyses_pred_2010_reg4",


  		




  149
  
  
    
                    "output_run_global_analyses_pred_2011_reg4","output_run_global_analyses_pred_2012_reg4",


  		




  150
  
  
    
                    "output_run_global_analyses_pred_2013_reg4","output_run_global_analyses_pred_2014_reg4")


  		




  151
  
  
    
  


  		




  
  154
  
    
  list_param_run_assessment_plotting$in_dir_list_filename #PARAM 0


  		




  152
  155
  
    
  in_dir <- list_param_run_assessment_plotting$in_dir #PARAM 1


  		




  153
  156
  
    
  y_var_name <- list_param_run_assessment_plotting$y_var_name #PARAM2


  		




  154
  157
  
    
  interpolation_method <- list_param_run_assessment_plotting$interpolation_method #c("gam_CAI") #PARAM3


  		




  ......




  943
  946
  
    
df_assessment_files <- read.table(df_assessment_files_name,stringsAsFactors=F,sep=",")


  		




  944
  947
  
    
#threshold_missing_day <- c(367,365,300,200) #PARM18


  		




  945
  948
  
    

  		




  946
  
  
    
list_param_run_assessment_plotting <-list(


  		




  
  949
  
    
list_param_run_assessment_plotting <-list(  list_param_run_assessment_plotting,


  		




  947
  950
  
    
    in_dir,y_var_name, interpolation_method, out_suffix,


  		




  948
  951
  
    
    out_dir, create_out_dir_param, mosaic_plot, proj_str, file_format, NA_flag_val,


  		




  949
  952
  
    
    multiple_region, countries_shp, plot_region, num_cores,


  		




  950
  953
  
    
    region_name, df_assessment_files_name, threshold_missing_day,year_predicted


  		




  951
  954
  
    
  )


  		




  952
  955
  
    

  		




  953
  
  
    
names(list_param_run_assessment_plotting) <- c(


  		




  
  956
  
    
names(list_param_run_assessment_plotting) <- c("list_param_run_assessment_plotting$in_dir_list_filename #PARAM 0"


  		




  954
  957
  
    
    "in_dir","y_var_name","interpolation_method","out_suffix",


  		




  955
  958
  
    
    "out_dir","create_out_dir_param","mosaic_plot","proj_str","file_format","NA_flag_val",


  		




  956
  959
  
    
    "multiple_region","countries_shp","plot_region","num_cores",


  		












Also available in:  Unified diff