Skip to content

Check calculations end reporter

Nach einem Modelrun (= 1 Jahr) möchte ich u.a. folgendes reporten

  1. yearly-biomass
  2. average-pesticide-load bzw. print-average-pesticide-load in relation to biomass

Für 1. habe ich momentan folgenden Vorgang:

to print-yearly-biomass ;Print total emerged biomass of aquatic emerging insects at the end of the model year in g
  let biomass-caddisfly (precision (total-caddisfly-emerged * 0.00794) 2) ;based on Sabo et al. (2002), Roodt et al. (2023) and Roodt (personal communciation, 2025)
  let biomass-mayfly (precision (total-mayfly-emerged * 0.00184) 2)
  let biomass-midge (precision (total-midge-emerged * 0.00068) 2)

  if spawn-caddisflies? [print (word "Total biomass caddisfly (g): " biomass-caddisfly)]
  if spawn-mayflies? [print (word "Total biomass mayfly (g): " biomass-mayfly)]
  if spawn-midges? [print (word "Total biomass midge (g): " biomass-midge)]
end

Für 2. habe ich:
```netlogo
to print-average-pesticide-load ;Print average pesticide-load per year
  if spawn-caddisflies? [print (word "Avg caddisfly pesticide-load per year: " (precision (total-caddisfly-load / total-caddisfly-emerged) 3) " (ng pesticide in individual insect) and " (precision ((total-caddisfly-load / total-caddisfly-emerged) / 0.00794) 2) " (ng/g dry weight)")]
  if spawn-mayflies? [print (word "Avg mayfly load per year: " (precision (total-mayfly-load / total-mayfly-emerged) 3) " (ng pesticide in individual insect) and " (precision ((total-mayfly-load / total-mayfly-emerged) / 0.00184) 2) " (ng/g dry weight)")]
  if spawn-midges? [print (word "Avg midge load per year (ng pesticide in individual insect): " (precision (total-midge-load / total-midge-emerged) 3) " (ng pesticide in individual insect) and " (precision ((total-midge-load / total-midge-emerged) / 0.00068) 2) " (ng/g dry weight)")]
end

Meine hauptsächliche Frage hierbei:

a) ist es richtig, wie ich momentan total-caddisfly-emerged; total-mayflies-emerged und total-midges-emerged tracke über

set total-caddisfly-emerged total-caddisfly-emerged + daily-caddisfly-emerged
  set total-mayfly-emerged total-mayfly-emerged + daily-mayflies-emerged
  set total-midge-emerged total-midge-emerged + daily-midges-emerged

in meinem to emerge-insect spawne?

b) ist es richtig, wie ich total-caddisfly-load; total-mayfly-load und total-midge load tracke, über:

to accumulate-pesticide-load ;used to track accumulated pesticide-load
  set total-caddisfly-load total-caddisfly-load + sum [pesticide-load] of caddisflies with [just-spawned? = false]
  set total-mayfly-load total-mayfly-load + sum [pesticide-load] of mayflies with [just-spawned? = false]
  set total-midge-load total-midge-load + sum [pesticide-load] of midges with [just-spawned? = false]
end

Das habe ich als eigenständigen primitive und nicht in to pesticide-uptake oder to update-pesticide integriert. Ist das ein Problem für NetLogo hinsichtlich der Berechnungen oder ist das egal?

Mit u.a. total-caddisfly-load kummuliere ich ja die pesticide-load aller Insekten einer Gruppe in ng. Wenn ich das dann durch die kummulierte Anzahl der gespawnten Insekten teile (total-caddisfly-emerged), sollte ich ja logischerweise die durchschnittliche pesticide-load errechnen. Da mich ebenfalls die pesticide-load in Abhängigkeit zur Biomasse interessiert, möchte ich das durch teilen der durchschnittlichen pesticide-load durch die Biomasse erreichen, also (total-caddisfly-load / total-caddisfly-emerged) / 0.00794 - Einheit wäre dann ng/g.

(Jede caddisfly hat die gleiche Biomasse, das ist der Durchschnitt, den ich aus meinen Daten habe und der auch für die Pesticide-Uptake Berechnung genutzt wird - eine range wäre auch möglich, aber das vielleicht später noch)

Natürlich sind durch meine kleineren Medianwerte, Durchschnittswerte etc. eine Abweichung zu den Ergebnissen von Roodt et al. logisch, trotzdem erwarte ich aber eigentlich, dass die ng/g Werte für die Midges am höchsten sind bzw. deutlich höher als caddisflies und mayflies (höchster BWAF, kleinstes Dry-Weight, siehe Bild unten) - ich sehe aber hier aber fast identische Werte wie für caddisflies und mayflies und frage mich, ob es an einem Fehler beim Tracking der Werte liegt

image