#callaway3 #reactions v1: $pool = {1.0}T lambda v2: {1.0}T = $pool d*T v3: {1.0}T = {1.0}Tstar k*(1 - kappa)*T*V v4: {1.0}Tstar = $pool deltaACC*Tstar**(1 + omega) v5: $pool = {1.0}V p*Tstar v6: {1.0}V = $pool c*V #parameters NT = 100.0 c = 13.0 d = 0.01 delta = 0.7 deltaACC = 0.7863 k = 8e-07 kappa = 0.788 lambda = 10000.0 omega = 0.01 p = 70.0 default = 1.0 #initial values T = 200000.0 Tstar = 9285.88875605 V = 50000.0 #assignment rules