Revision 79b398dd
Added by Benoit Parmentier about 9 years ago
| climate/research/oregon/interpolation/global_run_scalingup_mosaicing.R | ||
|---|---|---|
| 5 | 5 |
#Analyses, figures, tables and data are also produced in the script. |
| 6 | 6 |
#AUTHOR: Benoit Parmentier |
| 7 | 7 |
#CREATED ON: 04/14/2015 |
| 8 |
#MODIFIED ON: 10/28/2015
|
|
| 8 |
#MODIFIED ON: 11/09/2015
|
|
| 9 | 9 |
#Version: 5 |
| 10 | 10 |
#PROJECT: Environmental Layers project |
| 11 | 11 |
#COMMENTS: analyses run for reg4 1992 for test of mosaicing using 1500x4500km and other tiles |
| ... | ... | |
| 23 | 23 |
#MODULEPATH=$MODULEPATH:/nex/modules/files |
| 24 | 24 |
#module load pythonkits/gdal_1.10.0_python_2.7.3_nex |
| 25 | 25 |
# |
| 26 |
#setfacl -Rm user:aguzman4:rwx /nobackupp8/bparmen1/output_run10_1500x4500_global_analyses_pred_1992_10052015 |
|
| 26 | 27 |
# |
| 27 | 28 |
################################################################################################# |
| 28 | 29 |
|
| ... | ... | |
| 57 | 58 |
|
| 58 | 59 |
#### FUNCTION USED IN SCRIPT |
| 59 | 60 |
|
| 60 |
function_mosaicing <-"global_run_scalingup_mosaicing_function_10282015.R"
|
|
| 61 |
function_mosaicing <-"global_run_scalingup_mosaicing_function_11092015.R"
|
|
| 61 | 62 |
|
| 62 | 63 |
#in_dir_script <-"/home/parmentier/Data/IPLANT_project/env_layers_scripts" #NCEAS UCSB |
| 63 | 64 |
in_dir_script <- "/nobackupp8/bparmen1/env_layers_scripts" #NASA NEX |
| ... | ... | |
| 199 | 200 |
|
| 200 | 201 |
mosaic_method <- "use_edge_weights" #this is distance from edge |
| 201 | 202 |
out_suffix_tmp <- paste(interpolation_method,y_var_name,day_to_mosaic[i],out_suffix,sep="_") |
| 202 |
debug(mosaicFiles) |
|
| 203 |
#debug(mosaicFiles)
|
|
| 203 | 204 |
#can also loop through methods!!! |
| 204 | 205 |
#python_bin <- "/usr/bin/" #python gdal bin, on Atlas NCEAS |
| 205 | 206 |
#python_bin <- "/nobackupp6/aguzman4/climateLayers/sharedModules/bin" #on NEX |
| ... | ... | |
| 270 | 271 |
lf_mean_mosaic <- lf_plot |
| 271 | 272 |
|
| 272 | 273 |
|
| 273 |
list_param_plot_mosaic <- list(lf_raster_fname=unlist(lf_mean_mosaic[1:3]),
|
|
| 274 |
screenRast=FALSE,
|
|
| 274 |
list_param_plot_mosaic <- list(lf_raster_fname=unlist(lf_mean_mosaic[1:2]),
|
|
| 275 |
screenRast=TRUE,
|
|
| 275 | 276 |
l_dates=day_to_mosaic, |
| 276 | 277 |
out_dir_str=out_dir, |
| 277 | 278 |
out_prefix_str <- "dailyTmax_", |
| ... | ... | |
| 279 | 280 |
#plot_screen_raster_val(3,list_param_plot_mosaic) |
| 280 | 281 |
#debug(plot_screen_raster_val) |
| 281 | 282 |
|
| 282 |
num_cores <- 3
|
|
| 283 |
l_png_files <- mclapply(1:length(unlist(lf_mean_mosaic)[1:3]),FUN=plot_screen_raster_val,
|
|
| 283 |
num_cores <- 1
|
|
| 284 |
l_png_files <- mclapply(1:length(unlist(lf_mean_mosaic)[1:2]),FUN=plot_screen_raster_val,
|
|
| 284 | 285 |
list_param= list_param_plot_mosaic, |
| 285 | 286 |
mc.preschedule=FALSE,mc.cores = num_cores) |
| 286 | 287 |
|
Also available in: Unified diff
testing plotting of mosaics