[MAIN]
if (NEWIND != 2) NFOOD=FOOD
KAFD=NFOOD*THETA_NFOOD + 1
LAMB=THETA_LAMB
IMAX=THETA_IMAX
CLBW=pow((1.0/70.0)*BBW, THETA_CLBW)
VDBW=pow((1.0/70.0)*BBW, THETA_VDBW)
ST99=0
if (STUD == 99) ST99=1
CL99=ST99*THETA_CL99 + 1
VIS1=0
VIS2=0
VIS3=0
VIS4=0
VIS5=0
VIS6=0
VIS7=0
VIS8=0
VIS9=0
VIS10=0
if (OCC == 1) VIS1=1
if (OCC == 2) VIS2=1
if (OCC == 3) VIS3=1
if (OCC == 4) VIS4=1
if (OCC == 5) VIS5=1
if (OCC == 6) VIS6=1
if (OCC == 7) VIS7=1
if (OCC == 8) VIS8=1
if (OCC <= 14 && OCC > 8) VIS9=1
if (OCC > 14) VIS10=1
IOVF1A=ETA_5*VIS1 + ETA_6*VIS2 + ETA_7*VIS3 + ETA_8*VIS4 + ETA_9*VIS5
IOVF1B=ETA_10*VIS6 + ETA_11*VIS7 + ETA_12*VIS8 + ETA_13*VIS9 + ETA_14*VIS10
IOVF1=IOVF1A + IOVF1B
IOVKA1=ETA_15*VIS1 + ETA_16*VIS2 + ETA_17*VIS3 + ETA_18*VIS4 + ETA_19*VIS5
IOVKA2=ETA_20*VIS6 + ETA_21*VIS7 + ETA_22*VIS8 + ETA_23*VIS9 + ETA_24*VIS10
IOVKA=IOVKA1 + IOVKA2
CL=CL99*CLBW*THETA_CL*exp(ETA_1)
V2=THETA_V2*VDBW
Q3=CLBW*THETA_Q3
V3=THETA_V3*VDBW
KA=KAFD*THETA_KA*exp(ETA_3 + IOVKA)
F1TM=-IMAX + IMAX*exp(-LAMB*TIME) + 1
F1=F1TM*THETA_F1*exp(ETA_2 + IOVF1)
K20=CL/V2
K23=Q3/V2
K32=Q3/V3
KTR=KA
S2=V2
NFOOD=FOOD
DONO1=0
if (DONO == 1) DONO1=1

[ODE]
CP=A_2/V2
d/dt(A_1)=-A_1*KA
d/dt(A_2)=-A_2*K20 - A_2*K23 + A_3*K32 + A_7*KTR
d/dt(A_3)=A_2*K23 - A_3*K32
d/dt(A_4)=A_1*KA - A_4*KTR
d/dt(A_5)=A_4*KTR - A_5*KTR
d/dt(A_6)=A_5*KTR - A_6*KTR
d/dt(A_7)=A_6*KTR - A_7*KTR

[ERROR]
CONC=A_2/S2
IPRED=-9.2103403719761836
if (CONC > 0) IPRED=log(CONC)
W=THETA_W*exp(ETA_4)
Y=EPS_ADD*W + IPRED
IRES=-DV + IPRED
IWRES=IRES/W

