#hong1



#reactions
r1:
 $pool = {1.0}M
  system * function_4_r1
r10:
 $pool = {1.0}TF + {1.0}CP2
  system * kicd * IC
r11:
 {1.0}CP2 = $pool
  system * kcp2d * CP2
r12:
 {1.0}TF + {1.0}CP2 = $pool
  system * kica * CP2 * TF
r13:
 {1.0}CP2 = $pool
  system * function_4_r13
r14:
 {1.0}CP2 = $pool
  system * chk2 * CP2
r15:
 $pool = {1.0}TF
  system * kcp2d * IC
r16:
 $pool = {1.0}TF
  system * chk2c * IC
r17:
 $pool = {1.0}TF
  system * function_4_r17
r2:
 $pool = {1.0}M
  system * function_4_r2
r3:
 {1.0}M = $pool
  system * kmd * M
r4:
 $pool = {1.0}CP
  system * function_1
r5:
 {1.0}CP = $pool
  system * kcpd * CP
r6:
 {2.0}CP = {1.0}CP2
  system * ka * pow(CP, 2)
r7:
 {1.0}CP2 = {2.0}CP
  system * kd * CP2
r8:
 {1.0}CP = $pool
  system * function_4_r8
r9:
 {1.0}CP = $pool
  system * chk2 * CP


#parameters
Dex = 0.0
ICtot = 1.0
IR_pulse_duration = 2.0
J = 0.3
Jp = 0.05
TFtot = 0.5
chk2ToggleOnOff = 1.0
chk2_value = 0.25
chk2c_value = 0.1
chk2pulseactive = 1.0
chk2pulseactive_elim = 1.0
chk2pulseactive_elim_elim = 0.0
initialValue_CP = 0.039
initialValue_CP2 = 0.088
initialValue_M = 1.35
initialValue_TF = 0.06
ka = 100.0
kcp2d = 0.0525
kcpd = 0.525
kcps = 0.5
kd = 0.01
kica = 20.0
kicd = 0.01
kmd = 0.1
kms = 1.0
kp1 = 10.0
kp2 = 0.1
n = 2.0
trigONtime = 4.0
system = 1.0


#initial values
CP = 0.039
CP2 = 0.088
M = 1.35
TF = 0.06


#assignment rules
!F IC = TFtot - TF
!F function_1 = kcps * M
!F function_4_r17 = kp2 * IC / (Jp + CP + 2.0 * CP2 + 2.0 * IC) / system
!F function_4_r2 = kms * pow(TF, n) / (pow(J, n) + pow(TF, n)) / system
!F chk2pulseonoff = chk2pulseactive - (chk2pulseactive_elim_elim + chk2pulseactive_elim)
!F function_4_r8 = kp1 * CP / (Jp + CP + 2.0 * CP2 + 2.0 * IC) / system
!F function_4_r13 = kp2 * CP2 / (Jp + CP + 2.0 * CP2 + 2.0 * IC) / system
!F function_4_r1 = Dex / system
!F CPtot = CP + 2.0 * CP2 + 2.0 * IC
!F trigOFFtime = trigONtime + IR_pulse_duration
!F chk2c = chk2ToggleOnOff*(chk2pulseonoff * (chk2c_value / IR_pulse_duration))
!F chk2 = chk2ToggleOnOff*(chk2pulseonoff * (chk2_value / IR_pulse_duration))