# Plotting results from LPJ-GUESS versions >= 4

ggplot2 (R) code to plot the simulations results in txt format from LPJ-GUESS model output, for multiple taxa, for the analysis in the article Ekberzade et al., Latitude or altitude will save us? A centennial case for forests in Asia Minor and its immediate surroundings 

## LAI_Plots (LAI.R)

Input files which are produced as large spreadsheets/matrices with Lon Lat Date LAI as parameters, are turned into dataframes, filtered for individual taxa and their LAI for each year is printed as a map and written into their respective folders: 
The input files are read as a table, partitioned and for each taxa, plotted as a map through a loop for each simulation year separately, and printed as .png to a specified folder. 
The code is suitable for area, region and global model output (may be cumbersome for site specific research, line plots may be preferred in that case), printing them out as a map. 
For further analysis, the code can be changed for fixed scale on the Y axis (necessary note included inside the code as comment.)

## Density dominance, per stand (Density_Dominance.R) 

Input files which are produced as large spreadsheets/matrices with Lon Lat Date density as parameters, are turned into dataframes, filtered for dominant taxa per time period, and printed as a map and written into a specified folder: 
The input files are read as a table, partitioned and for each taxa, dominance per grid is calculated and plotted as a map for the mean of a time period, and additionally through a loop for each simulation year separately, and printed as .png to a specified folder.

The code is suitable for area, region and global model output (may be cumbersome for site specific research, line plots may be preferred in that case), printing them out as a map. 

## Building large dataframes for bivariate spatio-temporal analysis using GIS software (Raster4GIS_and_SAD.R)

Merged datasets of median (from individual GCM simulations) biomass and density per stand (calculated as density/ha) as well as density dominance/stand, and respective biomass for certain groups of taxa (Shrubs, Temperate, Boreal, etc.) calculated as means of timeperiods. Additionally, Sum of Absolute Deviations (SAD) calculated for ensemble median from certain period means.