#bosdriesz2

FIX: r

#reactions
vAAsynt1:
 $pool = {1.0}aa1
  bm * e1 * kn1 * (1 - r / rmax) / (nAmet * (1 + aa1 / kIa1))
vAAsynt10:
 $pool = {1.0}aa10
  bm * e10 * kn10 * (1 - r / rmax) / (nAmet * (1 + aa10 / kIa10))
vAAsynt11:
 $pool = {1.0}aa11
  bm * e11 * kn11 * (1 - r / rmax) / (nAmet * (1 + aa11 / kIa11))
vAAsynt12:
 $pool = {1.0}aa12
  bm * e12 * kn12 * (1 - r / rmax) / (nAmet * (1 + aa12 / kIa12))
vAAsynt13:
 $pool = {1.0}aa13
  bm * e13 * kn13 * (1 - r / rmax) / (nAmet * (1 + aa13 / kIa13))
vAAsynt14:
 $pool = {1.0}aa14
  bm * e14 * kn14 * (1 - r / rmax) / (nAmet * (1 + aa14 / kIa14))
vAAsynt15:
 $pool = {1.0}aa15
  bm * e15 * kn15 * (1 - r / rmax) / (nAmet * (1 + aa15 / kIa15))
vAAsynt16:
 $pool = {1.0}aa16
  bm * e16 * kn16 * (1 - r / rmax) / (nAmet * (1 + aa16 / kIa16))
vAAsynt17:
 $pool = {1.0}aa17
  bm * e17 * kn17 * (1 - r / rmax) / (nAmet * (1 + aa17 / kIa17))
vAAsynt18:
 $pool = {1.0}aa18
  bm * e18 * kn18 * (1 - r / rmax) / (nAmet * (1 + aa18 / kIa18))
vAAsynt19:
 $pool = {1.0}aa19
  bm * e19 * kn19 * (1 - r / rmax) / (nAmet * (1 + aa19 / kIa19))
vAAsynt2:
 $pool = {1.0}aa2
  bm * e2 * kn2 * (1 - r / rmax) / (nAmet * (1 + aa2 / kIa2))
vAAsynt20:
 $pool = {1.0}aa20
  bm * e20 * kn20 * (1 - r / rmax) / (nAmet * (1 + aa20 / kIa20))
vAAsynt3:
 $pool = {1.0}aa3
  bm * e3 * kn3 * (1 - r / rmax) / (nAmet * (1 + aa3 / kIa3))
vAAsynt4:
 $pool = {1.0}aa4
  bm * e4 * kn4 * (1 - r / rmax) / (nAmet * (1 + aa4 / kIa4))
vAAsynt5:
 $pool = {1.0}aa5
  bm * e5 * kn5 * (1 - r / rmax) / (nAmet * (1 + aa5 / kIa5))
vAAsynt6:
 $pool = {1.0}aa6
  bm * e6 * kn6 * (1 - r / rmax) / (nAmet * (1 + aa6 / kIa6))
vAAsynt7:
 $pool = {1.0}aa7
  bm * e7 * kn7 * (1 - r / rmax) / (nAmet * (1 + aa7 / kIa7))
vAAsynt8:
 $pool = {1.0}aa8
  bm * e8 * kn8 * (1 - r / rmax) / (nAmet * (1 + aa8 / kIa8))
vAAsynt9:
 $pool = {1.0}aa9
  bm * e9 * kn9 * (1 - r / rmax) / (nAmet * (1 + aa9 / kIa9))
