Change grid dimensions in example provided
Remove argument 'output_file' from function
Add several require statements
Replace hardcoded domain dimensions with variable
Replace code for determining extents with function calls
Modify indentation in ggplot code
Write formatted heights text to disk as 'heights.txt'
Create data frame for land use classifications
Create a RasterLayer object with lat/lon grid cell coords
Create a SpatialPixels object from the RasterLayer object
Extract lat/lon coords from SpatialPixels 'coords' slot
Create vectors for starting and ending dates
Column bind date vectors with 'modis_coordinates' df
Get all bands from the MODIS MCD12Q1 product
Extract MODIS 'MCD12Q1' satellite data for all coords
Generate a file list for all of the extract data
Generate a vector of IGBP Type 1 land cover classes
Get the corresponding CALMET categories through a join
Create a data frame for the LU categories, in row-major order
Generate a vector of space-delimited strings for LU codes