Project

General

Profile

« Previous | Next » 

Revision d622431a

Added by Benoit Parmentier over 9 years ago

assessment part2 1982 reg5 test

View differences:

climate/research/oregon/interpolation/global_run_scalingup_assessment_part2.R
5 5
#Analyses, figures, tables and data are also produced in the script.
6 6
#AUTHOR: Benoit Parmentier 
7 7
#CREATED ON: 03/23/2014  
8
#MODIFIED ON: 09/03/2015            
8
#MODIFIED ON: 09/21/2015            
9 9
#Version: 4
10 10
#PROJECT: Environmental Layers project     
11 11
#COMMENTS: analyses for run 10 global analyses,all regions 1500x4500km with additional tiles to increase overlap 
......
381 381
interpolation_method <- c("gam_CAI") #PARAM2
382 382
#out_suffix<-"run10_global_analyses_01282015"
383 383
#out_suffix <- "output_run10_1000x3000_global_analyses_02102015"
384
out_suffix <- "run10_1500x4500_global_analyses_pred_1992_09012015" #PARAM3
385
out_dir <- "/data/project/layers/commons/NEX_data/output_run10_1500x4500_global_analyses_pred_1992_09012015" #PARAM4
384
out_suffix <- "run10_1500x4500_global_analyses_pred_1982_09152015" #PARAM3
385
out_dir <- "/data/project/layers/commons/NEX_data/output_run10_1500x4500_global_analyses_pred_1982_09152015" #PARAM4
386 386
create_out_dir_param <- FALSE #PARAM 5
387 387

  
388 388
mosaic_plot <- FALSE #PARAM6
389 389

  
390 390
#if daily mosaics NULL then mosaicas all days of the year
391 391

  
392
day_to_mosaic <- c("19920101","19920102","19920103","19920104","19920105",
393
                   "19920106","19920107","19920108","19920109","19920110",
394
                   "19920111")
392
day_to_mosaic <- c("19820101","19820102","19820103","19820104","19820105",
393
                   "19820106","19820107","19820108","19820109","19820110",
394
                   "1982011")
395 395

  
396 396
CRS_WGS84 <- CRS("+proj=longlat +ellps=WGS84 +datum=WGS84 +towgs84=0,0,0") #Station coords WGS84 #CONSTANT1
397 397
CRS_locs_WGS84<-CRS("+proj=longlat +ellps=WGS84 +datum=WGS84 +towgs84=0,0,0") #Station coords WGS84
......
574 574
#centroids_pts <- tmp_pts 
575 575
  
576 576
#plot info: with labels
577
res_pix <- 1200
577
res_pix <-1200
578 578
col_mfrow <- 1 
579 579
row_mfrow <- 1
580 580

  
......
593 593
    text(coordinates(pt)[1],coordinates(pt)[2],labels=i,cex=1.3,font=2,col=c("red"))
594 594
  }
595 595
}
596
title(paste("Tiles ", tile_size,region_name,sep=""))
596
#title(paste("Tiles ", tile_size,region_name,sep=""))
597 597

  
598 598
dev.off()
599 599
      

Also available in: Unified diff