#muraro1



#reactions
r1:
 $pool = {1.0}IAAm
  phiIAAp * (lambda1 * F1 + F2 + lambda3 * F3) - IAAm
r10:
 $pool = {1.0}PINp
  1 / eps * (deltaPINp * PINm - PINp)
r11:
 $pool = {1.0}ARm
  phiARp * (lambda1 * F5a + F5b) - ARm
r12:
 $pool = {1.0}ARp
  1 / eps * (deltaARp * ARm - ARp)
r13:
 $pool = {1.0}CRm
  phiCRp * F4 - CRm
r14:
 $pool = {1.0}CRp
  1 / eps * (deltaCRp * CRm - CRp)
r15:
 $pool = {1.0}AHKph
  1 / eps * (rd * CkAHKph - ra * AHKph * Ck)
r16:
 $pool = {1.0}Ck
  muCk * (alphaCk - Ck) - etaCkPh / eps * (ra * AHKph * Ck - rd * CkAHKph)
r17:
 $pool = {1.0}ARRBph
  1 / eps * (ua * CkAHKph * ARRBp - ud * CkAHK * ARRBph)
r18:
 $pool = {1.0}ARRAph
  1 / eps * (sa * CkAHKph * ARRAp - sd * CkAHK * ARRAph)
r19:
 $pool = {1.0}ARRAm
  phiARRAp * F6 - ARRAm
r2:
 $pool = {1.0}IAAp
  1 / eps * (deltaIAAp * IAAm - la * IAAp * AuxTIR1 + ld * AuxTIAA) + etaARFIAA * (pd * ARFIAA - pa * IAAp * ARF)
r20:
 $pool = {1.0}ARRAp
  1 / eps * (deltaARRAp * ARRAm - ARRAp + etaAHKph * (sd * CkAHK * ARRAph - sa * CkAHKph * ARRAp))
r3:
 $pool = {1.0}AuxTIR1
  1 / eps * (ka * Aux * TIR1 - kd * AuxTIR1 + (ld + 1) * AuxTIAA - la * AuxTIR1 * IAAp)
r4:
 $pool = {1.0}AuxTIAA
  1 / eps * (la * IAAp * AuxTIR1 - (ld + 1) * AuxTIAA)
r5:
 $pool = {1.0}IAAs
  1 / eps * (AuxTIAA - muIAAs * IAAs)
r6:
 $pool = {1.0}ARFIAA
  pa * ARF * IAAp - pd * ARFIAA
r7:
 $pool = {1.0}ARF2
  qa * pow(ARF, 2) - qd * ARF2
r8:
 $pool = {1.0}Aux
  muAux * (alphaAux - Aux) - 1 / eps * etaAuxTIR1 * (ka * Aux * TIR1 - kd * AuxTIR1)
r9:
 $pool = {1.0}PINm
  phiPINp * (lambda1 * F5a + F5b) - PINm


#parameters
alphaAHK = 1.0
alphaARF = 1.0
alphaARRB = 2.0
alphaAux = 1.0
alphaCk = 1.0
alphaPH = 1.0
alphaTIR1 = 1.0
deltaARRAp = 1.0
deltaARp = 1.0
deltaCRp = 1.0
deltaIAAp = 1.0
deltaPINp = 1.0
eps = 0.01
etaAHKph = 1.0
etaARFIAA = 1.0
etaAuxTIR1 = 10.0
etaCkPh = 1.0
ka = 100.0
kd = 1.0
la = 0.5
lambda1 = 0.1
lambda3 = 0.02
ld = 0.1
muAux = 0.1
muCk = 0.1
muIAAs = 1.0
pa = 10.0
pd = 10.0
phiARRAp = 100.0
phiARp = 2.0
phiCRp = 2.0
phiIAAp = 100.0
phiPINp = 100.0
psiARF = 0.1
psiARFIAA = 0.1
qa = 1.0
qd = 1.0
ra = 1.0
rd = 1.0
sa = 1.0
sd = 1.0
thARFIAA = 0.1
thARRAph = 0.1
thARRBph = 0.1
thetaARF = 0.1
thetaARF2 = 0.01
thetaARp = 0.1
thetaIAAp = 0.1
ua = 1.0
ud = 1.0
cell = 1.0


#initial values
AHKph = 1.0
ARF2 = 0.0
ARFIAA = 0.0
ARRAm = 0.0
ARRAp = 0.0
ARRAph = 0.0
ARRBph = 0.0
ARm = 0.0
ARp = 0.0
Aux = 1.0
AuxTIAA = 0.0
AuxTIR1 = 0.0
CRm = 0.0
CRp = 0.0
Ck = 1.0
IAAm = 0.0
IAAp = 0.0
IAAs = 0.0
PINm = 0.0
PINp = 0.0


#assignment rules
!F CkAHK = alphaAHK - etaAHKph * (AHKph + CkAHKph)
!F CkAHKph = alphaPH - AHKph - ARRAph - ARRBph
!F ARRBp = alphaARRB - etaAHKph * ARRBph
!F ARF = alphaARF - 2.0 * ARF2 - ARFIAA
!F TIR1 = alphaTIR1 - AuxTIR1 - AuxTIAA
!F F6 = ARp / thetaARp / (1.0 + ARp / thetaARp)
!F F5b = (ARF2 / thetaARF2 + pow(ARF, 2.0) / psiARF) / (1.0 + ARF / thetaARF + ARF2 / thetaARF2 + ARFIAA / thARFIAA + ARF * IAAp / psiARFIAA + pow(ARF, 2.0) / psiARF)
!F F5a = ARF / thetaARF / (1.0 + ARF / thetaARF + ARF2 / thetaARF2 + ARFIAA / thARFIAA + ARF * IAAp / psiARFIAA + pow(ARF, 2.0) / psiARF)
!F F4 = ARRBph / thARRBph / (1.0 + ARRAph / thARRAph + ARRBph / thARRBph)
!F F3 = ARRBph / thARRBph / (1.0 + ARF / thetaARF + ARF2 / thetaARF2 + ARFIAA / thARFIAA + ARF * IAAp / psiARFIAA + pow(ARF, 2.0) / psiARF + ARRBph / thARRBph)
!F F2 = (ARF2 / thetaARF2 + pow(ARF, 2.0) / psiARF) / (1.0 + ARF / thetaARF + ARF2 / thetaARF2 + ARFIAA / thARFIAA + ARF * IAAp / psiARFIAA + pow(ARF, 2.0) / psiARF + ARRBph / thARRBph)
!F F1 = ARF / thetaARF / (1.0 + ARF / thetaARF + ARF2 / thetaARF2 + ARFIAA / thARFIAA + ARF * IAAp / psiARFIAA + pow(ARF, 2.0) / psiARF + ARRBph / thARRBph)