# Data sources: 
   - Francesca Jaroszynska, Climate and biotic interactions-drivers of plant community structure and ecosystemfunctioning in alpine grasslands, PhD Thesis, 2019, 
   - Klanderud K, Vandvik V, Goldberg D (2015) The Importance of Biotic vs. Abiotic Drivers of Local Plant Community Composition Along Regional Bioclimatic Gradients. PLoS ONE 10(6): e0130205. doi:10.1371/journal.pone.0130205

# Annual_average_Precip_mm:  use "RAIN_FROM_ATM" from monthly model output, need to do multi-year annual average (e.g., last 30 years simulation), convert the unit (mm/s) to (mm/year) by multiplying 3600*24*365

# Mean_temperature_6-9months_°C: use "TSA" from monthly model output, select month 6-9 to do multi-year average (e.g., last 30 years simulation), convert the unit (K) to (oC) by substract 273.15.

# Graminoid%, Forbs%, Grass(G+F)% : use "BA_SCPF" from monthly model output. 
                                    Note that this variable is not a default output variable, need to add "hist_fincl1 = 'BA_SCPF'" in "user_nl_clm"
                                    First need to do multi-year annual average (e.g., last 30 years simulation), then also do "PFT" averaing (using "fates_pftmap_levscpf" in the output)
                                    convert unit (m2/ha) to (%) by dividing 10000, grass pfts should have index 10, 11, 12, so compare values from these PFTs with observation value

# vegetation_height_cm : use "HTOP" from monthly model output. convert (m) to (cm) by multiply 100