[Frontpage] [Contents] [Previous] [Next]

Appendix 2: Oilseed rape model - equations and parameters

BiomBRACARR(t) = BiomBRACARR(t - dt) + (GrowBRACARR + SowBRACARR - HarvBRACARR - ControlBRACARR) * dt

INIT BiomBRACARR = 0

GrowBRACARR = if (MaxBRACA>BiomTotal) THEN (Hybridization*Winter*RGRBRACA*BiomBRACARR*(MaxBRACA-BiomTotal)/MaxBRACA)/dt else 0

SowBRACARR = Delay(BeginWeed,1)*SeedBRACARR*GermBRACA*TKVBRACA*10000/(1000*1000000)/dt

HarvBRACARR = if EndWeed=1 then (BiomBRACARR+0.1) else 0

ControlBRACARR = BiomBRACARR*EFBRARe

BiomBRACARS(t) = BiomBRACARS(t - dt) + (GrowBRACARS + SowBRACARS - HarvBRACARS - ControlBRACARS) * dt

INIT BiomBRACARS = 0

GrowBRACARS = if (MaxBRACA>BiomTotal) THEN (Hybridization*Winter*RGRBRACA*BiomBRACARS*(MaxBRACA-BiomTotal)/MaxBRACA)/dt else 0

SowBRACARS = Delay(BeginWeed,1)*SeedBRACARS*GermBRACA*TKVBRACA*10000/(1000*1000000)/dt

HarvBRACARS = if EndWeed=1 then (BiomBRACARS+0.1) else 0

ControlBRACARS = BiomBRACARS*EFBRARe

BiomBRACASS(t) = BiomBRACASS(t - dt) + (GrowBRACASS + SowBRACASS - HarvBRACASS - ControlBRACASS) * dt

INIT BiomBRACASS = 0

GrowBRACASS = if (MaxBRACA>BiomTotal) THEN (Hybridization*Winter*RGRBRACA*BiomBRACASS*(MaxBRACA-BiomTotal)/MaxBRACA)/dt else 0

SowBRACASS = Delay(BeginWeed,1)*SeedBRACASS*GermBRACA*TKVBRACA*10000/(1000*1000000)/dt

HarvBRACASS = if EndWeed=1 then (BiomBRACASS+0.1) else 0

ControlBRACASS = BiomBRACASS*EFBRASe

SeedBRACARR(t) = SeedBRACARR(t - dt) + (NewSeedBRACARR - DeadSeedBRACARR) * dt

INIT SeedBRACARR = 0

NewSeedBRACARR = ((HarvF1*0.25*HIBRACAF1)+((HarvBRACARR*HarvBRACARR+HarvBRACARS*HarvBRACARR+0.25*HarvBRACARS*HarvBRACARS)/BRACAtot)*HIBRACA)*NonPredBRACA*1000000*1000/(TKVBRACA*10000)/dt

DeadSeedBRACARR = (SeedBRACARR*(1-EXP(LOGN(1-DecayBRACA)/365))+BeginWeed*SeedBRACARR*GermBRACA)/dt

SeedBRACARS(t) = SeedBRACARS(t - dt) + (NewSeedBRACARS - DeadSeedBRACARS) * dt

INIT SeedBRACARS = 0

NewSeedBRACARS = ((HarvF1*0.5*HIBRACAF1)+((HarvBRACARR*HarvBRACARS+2*HarvBRACARR*HarvBRACASS+0.5*HarvBRACARS*HarvBRACARS+HarvBRACARS*HarvBRACASS)/BRACAtot)*HIBRACA)*NonPredBRACA*1000000*1000/(TKVBRACA*10000)/dt

DeadSeedBRACARS = (SeedBRACARS*(1-EXP(LOGN(1-DecayBRCA)/365))+BeginWeed*SeedBRACARS*GermBRACA)/dt

SeedBRACASS(t) = SeedBRACASS(t - dt) + (NewSeedBRACASS - DeadSeedBRACASS) * dt

INIT SeedBRACASS = 200

NewSeedBRACASS = ((HarvF1*0.25*HIBRACAF1)+((0.25*HarvBRACARS*HarvBRACARS+HarvBRACASS*HarvBRACARS+HarvBRACASS*HarvBRACASS)/BRACAtot))*HIBRACA*NonPredBRACA*1000000*1000/(TKVBRACA*10000)/dt

DeadSeedBRACASS = (SeedBRACASS*(1-EXP(LOGN(1-DecayBRACA)/365))+BeginWeed*SeedBRACASS*GermBRACA)/dt

HIBRACAF1 = 0.05

BiomWBarley(t) = BiomWBarley(t - dt) + (GrowWBarley + SowWBarley - HarvWBarley) * dt

INIT BiomWBarley = 0

GrowWBarley = if (MaxWBarley>BiomTotal) THEN (Winter*RGRWBarley*BiomWBarley*(MaxWBarley-BiomTotal)/MaxWBarley)/dt else 0

SowWBarley = if BeginWBarley=1 then SeedWBarley/dt else 0

HarvWBarley = if EndWBarley=1 then (BiomWBarley+0.1)/dt else 0

BiomWRape(t) = BiomWRape(t - dt) + (GrowWRape + SowWRape - HarvWRape) * dt

INIT BiomWRape = 0

GrowWRape = if (MaxWRape>BiomTotal) THEN (Winter*RGRWRape*BiomWRape*(MaxWRape-BiomTotal)/MaxWRape)/dt else 0

SowWRape = if BeginWRape=1 then SeedWRape/dt else 0

HarvWRape = if EndWRape=1 then (BiomWRape+0.1)/dt else 0

BiomWWheat(t) = BiomWWheat(t - dt) + (GrowWWheat + SowWWheat - HarvWWheat) * dt

INIT BiomWWheat = 0

GrowWWheat = if (MaxWWheat>BiomTotal) THEN (Winter*RGRWWheat*BiomWWheat*(MaxWWheat-BiomTotal)/MaxWWheat)/dt else 0

SowWWheat = IF(BeginWWheat)=1 THEN(SeedWWheat)/dt else 0

HarvWWheat = IF (EndWWheat)=1 THEN (BiomWWheat+0.1)/dt ELSE(0)

