# Simulation file for Lotka-Voltera reaction
# This file is used to test and demonstrate continuous-time command timing

graphics opengl
graphic_iter 5
accuracy 5

dim 3
species R F
max_mol 20000
molperbox 1
rand_seed 0

time_start 0
time_stop 10
time_step 0.001

output_files stdout

cmd b echo stdout "... before: b ...\n"
cmd b molcount
cmd @ 0.5 echo stdout "... at: @ 0.5 ...\n"
cmd @ 0.5 molcount
cmd @ 1 echo stdout "... interval: i 1 2 0.1 ...\n"
cmd i 1 2 0.1 molcount
cmd @ 3 echo stdout "... geometric: x 3 8 0.02 1.2 ...\n"
cmd x 3 8 0.02 1.2 molcount
cmd a echo stdout "... after: a ...\n"
cmd a molcount

cmd e ifno R stop
cmd e ifno F stop

difc R 100
difc F 100
color R 1 0 0
color F 0 1 0
display_size R 2
display_size F 3

low_wall 0 -100 p
high_wall 0 100 p
low_wall 1 -100 p
high_wall 1 100 p
low_wall 2 -10 p
high_wall 2 10 p
mol 1000 R u u u
mol 1000 F u u u

start_reaction
order 1
max_rxn 2
reactant R Rmultiply      # R -> 2R
rate Rmultiply 10
product Rmultiply R + R
reactant F Fdie          	# F -> 0
rate Fdie 10
end_reaction

start_reaction
order 2
max_rxn 1
reactant R + F Feat     # R+F -> 2F
rate Feat 8000
product Feat F + F
end_reaction

end_file