diff --git a/input/kinetic_example.cfg b/input/kinetic_example.cfg
new file mode 100644
index 0000000000000000000000000000000000000000..5109262450b59728429fded9e074a7ead7c57ef8
--- /dev/null
+++ b/input/kinetic_example.cfg
@@ -0,0 +1,266 @@
+filename="kinetic_example.cfg";
+
+#########################
+## Species definitions ##
+#########################
+    
+species = (
+  {
+    name="o2";
+    mass=0.0;
+    charge=0.0;
+  },
+  {
+    name="n2";
+    mass=0.0;
+    charge=0.0;
+  },
+  {
+    ## This is the generic neutral species used for compound reactions that take ##
+    ## any neutral as input. Make sure to set its density as the sum of the      ##
+    ## densities of all actual neutral species in the seed section               ##
+    name="neutral";
+    mass=0.0;
+    charge=0.0;
+  },
+  {
+    name="electrons";
+    mass=1.0;
+    charge=-1.0;
+  },
+  {
+    ## We don't distinguish between different positive ions, so one        ##
+    ## generic ion, X+, is used                                            ##
+    name="xplus";
+    mass=0.0;
+    charge=1.0;
+  },
+  {
+    name="o2minus";
+    mass=0.0;
+    charge=-1.0;
+  },
+  {
+    name="ominus";
+    mass=0.0;
+    charge=-1.0;
+  },
+    ## These dummy species are used for the photo-ionization code.          ##
+    ## D+ represents N2+, D- is just to preserve the charge-balance.        ##
+  {
+    name="dummyminus";
+    mass=0.0;
+    charge=-1.0;
+  },
+  {
+    name="dummyplus";
+    mass=0.0;
+    charge=1.0;
+  }
+  );
+    
+######################
+## Seed definitions ##
+######################
+    
+seed = (
+  {
+    species="electrons";
+    type="gaussian";
+    value=1.0e6;
+    x0=0.;
+    y0=0.;
+    z0=13044.;
+    sigma_x=32.0;
+    sigma_y=32.0;
+    sigma_z=32.0;
+  },
+  {
+    species="xplus";
+    type="gaussian";
+    value=1.0e6;
+    x0=0.;
+    y0=0.;
+    z0=13044.;
+    sigma_x=32.0;
+    sigma_y=32.0;
+    sigma_z=32.0;
+  },
+  {
+    species="xplus";
+    type="constant";
+    value=0.0e-8;
+  },
+  {
+    species="o2minus";
+    type="constant";
+    value=0.0e-8;
+  },
+  {
+    species="n2";
+    type="constant";
+    value=4.16e4;
+  },
+  {
+    species="o2";
+    type="constant";
+    value=1.04e4;
+  },
+  {
+    species="dummyminus";
+    type="constant";
+    value=0.0;
+  },
+  {
+    species="dummyplus";
+    type="constant";
+    value=0.0;
+  },
+  {
+    ## The value for this seed should be the sum of the values for all      ##
+    ## normal neutral seeds.                                                ##
+    species="neutral";
+    type="constant";
+    value=5.2e4;
+  }
+  );
+    ##########################
+    ## Reaction definitions ##
+    ##########################
+    
+    ## Impact ionization of N2: e + N2 -> 2 e + X+ + D+ + D-                ##
+    ## Since the rate of N2-impact-ionization is a source-term for the      ##
+    ## photoionization-equation, we use the dummy species to track this     ##
+    ## reaction                                                             ##
+
+reactions = (
+  {
+    nin=2;
+    nout=5;
+    specin0="n2";
+    specin1="electrons";
+    specout0="electrons";
+    specout1="electrons";
+    specout2="xplus";
+    specout3="dummyminus";
+    specout4="dummyplus";
+    reacttable="input/n2ion.reaction";
+  },
+  {
+    ## Impact ionization of O2:         e + O2 -> 2 e + X+                  ##
+    nin=2;
+    nout=3;
+    specin0="o2";
+    specin1="electrons";
+    specout0="electrons";
+    specout1="electrons";
+    specout2="xplus";
+    reacttable="input/o2ion.reaction";
+  },
+  {
+    ## Attachment to O2:                e + O2 -> O- + neutral              ##
+    nin=2;
+    nout=2;
+    specin0="o2";
+    specin1="electrons";
+    specout0="o2";
+    specout1="neutral";
+    reacttable="input/o2attach2.reaction";
+  },
+  {
+    ## Attachment to O2:                e + O2 + O2 -> O2- + O2             ##
+    nin=3;
+    nout=2;
+    specin0="o2";
+    specin1="o2";
+    specin2="electrons";
+    specout0="o2";
+    specout1="o2minus";
+    reacttable="input/o2attach3.reaction";
+  },
+  {
+    ## Attachment to O2:                e + O2 + N2 -> O2- + N2             ##
+    nin=3;
+    nout=2;
+    specin0="o2";
+    specin1="n2";
+    specin2="electrons";
+    specout0="o2minus";
+    specout1="n2";
+    reacttable="input/n2o2attach.reaction";
+  },
+  {
+    ## Detachment from O2-:             O2- + O2 -> e + 2 O2                ##
+    nin=2;
+    nout=3;
+    specin0="o2";
+    specin1="o2minus";
+    specout0="o2";
+    specout1="o2";
+    specout2="electrons";
+    reacttable="input/o2detach.reaction";
+  },
+  {
+    ## Detachment from O2-:             O2- + N2 -> e + O2 + N2             ##
+    nin=2;
+    nout=3;
+    specin0="n2";
+    specin1="o2minus";
+    specout0="n2";
+    specout1="electrons";
+    specout2="o2";
+    reacttable="input/o2n2detach.reaction";
+  },
+  {
+    ## Electron/ion recombination:      e + X+ -> X                         ##
+    nin=2;
+    nout=1;
+    specin0="xplus";
+    specin1="electrons";
+    specout0="neutral";
+    reacttable="input/elecionrecomb.reaction";
+  },
+  {
+    ## Ion/ion recombination:           O- + X+ -> X                        ##
+    nin=2;
+    nout=1;
+    specin0="xplus";
+    specin1="ominus";
+    specout0="neutral";
+    reacttable="input/ominusrecomb.reaction";
+  },
+  {
+    ## Ion/ion recombination:           O- + X+ + X -> 2 X + O              ##
+    nin=3;
+    nout=3;
+    specin0="xplus";
+    specin1="ominus";
+    specin2="neutral";
+    specout0="neutral";
+    specout1="neutral";
+    specout2="neutral";
+    reacttable="input/ominusneutralrecomb.reaction";
+  },
+  {
+    ## Ion/ion recombination:           O2- + X+ -> O2 + X                  ##
+    nin=2;
+    nout=2;
+    specin0="xplus";
+    specin1="o2minus";
+    specout0="o2";
+    specout1="neutral";
+    reacttable="input/o2minusrecomb.reaction";
+  },
+  {
+    ## Ion/ion recombination:           O2- + X+ + X -> O2 + X + X          ##
+    nin=3;
+    nout=3;
+    specin0="xplus";
+    specin1="o2minus";
+    specin2="xplus";
+    specout0="o2";
+    specout1="neutral";
+    specout2="neutral";
+    reacttable="input/o2minusneutralrecomb.reaction";
+  }
+  );