HIWBarley = 0.4

HIWRape = 0.25

HIWWheat = 0.40

MaxWBarley = 6.5*0.85/HIWBarley

MaxWRape = 3.8*0.92/HIWRape

MaxWWheat = 8.6*0.85/HIWWheat

RGRWBarley = 0.051

RGRWRape = 0.067

RGRWWheat = 0.042

SeedWBarley = 160*0.87/1000

SeedWRape = 5*0.92/1000

SeedWWheat = 180*0.87/1000

BiomF1(t) = BiomF1(t - dt) + (GrowF1 + SowF1 - HarvF1 - ControlF1) * dt

INIT BiomF1 = 0

GrowF1 = if (MaxWRape>BiomTotal) THEN (Hybridization*Winter*RGRWRape*BiomF1*(MaxWRape-BiomTotal)/MaxWRape)/dt else 0

SowF1 = (Delay(BeginWeed,1)*SeedF1*GermWRapeVol*TKWWRape/(1000*1000000))/dt

HarvF1 = if EndWeed=1 then (BiomF1+0.1)/dt ele 0

ControlF1 = BiomF1*EFBRARe

SeedF1(t) = SeedF1(t - dt) + (NewSeedF1 - DeadSeedF1) * dt

INIT SeedF1 = 0

NewSeedF1 = IF (Bn_Rr?=0) THEN (HarvBRACASS*1/(1.075+13.293*HarvBRACASS/(HarvBRACASS+HarvWRape))+(HarvWRape*LossWRape*1/(3.03+16.299*(HarvWRape/HarvBRACASS+HarvWRape))))*DormWRape*HIBRACA*1000000*1000/TKWWRape/dt

ELSE

0.5*(HarvBRACASS*1/(1.075+13.293*HarvBRACASS/(HarvBRACASS+HarvWRape))+(HarvWRape*LossWRape*1/(3.03+16.299*(HarvWRape/HarvBRACASS+HarvWRape))))*DormWRape*HIBRACA*1000000*1000/TKWWRape/dt

DeadSeedF1 = (SeedF1*(1-EXP(LOGN(1-DecayWRape)/365))+BeginWeed*SeedF1*GermWRapeVol)/dt

Hybridization = 1

BiomWRapeVolRe(t) = BiomWRapeVolRe(t - dt) + (GrowWRapeVolRe + SowWRapeVolRe - HarvWRapeVolRe - ControlWRapeVolRe) * dt

INIT BiomWRapeVolRe = 0

GrowWRapeVolRe = if (MaxWRape>BiomTotal) THEN (Winter*RGRWRape*BiomWRapeVolRe*(MaxWRape-BiomTotal)/MaxWRape)/dt else 0

SowWRapeVolRe = (Delay(BeginWeed,1)*(SeedBnReRe*GermWRapeVol*TKWWRape)/(1000*1000000))/dt

HarvWRapeVolRe = if EndWeed=1 then (BiomWRapeVolRe+0.1)/dt else 0

ControlWRapeVolRe = BiomWRapeVolRe*EFBRARe

BiomWRapeVolReSe(t) = BiomWRapeVolReSe(t - dt) + (GrowWRapeVolReSe + SowWRapeVolReSe - HarvWRapeVolReSe - ControlWRapeVolReSe) * dt

INIT BiomWRapeVolReSe = 0

GrowWRapeVolReSe = if (MaxWRape>BiomTotal) THEN (Winter*RGRWRape*BiomWRapeVolReSe*(MaxWRape-BiomTotal)/MaxWRape)/dt else 0

SowWRapeVolReSe = (Delay(BeginWeed,1)*(SeedBnReSe)*GermWRapeVol*TKWWRape/(1000*1000000))/dt

HarvWRapeVolReSe = if EndWeed=1 then (BiomWRapeVolReSe+0.1)/dt else 0

ControlWRapeVolReSe = BiomWRapeVolReSe*EFBRARe

BiomWRapeVolSe(t) = BiomWRapeVolSe(t - dt) + (GrowWRapeVolSe + SowWRapeVolSe - HarvWRapeVolSe - ControlWRapeVolSe) * dt

INIT BiomWRapeVolSe = 0

GrowWRapeVolSe = if (MaxWRape>BiomTotal) THEN (Winter*RGRWRape*BiomWRapeVolSe*(MaxWRape-BiomTotal)/MaxWRape)/dt else 0

SowWRapeVolSe = (Delay(BeginWeed,1)*SeedBnSeSe*GermWRapeVol*TKWWRape/(1000*1000000))/dt

HarvWRapeVolSe = if EndWeed=1 then (BiomWRapeVolSe+0.1)/dt else 0

ControlWRapeVolSe = BiomWRapeVolSe*EFBRASe