vbiomass:
 {0.05}taa9 + {0.05}taa8 + {0.05}taa7 + {0.05}taa6 + {0.05}taa5 + {0.05}taa4 + {0.05}taa3 + {0.05}taa20 + {0.05}taa2 + {0.05}taa19 + {0.05}taa18 + {0.05}taa17 + {0.05}taa16 + {0.05}taa15 + {0.05}taa14 + {0.05}taa13 + {0.05}taa12 + {0.05}taa11 + {0.05}taa10 + {0.05}taa1 = $pool
  krib * r / (1 + 0.05 * kapparta * (kappart + tau1 * r - taa1) / (kappart * taa1) + 0.05 * kapparta * (kappart + tau10 * r - taa10) / (kappart * taa10) + 0.05 * kapparta / taa11 + 0.05 * kapparta * (tau11 * r - taa11) / (kappart * taa11) + 0.05 * kapparta / taa12 + 0.05 * kapparta * (tau12 * r - taa12) / (kappart * taa12) + 0.05 * kapparta / taa13 + 0.05 * kapparta * (tau13 * r - taa13) / (kappart * taa13) + 0.05 * kapparta / taa14 + 0.05 * kapparta * (tau14 * r - taa14) / (kappart * taa14) + 0.05 * kapparta / taa15 + 0.05 * kapparta * (tau15 * r - taa15) / (kappart * taa15) + 0.05 * kapparta / taa16 + 0.05 * kapparta * (tau16 * r - taa16) / (kappart * taa16) + 0.05 * kapparta / taa17 + 0.05 * kapparta * (tau17 * r - taa17) / (kappart * taa17) + 0.05 * kapparta / taa18 + 0.05 * kapparta * (tau18 * r - taa18) / (kappart * taa18) + 0.05 * kapparta / taa19 + 0.05 * kapparta * (tau19 * r - taa19) / (kappart * taa19) + 0.05 * kapparta / taa2 + 0.05 * kapparta * (tau2 * r - taa2) / (kappart * taa2) + 0.05 * kapparta / taa20 + 0.05 * kapparta * (tau20 * r - taa20) / (kappart * taa20) + 0.05 * kapparta / taa3 + 0.05 * kapparta * (tau3 * r - taa3) / (kappart * taa3) + 0.05 * kapparta / taa4 + 0.05 * kapparta * (tau4 * r - taa4) / (kappart * taa4) + 0.05 * kapparta / taa5 + 0.05 * kapparta * (tau5 * r - taa5) / (kappart * taa5) + 0.05 * kapparta / taa6 + 0.05 * kapparta * (tau6 * r - taa6) / (kappart * taa6) + 0.05 * kapparta / taa7 + 0.05 * kapparta * (tau7 * r - taa7) / (kappart * taa7) + 0.05 * kapparta / taa8 + 0.05 * kapparta * (tau8 * r - taa8) / (kappart * taa8) + 0.05 * kapparta / taa9 + 0.05 * kapparta * (tau9 * r - taa9) / (kappart * taa9))
vtRNAchar1:
 {1.0}aa1 = {1.0}taa1
  ks1 * sTot1 * aa1 * (tau1 * r - taa1) / (kMaa1 * kMtf1 * (aa1 / kMaa1 + (tau1 * r - taa1) / kMtf1 + aa1 * (tau1 * r - taa1) / (kMaa1 * kMtf1)))
vtRNAchar10:
 {1.0}aa10 = {1.0}taa10
  ks10 * sTot10 * aa10 * (tau10 * r - taa10) / (kMaa10 * kMtf10 * (aa10 / kMaa10 + (tau10 * r - taa10) / kMtf10 + aa10 * (tau10 * r - taa10) / (kMaa10 * kMtf10)))
vtRNAchar11:
 {1.0}aa11 = {1.0}taa11
  ks11 * sTot11 * aa11 * (tau11 * r - taa11) / (kMaa11 * kMtf11 * (aa11 / kMaa11 + (tau11 * r - taa11) / kMtf11 + aa11 * (tau11 * r - taa11) / (kMaa11 * kMtf11)))
vtRNAchar12:
 {1.0}aa12 = {1.0}taa12
  ks12 * sTot12 * aa12 * (tau12 * r - taa12) / (kMaa12 * kMtf12 * (aa12 / kMaa12 + (tau12 * r - taa12) / kMtf12 + aa12 * (tau12 * r - taa12) / (kMaa12 * kMtf12)))
