BoundsError in ADBM rewiring model when using with nutrient intake model.

[evadelmas]

Seems that there is a problem in the ADBM's get_feeding_links function caused by the biomass object containing actual biomass AND nutrients concentration.

Expected behavior

• rewiring following extinction events

Actual behavior

• ERROR: BoundsError: attempt to access 4-element Array{Float64,1} at index [Base.LogicalIndex(Bool[false, true, false, false, false, false])] [...] at src/rewiring/ADBM.jl:94

Steps to reproduce

A = [0 1 1 0 ; 0 0 0 1 ; 0 0 0 1 ; 0 0 0 0]
p = model_parameters(A, Z = 10.0, rewire_method = :ADBM, Hmethod = :ratio, Nmethod = :biomass, h = 2.0, productivity = :nutrients)
biomass = [0.2, 0.0, 0.4, 0.7]
s = simulate(p, biomass)