SeedBnReRe(t = SeedBnReRe(t - dt) + (NewSeedBnReRe - DeadSeedBnReRe) * dt

INIT SeedBnReRe = 0

NewSeedBnReRe = IF (HarvWRape>0.01) THEN ((HarvWRapeVolRe+(HarvWRape+HarvWRapeVolReSe)*0.25)*Bn_Rr?*LossWRape*DormWRape*HIWRape*1000000*1000/TKWWRape)/dt

ELSE

((HarvWRapeVolRe+HarvWRapeVolReSe*0.25)*Bn_Rr?*DormWRape*HIWRape*1000000*1000/TKWWRape)/dt

DeadSeedBnReRe = (SeedBnReRe*(1-EXP(LOGN(1-DecayWRape)/365))+BeginWeed*SeedBnReRe*GermWRapeVol)/dt

SeedBnReSe(t) = SeedBnReSe(t - dt) + (NewSeedBnReSe - DeadSeedBn_ReSe) * dt

INIT SeedBnReSe = 0

NewSeedBnReSe = IF (HarvWRape>0.01) THEN ((HarvWRape+HarvWRapeVolReSe)*0.5*Bn_Rr?*LossWRape*DormWRape*HIWRape*1000000*1000/TKWWRape)/dt

ELSE

((HarvWRapeVolReSe*0.5*Bn_Rr?)*DormWRape*HIWRape*1000000*1000/TKWWRape)/dt

DeadSeedBn_ReSe = (SeedBnReSe*(1-EXP(LOGN(1-DecayWRape)/365))+BeginWeed*SeedBnReSe*GermWRapeVol)/dt

SeedBnSeSe(t) = SeedBnSeSe(t - dt) + (NewSeedBnSeSe - DeadSeedBnSeSe) * dt

INIT SeedBnSeSe = 0

NewSeedBnSeSe = IF (HarvWRape>0.01) THEN ((HarvWRapeVolSe+(HarvWRape+HarvWRapeVolReSe)*0.25)*Bn_Rr?*LossWRape*DormWRape*HIWRape*1000000*1000/TKWWRape)/dt

ELSE

((HarvWRapeVolSe+HarvWRapeVolReSe*0.25)*Bn_Rr?*DormWRape*HIWRape*1000000*1000/TKWWRape)/dt

DeadSeedBnSeSe = (SeedBnSeSe*(1-EXP(LOGN(1-DecayWRape)/365))+BeginWeed*SeedBnSeSe*GermWRapeVol)/dt

Bn_Rr? = 0

BiomCereals = BiomWBarley+BiomWWheat

BiomCrop = BiomWBarley+BiomWRape+BiomWWheat

BiomTotal = BiomCrop+BiomVol+BiomWeeds

BiomVol = BiomWBarleyVol+BiomWWheatVol+DicotVol

BiomWeeds = DicotWeed+BiomELYRE

BRACAhybtot = BiomF1+BRACAtot+BiomBRACA

BRACAtot = BiomBRACARR+BiomBRACARS+BiomBRACASS

DicotVol = BiomWRapeVolRe+BiomWRapeVolReSe+BiomWRapeVolSe+BiomWRapeVol

DicotWeed = BiomCABP+BiomCHEAL+BiomMYOAR+BiomSTEME+BRACAhybtot+BiomSTEMESU

DicotWeedVol = DicotWeed+DicotVol

MonocotWeedVol = BiomELYRE+BiomWBarleyVol+BiomWWheatVol

WeedsVolTotal = DicotWeedVol+MonocotWeedVol

WeedVolNoOSR = DicotWeed+MonocotWeedVol

BiomSTEMESU(t) = BiomSTEMESU(t - dt) + (GrowSTEMESU + SowSTEMESU - HarvSTEMESU - ControlSTEMESU) * dt

INIT BiomSTEMESU = 0

GrowSTEMEU = if (MaxSteme>BiomTotal) THEN ((Winter*RGRSTEME*BiomSTEMESU*(MaxSteme-BiomTotal))/MaxSteme)/dt else 0

SowSTEMESU = Delay(BeginWeed,1)*(SeedSTEME*Mutation_Frq+SeedSTEMESU)*10000*GermSTEME*TKVSTEME/(1000*1000000)/dt

HarvSTEMESU = if EndWeed=1 then (BiomSTEMESU+0.1) else 0

ControlSTEMESU = BiomSTEMESU*EFSTEMESU

SeedSTEMESU(t) = SeedSTEMESU(t - dt) + (NewSeedSTEMESU - DeadSeedSTEMESU) * dt

INIT SeedSTEMESU = 0

NewSeedSTEMESU = (HarvSTEMESU)*NonPredSTEME*HISTEME*1000000*1000/(10000*TKVSTEME)/dt

DeadSeedSTEMESU = (SeedSTEMESU*(1-EXP(LOGN(1-DecaySteme)/365))+BeginWeed*SeedSTEMESU*GermSTEME)/dt

Mutation_Frq = 1/10000

RotationYear(t) = RotationYear(t - dt) + (RotationYearRate - RotationYearReset) * dt

INIT RotationYear = 0

RotationYearRate = if DayInYear=365 then 1/DT else 0

RotationYearReset = IF(RotationYear)=4 and DayInYear=365 then (RotationYear)/DT else 0

BeginWBarley = if WBarleynr=1 and DayInYear=259 then 1 else 0

BeginWeed = if BeginWBarley=1 or BeginWRape=1 or BeginWWheat=1 then 1 else 0

BeginWRape = IF(WRapenr=1) and DayInYear=228 then 1 else 0

BeginWWheat = If (WWheatnr=1) and DayInYear=280 then 1 else 0

DayInYear = COUNTER(1,366)

EndWBarley = DELAY(BeginWBarley,308)

EndWeed = if EndWBarley=1 or EndWRape=1 or EndWWheat=1 then 1 else 0

EndWRape = DELAY(BeginWRape,351)

EndWWheat = DELAY(BeginWWheat,308)

SumAnn = If RotationYear>0 and DayInYear=76 then 1 else 0

WBarleynr = if (RotationYear=3) then 1 else 0

Winter = IF DayInYear<90 or DayInYear>320 then 0 else 1

WRapenr = IF (RotationYear=0 OR RotationYear=4) then 1 else 0

WWheatnr = if (RotationYear=1 or RotationYear=2) then 1 else 0

BiomWBarleyVol(t) = BiomWBarleyVol(t - dt) + (GrowWBarley_2 + SowWBarleyVol - HarvWBarleyVol - ControlWBarleyVol) * dt

INIT BiomWBarleyVol = 0

GrowWBarley_2 = if (MaxWBarley>BiomTotal) THEN (Winter*RGRWBarley*BiomWBarleyVol*(MaxWBarley-BiomTotal)/MaxWBarley)/dt else 0

SowWBarleyVol = (Delay(BeginWeed,1)*SeedWBarleyVol *GermWBarleyVol*TKWWBarley/(1000*1000000))/dt

HarvWBarleyVol = if EndWeed=1 then (BiomWBarleyVol+0.1)/dt ele 0

ControlWBarleyVol = BiomWBarleyVol*EFWBarley

BiomWRapeVol(t) = BiomWRapeVol(t - dt) + (GrowWRapeVol + SowWRapeVol - HarvWRapeVol - ControlWRapeVol) * dt

INIT BiomWRapeVol = 0

GrowWRapeVol = if (MaxWRape>BiomTotal) THEN (Winter*RGRWRape*BiomWRapeVol*(MaxWRape-BiomTotal)/MaxWRape)/dt else 0

SowWRapeVol = (Delay(BeginWeed,1)*SeedWRapeVol*GermWRapeVol*TKWWRape/(1000*1000000))/dt

HarvWRapeVol = if EndWeed=1 then (BiomWRapeVol+0.1)/dt else 0

ControlWRapeVol = BiomWRapeVol*EFBRARe

BiomWWheatVol(t) = BiomWWheatVol(t - dt) + (GrowWWheatVol + SowWWheatVol - HarvWWheatVol - ControlWWheatVol) * dt

INIT BiomWWheatVol = 0

GrowWWheatVol = if (MaxWWheat>BiomTotal) THEN (Winter*RGRWWheat*BiomWWheatVol*(MaxWWheat-BiomTotal)/MaxWWheat)/dt else 0

SowWWheatVol = (Delay(BeginWeed,1)*SeedWWheatVol *GermWWheatVol*TKWWWheat/(1000*1000000))/dt

HarvWWheatVol = IF (EndWeed)=1 THEN (BiomWWheatVol+0.1)/dt ELSE(0)

ControlWWheatVol = BiomWWheatVol*EFWWheat

SeedWBarleyVol(t) = SeedWBarleyVol(t - dt) + (NewSeedWBarleyVol - DeadSeedWBarleyVol) * dt

INIT SeedWBarleyVol = 0

NewSeedWBarleyVol = (((HarvWBarley+HarvWBarleyVol)*LossWBarley*DormWBarley*HIWBarley*1000000*1000)/(TKWWBarley))/dt

DeadSeedWBarleyVol = (SeedWBarleyVol*(1-EXP(LOGN(1-DecayWBarley)/365))+BeginWeed*SeedWBarleyVol*GermWBarleyVol)/dt

SeedWRapeVol(t) = SeedWRapeVol(t - dt) + (NewSeedWRapeVol - DeadSeedWRapeVol) * dt

INIT SeedWRapeVol = 0

NewSeedWRapeVol = IF (HarvWRape>0.01 AND Bn_Rr?=0) THEN ((HarvWRape+HarvWRapeVol)*LossWRape*DormWRape*HIWRape*1000000*1000/TKWWRape)/dt

ELSE

((-1)*(Bn_Rr?-1)* (HarvWRape*LossWRape+HarvWRapeVol)*DormWRape*HIWRape*1000000*1000/TKWWRape)/dt

DeadSeedWRapeVol = (SeedWRapeVol*(1-EXP(LOGN(1-DecayWRape)/365))+BeginWeed*SeedWRapeVol*GermWRapeVol)/dt

SeedWWheatVol(t) = SeedWWheatVol(t - dt) + (NewSeedWWheatVol - DeadSeedWWheatVol) * dt

INIT SeedWWheatVol = 0

NewSeedWWheatVol = ((HarvWWheat+HarvWWheatVol)*LossWwheat*DormWWheat*HIWWheat*1000000*1000)/(TKWWWheat)/dt

DeadSeedWWheatVol = (SeedWWheatVol*(1-EXP(LOGN(1-DecayWWheat)/365))+BeginWeed*SeedWWhetVol*GermWWheatVol)/dt

DecayWBarley = 0.99999

DecayWRape = 0.90

DecayWWheat = 0.99999

DormWBarley = 0.15

DormWRape = 0.025

DormWWheat = 0.15

GermWBarleyVol = 0.33

GermWRapeVol = 0.20

GermWWheatVol = 0.33

LossWBarley = 0.10

LossWRape = 0.10

LossWwheat = 0.1

TKWWBarley = 44

TKWWRape = 5.3

TKWWWheat = 44

AccCropYield(t) = AccCropYield(t - dt) + (AnnCropYield) * dt

INIT AccCropYield = 0

AnnCropYield = BiomCrop*EndWeed/dt

AccTotalYield(t) = AccTotalYield(t - dt) + (AnnTotalYield) * dt

INIT AccTotalYield = 0

AnnTotalYield = BiomTotal*EndWeed/dt

Scen125acc(t) = Scen125acc(t - dt) + (SprayScen125) * dt

INIT Scen125acc = 0

SprayScen125 = IF (BiomWRape>0 AND (Scenario?=1 OR Scenario?=2 OR Scenario?=5) AND Winter=1 AND MonocotWeedVol>0.5) THEN (WeedControl?*DoseFluazifop) ELSE (0)

Scen125TFq(t) = Scen125TFq(t - dt) + (Scen125NrRate) * dt

INIT Scen125TFq = 0

Scen125NrRate = If SprayScen125>0 THEN 1 ELSE 0

Scen1acc(t) = Scen1acc(t - dt) + (SprayScen1) * dt

INIT Scen1acc = 0

SprayScen1 = IF (BiomWRape>0 AND Scenario?=1 AND DayInYear=321 AND WeedVolNoOSR>0.5) THEN (WeedControl?*DosePropamid) ELSE (0)

Scen1TFq(t) = Scen1TFq(t - dt) + (Scen1NrRate) * dt

INIT Scen1TFq = 0

Scen1NrRate = If SprayScen1>0 THEN 1 ELSE 0

Scen2acc(t) = Scen2acc(t - dt) + (SprayScen2) * dt

INIT Scen2acc = 0

SprayScen2 = IF (BiomWRape>0 AND Scenario?=2 AND Winter=1 AND (DayInYear<165 OR DayInYear>228) AND DicotWeed>1) THEN (WeedControl?*DoseBenasa) ELSE (0)

Scen2TFq(t) = Scen2TFq(t - dt) + (Scen2NrRate) * dt

INIT Scen2TFq = 0

Scen2NrRate = If SprayScen2>0 THEN 1 ELSE 0

Scen3acc(t) = Scen3acc(t - dt) + (SprayScen3) * dt

INIT Scen3acc = 0

SprayScen3 = IF (BiomWRape>0 AND Scenario?=3 AND Winter=1 AND (DayInYear<165 OR DayInYear>228) AND WeedVolNoOSR>1.4) THEN (WeedControl?*DoseGlyphosate) ELSE (0)

Scen3TFq(t) = Scen3TFq(t - dt) + (Scen3NrRate) * dt

INIT Scen3TFq = 0

Scen3NrRate = If SprayScen3>0 THEN 1 ELSE 0

Scen4acc(t) = Scen4acc(t - dt) + (SprayScen4) * dt

INIT Scen4acc = 0

SprayScen4 = IF (BiomWRape>0 AND Scenario?=4 AND Winter=1 AND (DayInYea< 165 OR DayInYear>228) AND WeedVolNoOSR>1.4) THEN (WeedControl?*DoseGlufosinate) ELSE (0)

Scen4TFq(t) = Scen4TFq(t - dt) + (Scen4NrRate) * dt

INIT Scen4TFq = 0

Scen4NrRate = If SprayScen4>0 THEN 1 ELSE 0

Scen5acc(t) = Scen5acc(t - dt) + (SprayScen5) * dt

INIT Scen5acc = 0

SprayScen5 = IF (BiomWRape>0 AND Scenario?=5 AND Winter=1 AND Scen5TFq_2<2.9 AND (DayInYear<165 OR DayInYear>228) AND DicotWeed>0.8) THEN (WeedControl?*DoseSUWRape) ELSE (0)

Scen5TFq(t) = Scen5TFq(t - dt) + (Scen5NrRate) * dt

INIT Scen5TFq = 0

Scen5NrRate = If SprayScen5>0 THEN 1 ELSE 0

Scen5TFq_2(t) = Scen5TFq_2(t - dt) + (Scen5NrRate_2 - Reset5TFq) * dt

INIT Scen5TFq_2 = 0

Scen5NrRate_2 = If SprayScen5>0 THEN 1 ELSE 0

Reset5TFq = IF DayInYear=210 THEN Scen5TFq_2 ELSE 0

ScenCDicotacc(t) = ScenCDicotacc(t - dt) + (SprayCDicot) * dt

INIT ScenCDicotacc = 0

SprayCDicot = IF (BiomCereals>0 AND Winter=1 AND DicotWeedVol>0.1) THEN (WeedControl?*DoseSU) ELSE (0)

ScenCDicotTFq(t) = ScenCDicotTFq(t - dt) + (ScenCDicotNrRate) * dt

INIT ScenCDicotTFq = 0

ScenCDicotNrRate = If SprayCDicot>0 THEN 1 ELSE 0

ScenCELYRETFq(t) = ScenCELYRETFq(t - dt) + (ScenELYRENrRate) * dt

INIT ScenCELYRETFq = 0

ScenELYRENrRate = If SprayELYRE>0 THEN 1 ELSE 0

ScenCFLuMixacc(t) = ScenCFLuMixacc(t - dt) + (SprayCFluMix) * dt

INIT ScenCFLuMixacc = 0

SprayCFluMix = IF (BiomCereals>0 AND Winter=1 AND DicotWeedVol>0.1) THEN (WeedControl?*Dosefluroxypyr) ELSE (0)

ScenCFluroxypyrTFq(t) = ScenCFluroxypyrTFq(t - dt) + (ScenFluroxypyr) * dt

INIT ScenCFluroxypyrTFq = 0

ScenFluroxypyr = If SprayCFluMix>0 THEN 1 ELSE 0

ScenELYREacc(t) = ScenELYREacc(t - dt) + (SprayELYRE) * dt

INIT ScenELYREacc = 0

SprayELYRE = IF (BiomCereals>0 AND EndWeed=1 AND BiomELYRE>0.3) THEN (WeedControl?*DoseGlyphosatePH) ELSE (0)

BRACAResistance = IF ((SprayCDicot>0 AND Scenario?=5) OR SprayScen1>0 OR SprayScen2>0 OR SprayScen3>0 OR SprayScen4>0 OR SprayScen5>0) AND (BiomBRACARR+BiomBRACARS+BiomF1)/DicotWeedVol>0.25 THEN 1 ELSE 0

CerealTFq = ScenCDicotTFq+ScenCELYRETFq+ScenCFluroxypyrTFq

CerealUse = ScnCDicotacc+ScenCFLuMixacc+ScenELYREacc

Dose = 1

DoseBenasa = Dose*0.406

DoseFluazifop = 0.188

Dosefluroxypyr = Mix*0.144

DoseGlufosinate = Dose*0.600

DoseGlyphosate = Dose*0.445

DoseGlyphosatePH = 0.8

DosePropamid = 0.5

DoseSU = 0.0075*Dose

DoseSUWRape = 0.0075*Dose

dummy = IF (Scenario?=5) THEN 0 ELSE 1

Ef125ELYRE = IF SprayScen125>0 THEN 0.85 ELSE 0

Ef125WBarleyVol = IF SprayScen125>0 THEN 0.95 ELSE 0

Ef125WWheatVol = IF SprayScen125>0 THEN 0.95 ELSE 0

Ef1CAPB = IF SprayScen1>0 THEN 0.2 ELSE 0

Ef1CHEAL = IF SprayScen1>0 THEN 0.75 ELSE 0

Ef1ELYRE = IF SprayScen1>0 THEN 0.5 ELSE 0

Ef1MYOAR = IF SprayScen1>0 THEN 0.75 ELSE 0

Ef1STEME = IF SprayScen1>0 THEN 0.95 ELSE 0

Ef1WBarleyVol = IF SprayScen1>0 THEN 0.95 ELSE 0

Ef1WWheatVol = IF SprayScen1>0 THEN 0.95 ELSE 0

Ef2CAPB = IF SprayScen2>0 THEN (1-(0.28+(1-0.28)/(1+EXP(1.3*(LOGN(SprayScen2)-LOGN(0.600)))))) ELSE 0

Ef2CHEAL = IF SprayScen2>0 THEN (1-(0.25+(1-0.25)/(1+EXP(3.5*(LOGN(SprayScen2)-LOGN(0.3166)))))) ELSE 0

Ef2MYOAR = IF SprayScen2>0 THEN (1-(0.28+(1-0.28)/(1+EXP(1.3*(LOGN(SprayScen2)-LOGN(0.150)))))) ELSE 0

Ef2STEME = IF SprayScen2>0 THEN (1-(0.28+((1-0.28)/(1+EXP(1.3*(LOGN(SprayScen2)-LOGN(0.0839))))))) ELSE 0

Ef3BRASe = IF SprayScen3>0 THEN (1-(0.02+(1-0.02)/(1+EXP(3.1*(LOGN(SprayScen3)-LOGN(0.2067)))))) ELSE 0

Ef3CAPB = IF SprayScen3>0 THEN (1-(0.02+(1-0.02)/(1+EXP(1.9*(LOGN(SprayScen3)-LOGN(0.0485)))))) ELSE 0

Ef3CHEAL = IF SprayScen3>0 THEN (1-(0.02+(1-0.02)/(1+EXP(2.4*(LOGN(SprayScen3)-LOGN(0.1058)))))) ELSE 0

Ef3ELYRE = IF SprayScen3>0 THEN (1-(0.02+(1-0.02)/(1+EXP(1.9*(LOGN(SprayScen3)-LOGN(0.150)))))) ELSE 0

Ef3MYOAR = IF SprayScen3>0 THEN (1-(0.02+(1-0.02)/(1+EXP(1.2*(LOGN(SprayScen3)-LOGN(0.0824)))))) ELSE 0

Ef3STEME = IF SprayScen3>0 THEN (1-(0.02+(1-0.02)/(1+EXP(1.9*(LOGN(SprayScen3)-LOGN(0.073)))))) ELSE 0

Ef3WBarleyVol = IF SprayScen3>0 THEN (1-(0.02+(1-0.02)/(1+EXP(1.9*(LOGN(SprayScen3)-LOGN(0.073)))))) ELSE 0

Ef3WWheatVol = IF SprayScen3>0 THEN (1-(0.02+(1-0.02)/(1+EXP(1.9*(LOGN(SprayScen3)-LOGN(0.073)))))) ELSE 0

Ef4BRASe = IF SprayScen>0 THEN (1-(0.1+(1-0.1)/(1+EXP(3*(LOGN(SprayScen4)-LOGN(0.3715)))))) ELSE 0

Ef4CAPB = IF SprayScen4>0 THEN (1-(0.02+(1-0.02)/(1+EXP(1.7*(LOGN(SprayScen4)-LOGN(0.064)))))) ELSE 0

Ef4CHEAL = IF SprayScen4>0 THEN (1-(0.02+(1-0.02)/(1+EXP(3.5*(LOGN(SprayScen4)-LOGN(0.1613)))))) ELSE 0

Ef4ELYRE = IF SprayScen4>0 THEN (1-(0.25+(1-0.25)/(1+EXP(3*(LOGN(SprayScen4)-LOGN(0.3715)))))) ELSE 0

Ef4MYOAR = IF SprayScen4>0 THEN (1-(0.02+(1-0.02)/(1+EXP(2.7*(LOGN(SprayScen4)-LOGN(0.0573)))))) ELSE 0

Ef4STEME = IF SprayScen4>0 THEN (1-(0.02+(1-0.02)/(1+EXP(2.9*(LOGN(SprayScen4)-LOGN(0.0653)))))) ELSE 0

Ef4WBarleyVol = IF SprayScen4>0 THEN (1-(0.02+(1-0.02)/(1+EXP(2.9*(LOGN(SprayScen4)-LOGN(0.1706)))))) ELSE 0

Ef4WWheatVol = IF SprayScen4>0 THEN (1-(0.02+(1-0.02)/(1+EXP(2.9*(LOGN(SprayScen4)-LOGN(0.1706)))))) ELSE 0

EFBRARe = EfCRapeVol+EFCFluMix*(1-EfCRapeVol)

EFBRASe = Ef3BRASe+Ef4BRASe+Ef5BraSe+EfCBRASe+EFCFluMix*(1-Ef3BRASe+Ef4BRASe+Ef5BraSe+EfCBRASe)

EFCAPBP = Ef1CAPB+Ef2CAPB+Ef3CAPB+Ef4CAPB+Ef5CAPB+EfCCAPB+EFCFluMix*(1-Ef1CAPB+Ef2CAPB+Ef3CAPB+Ef4CAPB+Ef5CAPB+EfCCAPB)

EFCFluMix = IF SprayCFluMix>0 THEN 0.95 ELSE 0

EFCHEAL = Ef1CHEAL+Ef2CHEAL+Ef3CHEAL+Ef4CHEAL+Ef5CHEAL+EfCCHEAL+EFCFluMix*(1-Ef1CHEAL+Ef2CHEAL+Ef3CHEAL+Ef4CHEAL+Ef5CHEAL+EfCCHEAL)

EFELYRE = Ef125ELYRE+Ef1ELYRE+Ef3ELYRE+Ef4ELYRE+EfGlyPHELYRE

EfGlyPHELYRE = IF SprayELYRE>0 THEN 0.95 ELSE 0

EFMYOAR = Ef1MYOAR+Ef2MYOAR+Ef3MYOAR+Ef4MYOAR+Ef5MYOAR+EfCMYOAR+EFCFluMix*(1-Ef1MYOAR+Ef2MYOAR+Ef3MYOAR+Ef4MYOAR+Ef5MYOAR+EfCMYOAR)

EFSTEME = Ef1STEME+Ef2STEME+Ef3STEME+Ef4STEME+Ef5STEME+EfCSTEME

+EFCFluMix*(1-Ef1STEME+Ef2STEME+Ef3STEME+Ef4STEME+Ef5STEME+EfCSTEME)

EFSTEMESU = Ef1STEME+Ef2STEME+Ef3STEME+Ef4STEME+EFCFluMix

EFWBarley = Ef125WBarleyVol+Ef1WBarleyVol+Ef3WBarleyVol+Ef4WBarleyVol

EFWWheat = Ef125WWheatVol+Ef1WWheatVol+Ef3WWheatVol+Ef4WWheatVol

Mix = IF (BRACAResistance=1 OR SUResistance=1 OR WRapeVolResistance=1) THEN 1 ELSE 0

RapeUse = Scen125acc+Scen1acc+Scen2acc+Scen3acc+Scen4acc+Scen5acc

Scenario? = 2

SUResistance = IF (SprayCDicot>0 OR SprayScen1>0 OR SprayScen2>0OR SprayScen3>0 OR SprayScen4>0 OR SprayScen5>0) AND BiomSTEMESU/DicotWeed>0.25 THEN 1 ELSE 0

TotalTFq = CerealTFq+WRapeTFq

TotalUse = CerealUse+RapeUse

WeedControl? = 1

WRapeTFq = Scen125TFq+Scen1TFq+Scen2TFq+Scen3TFq+Scen4TFq+Scen5TFq

WRapeVolResistance = IF ((SprayCDicot>0 AND Scenario?=5) OR SprayScen1>0 OR SprayScen2>0 OR SprayScen3>0 OR SprayScen4>0 OR SprayScen5>0) AND (BiomWRapeVol+BiomWRapeVolReSe+BiomWRapeVolRe)/DicotWeedVol>0.25 THEN 1 ELSE 0

Yieldreduction = IF AccCropYield>0 THEN (AccTotalYield-AccCropYield)/AccTotalYield ELSE 0

Ef5BraSe = GRAPH(SprayScen5)

(0.00, 0.00), (0.0015, 0.895), (0.003, 0.93), (0.0045, 0.95), (0.006, 0.96), (0.0075, 0.97), (0.009, 0.975), (0.0105, 0.975), (0.012, 0.98), (0.0135, 0.98), (0.015, 0.98)

Ef5CAPB = GRAPH(SprayScen5)

(0.00, 0.00), (0.0015, 0.575), (0.003, 0.7), (0.0045, 0.755), (0.006, 0.8), (0.0075, 0.83), (0.009, 0.86), (0.0105, 0.875), (0.012, 0.89), (0.0135, 0.89), (0.015, 0.89)

Ef5CHEAL = GRAPH(SprayScen5)

(0.00, 0.00), (0.0015, 0.615), (0.003, 0.725), (0.0045, 0.775), (0.006, 0.82), (0.0075, 0.85), (0.009, 0.875), (0.0105, 0.895), (0.012, 0.905), (0.0135, 0.91), (0.015, 0.91)

Ef5MYOAR = GRAPH(SprayScen5)

(0.00, 0.00), (0.0015, 0.625), (0.003, 0.73), (0.0045, 0.77), (0.006, 0.815), (0.0075, 0.85), (0.009, 0.87), (0.0105, 0.885), (0.012, 0.9), (0.0135, 0.91), (0.015, 0.91)

Ef5STEME = GRAPH( SprayScen5)

(0.00, 0.00), (0.0015, 0.835), (0.003, 0.895), (0.0045, 0.925), (0.006, 0.945), (0.0075, 0.95), (0.009, 0.96), (0.0105, 0.965), (0.012, 0.965), (0.0135, 0.97), (0.015, 0.97)

EfCBRASe = GRAPH( SprayCDicot)

(0.00, 0.00), (0.0015, 0.895), (0.003, 0.93), (0.0045, 0.95), (0.006, 0.965), (0.0075, 0.97), (0.009, 0.97), (0.0105, 0.975), (0.012, 0.98), (0.0135, 0.98), (0.015, 0.98)

EfCCAPB = GRAPH( SprayCDicot)

(0.00, 0.00), (0.0015, 0.575), (0.003, 0.7), (0.0045, 0.74), (0.006, 0.78), (0.0075, 0.82), (0.009, 0.845), (0.0105, 0.865), (0.012, 0.885), (0.0135, 0.89), (0.015, 0.89)

EfCCHEAL = GRAPH( SprayCDicot)

(0.00, 0.00), (0.0015, 0.6), (0.003, 0.705), (0.0045, 0.77), (.006, 0.83), (0.0075, 0.85), (0.009, 0.87), (0.0105, 0.885), (0.012, 0.895), (0.0135, 0.91), (0.015, 0.91)

EfCMYOAR = GRAPH( SprayCDicot)

(0.00, 0.00), (0.0015, 0.61), (0.003, 0.72), (0.0045, 0.77), (0.006, 0.81), (0.0075, 0.85), (0.009, 0.88), (0.0105, 0.89), (0.012, 0.9), (0.0135, 0.91), (0.015, 0.91)

EfCRapeVol = GRAPH(dummy*SprayCDicot)

(0.00, 0.00), (0.0015, 0.89), (0.003, 0.925), (0.0045, 0.95), (0.006, 0.96), (0.0075, 0.97), (0.009, 0.97), (0.0105, 0.975), (0.012, 0.975), (0.0135, 0.98), (0.015, 0.98)

EfCSTEME = GRAPH(SprayCDicot)

(0.00, 0.00), (0.0015, 0.825), (0.003, 0.895), (0.0045, 0.925), (0.006, 0.95), (0.0075, 0.955), (0.009, 0.96), (0.0105, 0.965), (0.012, 0.97), (0.0135, 0.975), (0.015, 0.975)

BiomBRACA(t) = BiomBRACA(t - dt) + (GrowBRACA + SowBRACA - HarvBRACA - ControlBRACA) * dt

INIT BiomBRACA = 0

GrowBRACA = if (MaxBRACA>BiomTotal AND Hybridization=0) THEN (Winter*RGRBRACA*BiomBRACA*(MaxBRACA-BiomTotal)/MaxBRACA)/dt else 0

SowBRACA = Delay(BeginWeed,1)*SeedBRACA*GermBRACA*TKVBRACA*10000/(1000*1000000)/dt

HarvBRACA = if EndWeed=1 then (BiomBRACA+0.1) else 0

ControlBRACA = BiomBRACA*EFBRASe

BiomCABP(t) = BiomCABP(t - dt) + (GrowCABP + SowCABP - HarvCABP - ControlCABP) * dt

INIT BiomCABP = 0

GrowCABP = if (MaxCABP>BiomTotal) THEN (Winter*RGRCABP*BiomCABP*(MaxCABP-BiomTotal)/MaxCABP)/dt else 0

SowCABP = Delay(BeginWeed,1)*SeedCABP*GermCABP*TKVCABP*10000/(1000*1000000)/dt

HarvCABP = if EndWeed=1 then (BiomCABP+0.1) else 0

ControlCABP = BiomCABP*EFCAPBP

BiomCHEAL(t) = BiomCHEAL(t - dt) + (GrowCHEAL + SowCHEAL - HarvCHEAL - ControlCHEAL) * dt

INIT BiomCHEAL = 0

GrowCHEAL = if (MaxCHEAL>BiomTotal) THEN (Winter*RGRCHEAL*BiomCHEAL*(MaxCHEAL-BiomTotal)/MaxCHEAL)/dt else 0

SowCHEAL = if SumAnn=1 then SeedCHEAL*10000*GermCHEAL*TKVCHEAL/(1000*1000000)/dt else 0

HarvCHEAL = if EndWeed=1 then (BiomCHEAL+0.1) else 0

ControlCHEAL = BiomCHEAL*EFCHEAL

BiomELYRE(t) = BiomELYRE(t - dt) + (GrowELYRE + InitELYRE - HarvELYRE - ControlELYRE) * dt

INIT BiomELYRE = 0

GrowELYRE = if (MaxELYRE>BiomTotal) THEN (Winter*GRELYRE*BiomELYRE*(MaxELYRE-BiomTotal)/MaxELYRE)/dt else 0

InitELYRE = IF (TIME<365) THEN (ELYRE?*BeginWeed*0.005) ELSE (0)

HarvELYRE = DELAY(EndWeed,1)*BiomELYRE*0.975

ControlELYRE = BiomELYRE*EFELYRE

BiomMYOAR(t) = BiomMYOAR(t - dt) + (GrowMYOAR + SowMYOAR - HarvMYOAR - ControlMYOAR) * dt

INIT BiomMYOAR = 0

GrowMYOAR = if (MaxMYOAR>BiomTotal) THEN (Winter*RGRMYOAR*BiomMYOAR*(MaxMYOAR-BiomTotal)/MaxMYOAR)/dt else 0

SowMYOAR = Delay(BeginWeed,1)*SeedMYOAR*GermMYOAR*TKVMYOAR*10000/(1000*1000000)/dt

HarvMYOAR = if EndWeed=1 then (BiomMYOAR+0.1) else 0

ControlMYOAR = BiomMYOAR*EFMYOAR

BiomSTEME(t) = BiomSTEME(t - dt) + (GrowSTEME + SowSTEME - HarvSTEME - ControlSTEME) * dt

INIT BiomSTEME = 0

GrowSTEME = if (MaxSteme>BiomTotal) THEN (Winter*RGRSTEME*BiomSTEME*(MaxSteme-BiomTotal)/MaxSteme)/dt else 0

SowSTEME = Delay(BeginWeed,1)*SeedSTEME*10000*GermSTEME*TKVSTEME

/(1000*1000000)/dt

HarvSTEME = if EndWeed=1 then (BiomSTEME+0.1) else 0

ControlSTEME = BiomSTEME*EFSTEME

SeedBRACA(t) = SeedBRACA(t - dt) + (NewSeedBRACA - DeadSeedBRACA) * dt

INIT SeedBRACA = 500

NewSeedBRACA = (HarvBRACA)*NonPredBRACA*HIBRACA*1000000*1000/(TKVBRACA*10000)/dt

DeadSeedBRACA = If Hybridization=0 THEN (SeedBRACA*(1-EXP(LOGN(1-DecayBRACA)/365))+BeginWeed*SeedBRACA*GermBRACA)/dt

ELSE SeedBRACA

SeedCABP(t) = SeedCABP(t - dt) + (NewSeedCABP - DeadSeedCABP) * dt

INIT SeedCABP = 3600

NewSeedCABP = (HarvCABP)*NonPredCABP*HICABP*1000000*1000/(TKVCABP*10000)/dt

DeadSeedCABP = (SeedCABP*(1-EXP(LOGN(1-DecayCABP)/365))+BeginWeed*SeedCABP*GermCABP)/dt

SeedCHEAL(t) = SeedCHEAL(t - dt) + (NewSeedCHEAL - DeadSeedCHEAL) * dt

INIT SeedCHEAL = 1000

NewSeedCHEAL = (HarvCHEAL)*NonPredCHEAL*HICHEAL*1000000*1000/(TKVCHEAL*10000*dt)

DeadSeedCHEAL = (SeedCHEAL*(1-EXP(LOGN(1-DecayCHEAL)/365))+BeginWeed*SeedCHEAL*GermCHEAL)/dt

SeedMYOAR(t) = SeedMYOAR(t - dt) + (NewSeedMYOAR - DeadSeedMYOAR) * dt

INIT SeedMYOAR = 3400

NewSeedMYOAR = (HarvMYOAR)*NonPredMYOAR*HIMYOAR*1000000*1000/(TKVMYOAR*10000)/dt

DeadSeedMYOAR = (SeedMYOAR*(1-EXP(LOGN(1-DecayMYOAR)/365))+Begineed*SeedMYOAR*GermMYOAR)/dt

SeedSTEME(t) = SeedSTEME(t - dt) + (NewSeedSTEME - DeadSeedSTEME) * dt

INIT SeedSTEME = 10800

NewSeedSTEME = (HarvSTEME)*NonPredSTEME*HISTEME*1000000*1000/(10000*TKVSTEME)/dt

DeadSeedSTEME = (SeedSTEME*(1-EXP(LOGN(1-DecaySteme)/365))+BeginWeed*SeedSTEME*GermSTEME)/dt

DecayBRACA = 0.33

DecayCABP = 0.33

DecayCHEAL = 0.33

DecayMYOAR = 0.33

DecaySteme = 0.33

ELYRE? = 1

GermBRACA = 0.05

GermCABP = 0.05

GermCHEAL = 0.05

GermMYOAR = 0.03

GermSTEME = 0.05

HIBRACA = 0.25

HICABP = 0.25

HICHEAL = 0.5

HIMYOAR = 0.25

HISTEME = 0.5

MaxBRACA = 8

MaxCABP = 5

MaxCHEAL = 8

MaxELYRE = 8

MaxMYOAR = 5

MaxSteme = 5

NonPredBRACA = 0.8

NonPredCABP = 0.85

NonPredCHEAL = 0.5

NonPredMYOAR = 0.6

NonPredSTEME = 0.5

RGRBRACA = 0.061

RGRCABP = 0.07

RGRCHEAL = 0.15

RGRELYRE = 0.028

RGRMYOAR = 0.08

RGRSTEME = 0.067

TKVBRACA = 2

TKVCABP = 0.1

TKVCHEAL = 1.2

TKVMYOAR = 0.3