vtRNAchar13:
 {1.0}aa13 = {1.0}taa13
  ks13 * sTot13 * aa13 * (tau13 * r - taa13) / (kMaa13 * kMtf13 * (aa13 / kMaa13 + (tau13 * r - taa13) / kMtf13 + aa13 * (tau13 * r - taa13) / (kMaa13 * kMtf13)))
vtRNAchar14:
 {1.0}aa14 = {1.0}taa14
  ks14 * sTot14 * aa14 * (tau14 * r - taa14) / (kMaa14 * kMtf14 * (aa14 / kMaa14 + (tau14 * r - taa14) / kMtf14 + aa14 * (tau14 * r - taa14) / (kMaa14 * kMtf14)))
vtRNAchar15:
 {1.0}aa15 = {1.0}taa15
  ks15 * sTot15 * aa15 * (tau15 * r - taa15) / (kMaa15 * kMtf15 * (aa15 / kMaa15 + (tau15 * r - taa15) / kMtf15 + aa15 * (tau15 * r - taa15) / (kMaa15 * kMtf15)))
vtRNAchar16:
 {1.0}aa16 = {1.0}taa16
  ks16 * sTot16 * aa16 * (tau16 * r - taa16) / (kMaa16 * kMtf16 * (aa16 / kMaa16 + (tau16 * r - taa16) / kMtf16 + aa16 * (tau16 * r - taa16) / (kMaa16 * kMtf16)))
vtRNAchar17:
 {1.0}aa17 = {1.0}taa17
  ks17 * sTot17 * aa17 * (tau17 * r - taa17) / (kMaa17 * kMtf17 * (aa17 / kMaa17 + (tau17 * r - taa17) / kMtf17 + aa17 * (tau17 * r - taa17) / (kMaa17 * kMtf17)))
vtRNAchar18:
 {1.0}aa18 = {1.0}taa18
  ks18 * sTot18 * aa18 * (tau18 * r - taa18) / (kMaa18 * kMtf18 * (aa18 / kMaa18 + (tau18 * r - taa18) / kMtf18 + aa18 * (tau18 * r - taa18) / (kMaa18 * kMtf18)))
vtRNAchar19:
 {1.0}aa19 = {1.0}taa19
  ks19 * sTot19 * aa19 * (tau19 * r - taa19) / (kMaa19 * kMtf19 * (aa19 / kMaa19 + (tau19 * r - taa19) / kMtf19 + aa19 * (tau19 * r - taa19) / (kMaa19 * kMtf19)))
vtRNAchar2:
 {1.0}aa2 = {1.0}taa2
  ks2 * sTot2 * aa2 * (tau2 * r - taa2) / (kMaa2 * kMtf2 * (aa2 / kMaa2 + (tau2 * r - taa2) / kMtf2 + aa2 * (tau2 * r - taa2) / (kMaa2 * kMtf2)))
vtRNAchar20:
 {1.0}aa20 = {1.0}taa20
  ks20 * sTot20 * aa20 * (tau20 * r - taa20) / (kMaa20 * kMtf20 * (aa20 / kMaa20 + (tau20 * r - taa20) / kMtf20 + aa20 * (tau20 * r - taa20) / (kMaa20 * kMtf20)))
vtRNAchar3:
 {1.0}aa3 = {1.0}taa3
  ks3 * sTot3 * aa3 * (tau3 * r - taa3) / (kMaa3 * kMtf3 * (aa3 / kMaa3 + (tau3 * r - taa3) / kMtf3 + aa3 * (tau3 * r - taa3) / (kMaa3 * kMtf3)))
