Fixed the suitability bug. Now current and suitable forage grass and forage crop are having the same yields.

Did some validation. Still there is something weird with the forage crops. They should be identical to the current yield per country but they are not. Often in the suitability yield dataset they are less. Don't know why yet...

Finalized the validation of the EU28 dataset with bar charts and summary tables.

Did some debugging, removed very small numbers in area and production and did some validation (visual, summary)

Created an R script for manipulating the GHG EFs and to adjust them to the new zones

Created a new folder with R script to calculate NP contents in R for forage crops from Nutrient tool database https://plantsorig.sc.egov.usda.gov/npk/main

external function script

Added the suitability crops to the excel > use the distinct crop list -> the other is not updated with suitability crops.

Creating a dataset with rangeland, arable land, arable pasture, natural pasture per zone (adm0, aez, soilagri). draft, not finished.

Final dataset for the globe and for EU28. Also cifos crop names are added. Issues are the not available disaggregation and the impication on nutrient intake (Rest, tnut, Veggie). Which nutrient content should be taken as a proxy?

all crops done. treenuts, veggies, EU subsetting. Function R script created. Suitability file created and transformed, linked to zones AEZ ADM0 Soilagri.

New zonation based on ADM0 - AEZ - DOMSOIL_IPCC(mineral and Organic) - Validation with EUROSTAT and FAO. Data transformation. Different types of agricultural land extractions

Function, defining the different surogate crop species for aggregated groups.

Peat done (except of merging the yield and area of grass >> still some issues). Single files are fine for crop and grassland area now with peat and nonpeat.
SPAM extraction from GEOTIFF file >> New zonal statistics function > much more efficient now. Replaced left_join with bind function. Now much less memory demanding.
Extracted factors from columnames.

Extraction of SPAM geotiff files based on ULrikes recommendation. Then I can use GAUL as a consistent dataset.

Deleted all old carbage from the script.

Started with the creation of the EU crop area validation data >> still needs finishing (see issues on git)

> No NA's > Only the **ratio of harvested area** for disaggregation. Like this, there is more consistency with the presence of nuts/veggies in adm1 zones. Otherwise, it happened that there is area but no production vice versa. Then you divide by 0 = inf.
> No inf values: The inf values came from the division of P/H=Y calculation. The P was in all cases very low so I decided to simply recode them to 0. Then Yields would be also 0.
> rounding fixed > 1 decimal like in the original dataset so we do not lose the level of detail.
> Yield conversion (production in mt) checked and correct.
> Aggregation to adm1 checked and correct.