This [script](https://osf.io/ur3bz) prints the creation and destruction integrated currents for a species in zones in a Webnucleo XML file. To use, first download the file by typing **curl -o print_currents.py -J -L https://osf.io/b85k6/download** Use the script on an output Webnucleo XML file (in this example, called *out.xml*). This will be an output from a Webnucleo network calculation, but, to obtain an [example](https://osf.io/zhrfn/), type **curl -o out.xml -J -L https://osf.io/ytud2/download** Then print the integrated currents from *out.xml* for a species (in this example, *si28*) in particular zones (in this example, the first 10 zones) by typing **python print_currents.py out.xml si28 "[position() <= 10]"** Note that negative currents indicate that the reverse direction for the reaction dominates. To restrict the print out to only currents with absolute magnitude greater than some value (in this example, 1.e-10), use the cutoff option: **python print_currents.py out.xml si28 "[position() <= 10]" --cutoff 1.e-10**