Project

General

Profile

« Previous | Next » 

Revision c1bdf4b3

Added by Benoit Parmentier over 11 years ago

run prediction follwing modifications for year 2010 Venezuela with monthly 1981-2010 averages

View differences:

climate/research/oregon/interpolation/master_script_temp.R
54 54 

  		




  55
  55
  
    
#/home/parmentier/Data/IPLANT_project/Venezuela_interpolation/Venezuela_01142013/climatology_03192013.py


  		




  56
  56
  
    
source(file.path(script_path,"covariates_production_temperatures_03212013.R"))


  		




  57
  
  
    
source(file.path(script_path,"Database_stations_covariates_processing_function_03272013.R"))


  		




  
  57
  
    
source(file.path(script_path,"Database_stations_covariates_processing_function_03242013.R"))


  		




  58
  58
  
    
source(file.path(script_path,"GAM_fusion_analysis_raster_prediction_multisampling_03272013.R"))


  		




  59
  59
  
    
source(file.path(script_path,"results_interpolation_date_output_analyses_03272013.R"))


  		




  60
  60
  
    
#source(file.path(script_path,"results_covariates_database_stations_output_analyses_03272013.R"))


  		




  ......




  75
  75
  
    
############ STAGE 2: Covariate production ################


  		




  76
  76
  
    

  		




  77
  77
  
    
##Paths to inputs and output


  		




  78
  
  
    
var<-"TMIN"


  		




  
  78
  
    
#var<-"TMIN"


  		




  
  79
  
    
var<-"TMAX"


  		




  79
  80
  
    
in_path <- "/home/parmentier/Data/IPLANT_project/Venezuela_interpolation/Venezuela_01142013/input_data/"


  		




  80
  81
  
    
out_path<- "/home/parmentier/Data/IPLANT_project/Venezuela_interpolation/Venezuela_01142013/output_data/"


  		




  81
  82
  
    

  		




  ......




  116
  117
  
    

  		




  117
  118
  
    
#Setting up input argurments for script function...


  		




  118
  119
  
    
#set up earlier


  		




  119
  
  
    
#var <- "TMIN"           # name of the variables to keep: TMIN, TMAX or PRCP --already set up earlier


  		




  
  120
  
    
var <- "TMAX"           # name of the variables to keep: TMIN, TMAX or PRCP --already set up earlier


  		




  120
  121
  
    

  		




  121
  
  
    
#infile_covariates<-"covariates__venezuela_region__VE_01292013.tif" #this is an output from covariate script and used in stage 3 and stage 4


  		




  122
  
  
    
infile_covariates<-"covariates__venezuela_region_TMIN__VE_03192013.tif" #covariates stack for TMIN


  		




  
  122
  
    
infile_covariates<-"covariates__venezuela_region__VE_01292013.tif" #this is an output from covariate script and used in stage 3 and stage 4


  		




  
  123
  
    
#infile_covariates<-"covariates__venezuela_region_TMIN__VE_03192013.tif" #covariates stack for TMIN


  		




  123
  124
  
    

  		




  124
  125
  
    
CRS_locs_WGS84<-CRS("+proj=longlat +ellps=WGS84 +datum=WGS84 +towgs84=0,0,0") #Station coords WGS84: same as earlier


  		




  125
  126
  
    
infile1<- "outline_venezuela_region__VE_01292013.shp"      #This is the shape file of outline of the study area                                                      #It is an input/output of the covariate script


  		




  ......




  127
  128
  
    

  		




  128
  129
  
    
#specific to this stage


  		




  129
  130
  
    
db.name <- "ghcn"       # name of the Postgres database


  		




  130
  
  
    
range_years<-c("2000","2001") #right bound not included in the range!!


  		




  
  131
  
    
range_years<-c("2010","2011") #right bound not included in the range!!


  		




  131
  132
  
    
range_years_clim<-c("1981","2011") #right bound not included in the range!!


  		




  132
  133
  
    
infile2<-"/home/layers/data/climate/ghcn/v2.92-upd-2012052822/ghcnd-stations.txt"                              #This is the textfile of station locations from GHCND


  		




  133
  134
  
    
in_path <- "/home/parmentier/Data/IPLANT_project/Venezuela_interpolation/Venezuela_01142013/input_data/"


  		




  134
  
  
    
out_prefix<-"_365d_GAM_fus5_all_lstd_03292013"                #User defined output prefix


  		




  
  135
  
    
out_prefix<-"_365d_GAM_fus5_all_lstd_03282013c"                #User defined output prefix


  		




  135
  136
  
    
qc_flags_stations<-c("0","S")    #flags allowed for screening after the query from the GHCND??


  		




  136
  137
  
    

  		




  137
  138
  
    
#list of 12 parameters for input in the function...


  		












Also available in:  Unified diff