dt = 1e-6
TIME_STEP = dt
TIME_STEP_MAX = dt
ITERATIONS = 101
EFFECTOR_GRID_DENSITY = 10000

DEFINE_MOLECULES {
  A { DIFFUSION_CONSTANT_3D = 2e-6 }
  B { DIFFUSION_CONSTANT_3D = 2e-6 }
  A.B { DIFFUSION_CONSTANT_3D = 0 }
}


DEFINE_REACTIONS {
  A + B -> A.B [2e8]
}

A_site SPHERICAL_RELEASE_SITE {
  LOCATION = [-0.1,0,0]
  MOLECULE = A
  NUMBER_TO_RELEASE = 1000
  SITE_DIAMETER = 0.05
}

B_site SPHERICAL_RELEASE_SITE {
  LOCATION = [0.1,0,0]
  MOLECULE = B
  NUMBER_TO_RELEASE = 1000
  SITE_DIAMETER = 0.05
}


INSTANTIATE world OBJECT {
  A_site OBJECT A_site {}
  B_site OBJECT B_site {}
}

VIZ_DATA_OUTPUT {
  MODE = DX
  STATE_VALUES {
    A = 1
    B = 2
    A.B = 3
  }
  MOLECULE_FILE_PREFIX = "./viz/unbounded_mmi"
  ITERATION_LIST = [0,1,10,100]
}

REACTION_DATA_OUTPUT {
  STEP = 1*dt
  {COUNT[A,WORLD]} => "./dat/A.dat"
  {COUNT[B,WORLD]} => "./dat/B.dat"
  {COUNT[A.B,WORLD]} => "./dat/A.B.dat"
  {COUNT[A,WORLD]+COUNT[A.B,WORLD]} => "./dat/A_A.B_tot.dat"
}