vtRNAchar4:
 {1.0}aa4 = {1.0}taa4
  ks4 * sTot4 * aa4 * (tau4 * r - taa4) / (kMaa4 * kMtf4 * (aa4 / kMaa4 + (tau4 * r - taa4) / kMtf4 + aa4 * (tau4 * r - taa4) / (kMaa4 * kMtf4)))
vtRNAchar5:
 {1.0}aa5 = {1.0}taa5
  ks5 * sTot5 * aa5 * (tau5 * r - taa5) / (kMaa5 * kMtf5 * (aa5 / kMaa5 + (tau5 * r - taa5) / kMtf5 + aa5 * (tau5 * r - taa5) / (kMaa5 * kMtf5)))
vtRNAchar6:
 {1.0}aa6 = {1.0}taa6
  ks6 * sTot6 * aa6 * (tau6 * r - taa6) / (kMaa6 * kMtf6 * (aa6 / kMaa6 + (tau6 * r - taa6) / kMtf6 + aa6 * (tau6 * r - taa6) / (kMaa6 * kMtf6)))
vtRNAchar7:
 {1.0}aa7 = {1.0}taa7
  ks7 * sTot7 * aa7 * (tau7 * r - taa7) / (kMaa7 * kMtf7 * (aa7 / kMaa7 + (tau7 * r - taa7) / kMtf7 + aa7 * (tau7 * r - taa7) / (kMaa7 * kMtf7)))
vtRNAchar8:
 {1.0}aa8 = {1.0}taa8
  ks8 * sTot8 * aa8 * (tau8 * r - taa8) / (kMaa8 * kMtf8 * (aa8 / kMaa8 + (tau8 * r - taa8) / kMtf8 + aa8 * (tau8 * r - taa8) / (kMaa8 * kMtf8)))
vtRNAchar9:
 {1.0}aa9 = {1.0}taa9
  ks9 * sTot9 * aa9 * (tau9 * r - taa9) / (kMaa9 * kMtf9 * (aa9 / kMaa9 + (tau9 * r - taa9) / kMtf9 + aa9 * (tau9 * r - taa9) / (kMaa9 * kMtf9)))


