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";

  }

  );