Set the soil type of a FarmPlot to the most frequent soil type of its pixels

src/crop/cropmodels.jl

@@ -3,6 +3,8 @@
 ### Functions for switchable crop models.
 ###
 
+using StatsBase: mode
+
 """
     initfields!(cropmodel, cropfile, growthfile)
 
@@ -61,6 +63,12 @@ function initfields!(model::SimulationModel)
             end
         end
     end
+
+    # Adjust soil type to most common soil type of the pixels
+    # belonging to a farmplot
+    for fp in model.farmplots
+        fp.soiltype = mode(map(p -> model.landscape[p[1], p[2]].soiltype, fp.pixels))
+    end
     @info "Initialised $(length(model.farmplots)) farm plots."
 end