Project

General

Profile

« Previous | Next » 

Revision 09d89ec0

Added by Benoit Parmentier over 9 years ago

generating mosaicing for reg4, test with automated code

View differences:

climate/research/oregon/interpolation/global_run_scalingup_mosaicing.R
46 46 

  		




  47
  47
  
    
#### FUNCTION USED IN SCRIPT


  		




  48
  48
  
    

  		




  49
  
  
    
#function_analyses_paper1 <-"contribution_of_covariates_paper_interpolation_functions_07182014.R" #first interp paper


  		




  50
  
  
    
#function_analyses_paper2 <-"multi_timescales_paper_interpolation_functions_08132014.R"


  		




  51
  
  
    

  		




  52
  49
  
    
function_mosaicing <-"global_run_scalingup_mosaicing_function_06212015.R"


  		




  53
  50
  
    

  		




  54
  51
  
    
in_dir_script <-"/home/parmentier/Data/IPLANT_project/env_layers_scripts"


  		




  ......




  59
  56
  
    

  		




  60
  57
  
    
#Data is on ATLAS: reg4 (South America)


  		




  61
  58
  
    

  		




  62
  
  
    
in_dir <- "/data/project/layers/commons/NEX_data/mosaicing_data_test" #reg1 is North America and reg5 is Africa


  		




  63
  
  
    
#in_dir <- "/data/project/layers/commons/NEX_data/mosaicing_data_test/reg1"


  		




  
  59
  
    
in_dir <- "/data/project/layers/commons/NEX_data/mosaicing_data_test" 


  		




  
  60
  
    
#in_dir <- "/data/project/layers/commons/NEX_data/mosaicing_data_test/reg1" #North America


  		




  64
  61
  
    
#in_dir <- "/data/project/layers/commons/NEX_data/mosaicing_data_test/reg2" #Europe


  		




  
  62
  
    
#in_dir <- "/data/project/layers/commons/NEX_data/mosaicing_data_test/reg4" #South America


  		




  65
  63
  
    
#in_dir <- "/data/project/layers/commons/NEX_data/mosaicing_data_test/reg5" #Africa


  		




  66
  64
  
    

  		




  67
  65
  
    
y_var_name <- "dailyTmax" #PARAM1


  		




  68
  66
  
    
interpolation_method <- c("gam_CAI") #PARAM2


  		




  69
  
  
    
region_name <- "reg2" #PARAM 13 #reg4 South America, Africa reg5,Europe reg2, North America reg1, Asia reg3


  		




  
  67
  
    
region_name <- "reg4" #PARAM 13 #reg4 South America, Africa reg5,Europe reg2, North America reg1, Asia reg3


  		




  70
  68
  
    
mosaicing_method <- c("unweighted","use_edge_weights")


  		




  71
  69
  
    
out_suffix <- paste(region_name,"_","mosaic_run10_1500x4500_global_analyses_06212015",sep="") 


  		




  72
  70
  
    
#PARAM3


  		




  ......




  162
  160
  
    
#undebug(plot_mosaic)


  		




  163
  161
  
    
#plot_mosaic(1,list_param=list_param_plot_mosaic)


  		




  164
  162
  
    
num_cores <- 4


  		




  165
  
  
    
l_png_files <- mclapply(1:length(lf_mean_mosaic),FUN=plot_mosaic,list_param= list_param_plot_mosaic,


  		




  
  163
  
    
l_png_files <- mclapply(1:length(unlist(lf_mean_mosaic)),FUN=plot_mosaic,


  		




  
  164
  
    
                        list_param= list_param_plot_mosaic,


  		




  166
  165
  
    
                        mc.preschedule=FALSE,mc.cores = num_cores)


  		




  167
  166
  
    

  		




  168
  167
  
    
####################


  		












Also available in:  Unified diff