Project

General

Profile

« Previous | Next » 

Revision 4a25edab

Added by Benoit Parmentier over 8 years ago

global product assessment script part 1, time series profile at stations

View differences:

climate/research/oregon/interpolation/global_product_assessment_part1.R
4 4 
#Combining tables and figures for individual runs for years and tiles.
5 5 
#AUTHOR: Benoit Parmentier
6 6 
#CREATED ON: 05/15/2016
7 
#MODIFIED ON: 05/23/2016
7 
#MODIFIED ON: 05/24/2016
8 8 
#Version: 1
9 9 
#PROJECT: Environmental Layers project
10 10 
#COMMENTS: Initial commit, script based on part NASA biodiversity conferenc
......
80 80 
source(file.path(script_path,function_assessment_part3)) #source all functions used in this script
81 81 

  		




  82
  82
  
    
#Product assessment


  		




  83
  
  
    
function_product_assessment_part1_functions <- "global_product_assessment_part1_05242016.R"


  		




  
  83
  
    
function_product_assessment_part1_functions <- "global_product_assessment_part1_functions_05242016b.R"


  		




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


  		




  85
  85
  
    

  		




  86
  86
  
    
###############################


  		




  ......




  687
  687
  
    

  		




  688
  688
  
    
data_stations_tmp$date <- as.Date(strptime(data_stations_tmp$date,"%Y%m%d"))


  		




  689
  689
  
    
#data_stations_tmp$date <- as.Date(strptime(data_stations_tmp$date_str,"%Y%m%d"))


  		




  690
  
  
    
d_z4 <- 


  		




  
  690
  
    
#d_z4 <- 


  		




  691
  691
  
    
d_z_tmp4 <- zoo(data_stations_tmp$dailyTmax,data_stations_tmp$date)


  		




  692
  
  
    
plot(d_z_tmp4,cex.main=1.5,cex.lab=1.5,font.lab=2,


  		




  
  692
  
    
plot(d_z_tmp,cex.main=1.5,cex.lab=1.5,font.lab=2,


  		




  693
  693
  
    
     lty=3)


  		




  694
  
  
    
lines(d_z_tmp,ylab="tmax in deg C",xlab="Daily time steps",


  		




  
  694
  
    
lines(d_z_tmp4,ylab="tmax in deg C",xlab="Daily time steps",


  		




  695
  695
  
    
     main=title_str,cex=3,font=2,


  		




  696
  696
  
    
     col="red",


  		




  697
  697
  
    
     cex.main=1.5,cex.lab=1.5,font.lab=2,


  		




  ......




  706
  706
  
    
end_year <- year(end_date)


  		




  707
  707
  
    
d_z4 <- window(d_z_tmp4,start=start_date,end=end_date)


  		




  708
  708
  
    

  		




  709
  
  
    
###


  		




  
  709
  
    
data_stations_tmp$date[7190]


  		




  
  710
  
    

  		




  
  711
  
    
###TO DO:


  		




  
  712
  
    
#1) compute correlation r, RMSE, MAE etc by station and profile


  		




  
  713
  
    
#2) Compute temporal autocorrelation by station and profile


  		




  
  714
  
    
#3) 


  		




  
  715
  
    

  		




  
  716
  
    

  		




  710
  717
  
    
############################ END OF SCRIPT ##################################


  		












Also available in:  Unified diff