#liu1

FIX: ACC CK_ex CTR1T IAA RaT ReT

#reactions
reaction_1:
 $pool = {1.0}Auxin
  compartment_1 * function_2(k1a, X, k1)
reaction_10:
 $pool = {1.0}Re_star
  compartment_1 * function_10(k10, PLSp, k10a, Re)
reaction_11:
 {1.0}Re_star = $pool
  compartment_1 * function_12(Re_star, ET, k11)
reaction_12:
 $pool = {1.0}ET
  compartment_1 * function_13(Auxin, CK, k12, k12a)
reaction_13:
 {1.0}ET = $pool
  compartment_1 * function_15(ET, k13)
reaction_14:
 $pool = {1.0}CTR1_star
  compartment_1 * function_16(Re_star, k14, CTR1)
reaction_15:
 {1.0}CTR1_star = $pool
  compartment_1 * function_17(CTR1_star, k15)
reaction_16:
 $pool = {1.0}X
  compartment_1 * function_18(CTR1_star, k16, k16a)
reaction_17:
 {1.0}X = $pool
  compartment_1 * function_19(X, k17)
reaction_18:
 $pool = {1.0}CK
  compartment_1 * function_20(Auxin, k18a, k18)
reaction_19:
 {1.0}CK = $pool
  compartment_1 * function_21(CK, k19)
reaction_2:
 $pool = {1.0}Auxin
  compartment_1 * function_3(k2, k2a, ET, CK, k2b, PLSp, k2c)
reaction_3:
 {1.0}Auxin = $pool
  compartment_1 * function_1(k3, k3a, X, Auxin)
reaction_4:
 $pool = {1.0}Ra_star
  compartment_1 * function_4(k4, Auxin, Ra)
reaction_5:
 {1.0}Ra_star = $pool
  compartment_1 * k5 * Ra_star
reaction_6:
 $pool = {1.0}PLSm
  compartment_1 * function_6(k6, Ra_star, ET, k6a)
reaction_7:
 {1.0}PLSm = $pool
  compartment_1 * function_7(k7, PLSm)
reaction_8:
 $pool = {1.0}PLSp
  compartment_1 * function_8(k8, PLSm)
reaction_9:
 {1.0}PLSp = $pool
  compartment_1 * function_9(k9, PLSp)
v_Auxin:
 {1.0}IAA = {1.0}Auxin
  compartment_1 * k_auxin * IAA
v_Cytokinin:
 {1.0}CK_ex = {1.0}CK
  compartment_1 * k_cytokinin * CK_ex
v_Ethylene:
 {1.0}ACC = {1.0}ET
  compartment_1 * k_ethylene * ACC


#parameters
k1 = 1.0
k10 = 0.0003
k10a = 0.5
k11 = 5.0
k12 = 0.1
k12a = 0.1
k13 = 1.0
k14 = 3.0
k15 = 0.085
k16 = 0.3
k16a = 1.0
k17 = 0.1
k18 = 0.1
k18a = 1.0
k19 = 1.0
k1a = 1.0
k2 = 0.2
k2a = 2.8
k2b = 1.0
k2c = 0.01
k3 = 2.0
k3a = 0.45
k4 = 1.0
k5 = 1.0
k6 = 0.3
k6a = 0.2
k7 = 1.0
k8 = 1.0
k9 = 1.0
k_auxin = 70.0
k_cytokinin = 10.0
k_ethylene = 0.5
ACC = 0.0
CK_ex = 0.0
CTR1T = 0.3
IAA = 0.0
RaT = 1.0
ReT = 0.3
compartment_1 = 1.0
extra = 1.0


#initial values
Auxin = 0.1
CK = 0.1
CTR1_star = 0.3
ET = 0.1
PLSm = 0.1
PLSp = 0.1
Ra_star = 1.0
Re_star = 0.3
X = 0.1


#assignment rules
!F CTR1 = CTR1T - CTR1_star
!F Re = ReT - Re_star
!F Ra = RaT - Ra_star