#rong3



#reactions
v1:
 $pool = {1.0}T
  lambda
v10:
 {1.0}T1 = {1.0}L
  alphaL*k*T*Vi
v11:
 {1.0}L = {1.0}T1
  alphaL*epsilonRT*k*T*Vi
v12:
 $pool = {1.0}V + {1.0}Vi
  delta*NN*T1
v13:
 {1.0}Vi = {1.0}Vni
  delta*epsilonPI*NN*T1
v2:
 {1.0}T = $pool
  dT*T
v3:
 {1.0}T = {1.0}T1
  k*T*Vi
v4:
 {1.0}L = $pool
  dL*L
v5:
 {1.0}L = {1.0}T1
  a*L
v6:
 {1.0}T1 = $pool
  delta*T1
v7:
 {1.0}V + {1.0}Vi = $pool
  c*Vi
v8:
 {1.0}V + {1.0}Vni = $pool
  c*Vni
v9:
 {1.0}T1 = {1.0}T
  epsilonRT*k*T*Vi


#parameters
NN = 4000.0
a = 0.1
alphaL = 0.001
c = 23.0
dL = 0.004
dT = 0.01
delta = 1.0
epsilonPI = 0.65
epsilonRT = 0.0
k = 2.4e-08
lambda = 10000.0
default = 1.0


#initial values
L = 2.0
T = 600000.0
T1 = 0.3
V = 50000.0
Vi = 50000.0
Vni = 0.0


#assignment rules
!F ratio = Log(Vni/Vi)/Log(10.0)