#parameters
EXTERNAL = 1.0
bm = 1000000.0
e1 = 0.05
e10 = 0.05
e11 = 0.05
e12 = 0.05
e13 = 0.05
e14 = 0.05
e15 = 0.05
e16 = 0.05
e17 = 0.05
e18 = 0.05
e19 = 0.05
e2 = 0.05
e20 = 0.05
e3 = 0.05
e4 = 0.05
e5 = 0.05
e6 = 0.05
e7 = 0.05
e8 = 0.05
e9 = 0.05
kIa1 = 100.0
kIa10 = 100.0
kIa11 = 100.0
kIa12 = 100.0
kIa13 = 100.0
kIa14 = 100.0
kIa15 = 100.0
kIa16 = 100.0
kIa17 = 100.0
kIa18 = 100.0
kIa19 = 100.0
kIa2 = 100.0
kIa20 = 100.0
kIa3 = 100.0
kIa4 = 100.0
kIa5 = 100.0
kIa6 = 100.0
kIa7 = 100.0
kIa8 = 100.0
kIa9 = 100.0
kMaa1 = 100.0
kMaa10 = 100.0
kMaa11 = 100.0
kMaa12 = 100.0
kMaa13 = 100.0
kMaa14 = 100.0
kMaa15 = 100.0
kMaa16 = 100.0
kMaa17 = 100.0
kMaa18 = 100.0
kMaa19 = 100.0
kMaa2 = 100.0
kMaa20 = 100.0
kMaa3 = 100.0
kMaa4 = 100.0
kMaa5 = 100.0
kMaa6 = 100.0
kMaa7 = 100.0
kMaa8 = 100.0
kMaa9 = 100.0
kMtf1 = 1.0
kMtf10 = 1.0
kMtf11 = 1.0
kMtf12 = 1.0
kMtf13 = 1.0
kMtf14 = 1.0
kMtf15 = 1.0
kMtf16 = 1.0
kMtf17 = 1.0
kMtf18 = 1.0
kMtf19 = 1.0
kMtf2 = 1.0
kMtf20 = 1.0
kMtf3 = 1.0
kMtf4 = 1.0
kMtf5 = 1.0
kMtf6 = 1.0
kMtf7 = 1.0
kMtf8 = 1.0
kMtf9 = 1.0
kappart = 500.0
kapparta = 1.0
kn1 = 0.176694097633303
kn10 = 0.175155965973061
kn11 = 0.21893081294363
kn12 = 0.189250477436635
kn13 = 0.182501157889419
kn14 = 0.168797176256554
kn15 = 0.178172813430161
kn16 = 0.163557478413
kn17 = 0.240562151147302
kn18 = 0.234708428287117
kn19 = 0.238857324820486
kn2 = 0.170472469931078
kn20 = 0.192552580536003
kn3 = 0.215014924490181
kn4 = 0.160529780205187
kn5 = 0.154007815602573
kn6 = 0.232252087628394
kn7 = 0.188971693282943
kn8 = 0.202409099613367
kn9 = 0.221447303290898
krib = 20.0
ks1 = 100.0
ks10 = 100.0
ks11 = 100.0
ks12 = 100.0
ks13 = 100.0
ks14 = 100.0
ks15 = 100.0
ks16 = 100.0
ks17 = 100.0
ks18 = 100.0
ks19 = 100.0
ks2 = 100.0
ks20 = 100.0
ks3 = 100.0
ks4 = 100.0
ks5 = 100.0
ks6 = 100.0
ks7 = 100.0
ks8 = 100.0
ks9 = 100.0
nAmet = 300.0
rmax = 81.2522598284765
sTot1 = 1.0
sTot10 = 1.0
sTot11 = 1.0
sTot12 = 1.0
sTot13 = 1.0
sTot14 = 1.0
sTot15 = 1.0
sTot16 = 1.0
sTot17 = 1.0
sTot18 = 1.0
sTot19 = 1.0
sTot2 = 1.0
sTot20 = 1.0
sTot3 = 1.0
sTot4 = 1.0
sTot5 = 1.0
sTot6 = 1.0
sTot7 = 1.0
sTot8 = 1.0
sTot9 = 1.0
tau1 = 0.5
tau10 = 0.5
tau11 = 0.5
tau12 = 0.5
tau13 = 0.5
tau14 = 0.5
tau15 = 0.5
tau16 = 0.5
tau17 = 0.5
tau18 = 0.5
tau19 = 0.5
tau2 = 0.5
tau20 = 0.5
tau3 = 0.5
tau4 = 0.5
tau5 = 0.5
tau6 = 0.5
tau7 = 0.5
tau8 = 0.5
tau9 = 0.5
r = 24.3757
default = 1.0


#initial values
aa1 = 1.0
aa10 = 1.0
aa11 = 1.0
aa12 = 1.0
aa13 = 1.0
aa14 = 1.0
aa15 = 1.0
aa16 = 1.0
aa17 = 1.0
aa18 = 1.0
aa19 = 1.0
aa2 = 1.0
aa20 = 1.0
aa3 = 1.0
aa4 = 1.0
aa5 = 1.0
aa6 = 1.0
aa7 = 1.0
aa8 = 1.0
aa9 = 1.0
taa1 = 7.3127
taa10 = 7.3127
taa11 = 7.3127
taa12 = 7.3127
taa13 = 7.3127
taa14 = 7.3127
taa15 = 7.3127
taa16 = 7.3127
taa17 = 7.3127
taa18 = 7.3127
taa19 = 7.3127
taa2 = 7.3127
taa20 = 7.3127
taa3 = 7.3127
taa4 = 7.3127
taa5 = 7.3127
taa6 = 7.3127
taa7 = 7.3127
taa8 = 7.3127
taa9 = 7.3127


#assignment rules