Added units to insectbiomass()

src/Persefone.jl

@@ -100,12 +100,13 @@ export
     loadmodelobject
 
 ## Import and define units and dimensions
-import Unitful: cm, m, km, ha, Length, Area
+import Unitful: cm, m, km, ha, mg, g, kg, Length, Area, Mass
 const m² = m^2
 const km² = km^2
 import Base./ # enable division with different length/area unit types
 /(x::S,y::T) where {S<:Length, T<:Length} = (upreferred(x)/m) / (upreferred(y)/m)
 /(x::S,y::T) where {S<:Area, T<:Area} = (upreferred(x)/m²) / (upreferred(y)/m²)
+/(x::S,y::T) where {S<:Mass, T<:Mass} = (upreferred(x)/g) / (upreferred(y)/g)
 
 """
     SimulationModel

src/nature/insects.jl

@@ -28,8 +28,8 @@ function insectbiomass(pixel::Pixel, model::SimulationModel)::Float64
 
     ## if no factors are configured, insect abundance defaults to 300 mg/m²,
     ## a value in the upper range of insect biomass density in agricultural landscapes
-    baseline = 300
-    seasonfactor = 0.0
+    baseline = 300mg/m²
+    seasonfactor = 0.0mg/m²
     habitatfactor = 1.0
     weatherfactor = 1.0
     pesticidefactor = 1.0
@@ -73,5 +73,5 @@ function insectbiomass(pixel::Pixel, model::SimulationModel)::Float64
 
     ## calculate biomass using a parabolic equation in the vertex form
     biomass = seasonfactor+baseline*habitatfactor*pesticidefactor*weatherfactor
-    biomass > 0 ? biomass/1000 : 0.0 # convert mg to g
+    biomass > 0mg/m² ? biomass |> g/m² : 0.0g/m² # convert mg to g
 end