get.site.data.R

get.site.data <- function(site.number) {
  # The function get.site.data takes an integer input corresponding to the
  # sites listed in the file site_metadata.csv and creates/updates the GCC and
  # MODIS data files for that site. 
  
  gcc_filename <- sprintf("gcc_data_site%i.csv",site.number)
  ndvi_filename <- sprintf("ndvi_data_site%i.csv",site.number)
  
  # Check to see if some data has already been downloaded (i.e. if the files
  # already exist):
  some.data.downloaded <- file.exists(gcc_filename) &&  file.exists(ndvi_filename) # some.data.downloaded is TRUE/FALSE
  
  # If there isn't at least SOME GCC and MODIS data, then download all
  # of the gcc and MODIS data from 2000-present, and process it. Save the 
  # data in the file SOMETHING2.CSV:
  if(!some.data.downloaded){
    
    # Download phenocam data and save (creates file named "phenocam_data_siteX.csv", 
    # where X is the site number):
    print(sprintf("Getting phenocam data for site %i...",site.number))
    source("download.phenocam.data.R")
    download.phenocam.data(site.number)
    
    # Process phenocam data into gcc data:
    print(sprintf("Processing phenocam data for site %i...",site.number))
    source("create.gcc.data.R")
    create.gcc.data(site.number) # creates file gcc_data_siteX.csv, where X = site.number
    
    # Before we download the MODIS data, we'll remove the file that we'll 
    # get a filename from:
    unlink("Subset Download*.csv")
    
    # Download ALL MODIS data:
    print(sprintf("Downloading MODIS data for site %i...",site.number))
    source("download.all.modis.data.R")
    download.all.modis.data(site.number)
    
    # Process MODIS data into NDVI data:
    print(sprintf("Processing MODIS data for site %i...",site.number))
    source("create.ndvi.data.R")
    create.ndvi.data(site.number)    
  }
  else{ #ie if some.data.downloaded is TRUE, just need to update:
    
    # Probably simplest to just re-download all of the gcc data:
    print(sprintf("Getting phenocam data for site %i...",site.number))
    source("download.phenocam.data.R")
    download.phenocam.data(site.number)
    
    print(sprintf("Processing phenocam data for site %i...",site.number))
    source("create.gcc.data.R")
    create.gcc.data(site.number) # creates file gcc_data_siteX.csv, where X = site.number
    
    
    # Just need to download the last year of MODIS data (SUPER SLOW!!), and then add
    # it with the existing data:
    unlink("Subset Download*.csv")
    print(sprintf("Downloading MODIS data for site %i...",site.number))
    source("download.new.modis.data.R")
    download.new.modis.data(site.number)
    
    print(sprintf("Processing MODIS data for site %i...",site.number))
    source("update.ndvi.data.R")
    update.ndvi.data(site.number)
  }
  
}