diff --git a/siesta.pdf b/siesta.pdf
new file mode 100644
index 0000000..d2c3642
Binary files /dev/null and b/siesta.pdf differ
diff --git a/tbtrans.pdf b/tbtrans.pdf
new file mode 100644
index 0000000..114d994
Binary files /dev/null and b/tbtrans.pdf differ
diff --git a/tselecs.sh b/tselecs.sh
new file mode 100755
index 0000000..7d48a5d
--- /dev/null
+++ b/tselecs.sh
@@ -0,0 +1,659 @@
+#!/bin/bash
+
+# Script to create output for multiple electrodes
+# and/or multiple chemical potentials
+
+_this=$(basename $0)
+
+function error {
+    echo "ERROR: something went wrong... Please see the following output:"
+    while [ $# -gt 0 ]; do
+	echo "$1"
+	shift
+    done
+    exit 1
+}
+
+
+# We must require to be able to use the hash-able arrays
+# This allows us to reuse the same input twice more easily
+help=0
+def=2
+tbt=0
+# Options about what to print
+print_mu=0
+print_el=0
+print_c=0
+declare -a opts
+_N=0
+volt=0
+while [ $# -gt 0 ]; do
+    opt=$1
+    case $opt in
+	--*)
+	    opt=${1:1} ;;
+    esac
+    case $opt in
+	-def|-orig|-original)
+	    # Will produce output for standard TS
+	    def=2 ; shift
+	    ;;
+	-V|-bias)
+	    volt=$1 ; shift ; shift
+	    ;;
+	-1|-2|-3|-4|-5|-6|-7|-8|-9)
+	    # Will produce output for standard TS
+	    def=${opt:1} ; shift
+	    ;;
+	-only-el|-only-elec|-only-electrode)
+	    print_el=1 ; shift
+	    ;;
+	-only-mu|-only-chem|-only-chemical)
+	    print_mu=1 ; shift
+	    ;;
+	-only-c|-only-contour)
+	    print_c=1 ; shift
+	    ;;
+	-T|-tbtrans)
+	    # Prefix with TBT instead of TS
+	    tbt=1 ; shift
+	    ;;
+	-help|-h)
+	    # We need to print help
+	    help=1 ; shift
+	    ;;
+	-*)
+	    error "Unknown option: $opt"
+	    ;;
+	*)
+	    break
+    esac
+done
+
+# Correct what to print (if the user)
+# has not specified anything
+if [ $print_el -eq 0 ] && [ $print_mu -eq 0 ] && [ $print_c -eq 0 ]; then
+    print_el=1
+    print_mu=1
+    print_c=1
+fi
+
+if [ $help -eq 1 ]; then
+    fmt="      %20s : %s\n"
+    echo "Usage:"
+    echo "    $_this <options>"
+    echo ""
+    echo "For intrinsic transiesta runs call this:"
+    echo "    $_this -2 -V <bias>"
+    echo ""
+    echo "If one wishes to divide it in several files do:"
+    echo "    $_this -only-el > ELECS.fdf"
+    echo "    $_this -only-mu > CHEMPOTS.fdf"
+    echo "    $_this -only-c  > CONTOURS.fdf"
+    echo ""
+    echo "Only print out tbtrans related options:"
+    echo "    $_this -T|--tbtrans"
+    echo ""
+    echo "For specifying the contour energy levels you can use these options:"
+    printf "$fmt" "-Emin <val>" "band bottom energy in eV (-40. eV)"
+    printf "$fmt" "-dE <val>" "real-axis distance between integration points (0.01 eV)"
+    echo ""
+    echo "There are preset creations of 2,3 and 4 electrodes, call:"
+    echo "    $_this -[2-9]"
+    echo "to create the equivalent systems, note that these options can be accompanied by further setup."
+    echo ""
+    echo "Available options for the chemical potentials:"
+    echo "    -mu<idx>-X where X is one of the following:"
+    printf "$fmt" "mu <val>" "the chemical potential value in eV, or in fractions of V (V/2, |V|/3, etc.)"
+    printf "$fmt" "name <name>" "the name of the chemical potential (to be referenced by the electrodes)"
+    echo ""
+    echo "A shorthand notation for supplying the options is:"
+    echo "    -mu<idx> Y=<val>,Z=<val> where Y,Z is one of X above"
+    echo ""
+    echo ""
+    echo "Available options for the electrodes:"
+    echo "    -el<idx>-X where X is one of the following:"
+    printf "$fmt" "mu <name>" "the chemical potential name as given by -mu-name"
+    printf "$fmt" "name <name>" "the name of the electrode"
+    printf "$fmt" "tshs <file>" "the electrode TSHS-file"
+    printf "$fmt" "semi-inf <val>" "the semi-infinite direction [-+][a-c|a[1-3]]"
+    printf "$fmt" "pos <val>" "the first atom of the electrode in the structure"
+    printf "$fmt" "end-pos <val>" "the last atom of the electrode in the structure"
+    printf "$fmt" "used-atoms|ua <val>" "number of atoms used in the electrode (all)"
+    printf "$fmt" "bloch-a[1-3] <val>" "the Bloch expansion in the equivalent direction (1)"
+    printf "$fmt" "bulk <T|F>" "whether the electrode is a bulk electrode (true)"
+    printf "$fmt" "gf <file>" "name of the Greens function file (TSGF[el<idx>-name])"
+    echo ""
+    echo "A shorthand notation for supplying the options is:"
+    echo "    -el<idx> Y=<val>,Z=<val> where Y,Z is one of X above"
+    echo ""
+    echo "You can combine -2 with additional electrode options to easily generate different schemes"
+    echo "To generate 3 electrodes and 2 chemical potentials you can do this:"
+    echo "    $_this -2 -el3 name=Top,mu=Left,inf-dir=+a2,end-pos=-1"
+    exit 1
+fi
+
+prefix=TS
+# If we are doing tbtrans output
+# we should not print the contour information
+if [ $tbt -eq 1 ]; then
+    prefix=TBT
+    print_c=0
+fi
+
+# Add an option to the option array, enables easy creations
+function add_opt {
+    local rem=0
+    [ "x$1" == "x-r" ] && rem=1 && shift
+    local opt=$1 ; shift
+    local val=$1 ; shift
+    if [ $rem -eq 1 ]; then
+	rem_opt $opt
+    fi
+    # array indices are zero-based
+    local tmp=$(get_opt $opt 1)
+    [ ${#tmp} -gt 0 ] && return 0
+    opts+=($opt)
+    opts+=($val)
+    let _N++
+    let _N++
+}
+
+# Removes an option from the option array, enables customization of the options in-code
+function rem_opt {
+    local opt=
+    local i=
+    while [ $# -gt 0 ]; do
+	opt=$1 ; shift
+	i=0
+	while [ $i -le $_N ]; do
+            local val=${opts[$i]}
+            if [ "x$val" == "x$opt" ]; then
+                #echo "Removing ($val|$opt: ${opts[$i]}"
+                unset opts[$i]
+                let i++
+                #echo "Removing: ${opts[$i]}"
+                unset opts[$i]
+		break
+            fi
+            let i++
+        done
+    done
+    opts=( "${opts[@]}" )
+    _N=${#opts[@]}
+}
+
+
+# Lower-case function
+function lc { printf "%b" "$@" | tr '[A-Z]' '[a-z]' ; }
+
+
+function get_opt {
+    # Returns the option from the array input
+    # get_opt <opt-name (with ONE dash)> [number of arguments gobbled]
+    # If [...] is 0:
+    #   0) will return 0 or 1 whether or not the option exists
+    #   1-9) will return the consecutive [...] number of elements
+    local opt="$1" ; shift
+    local count=0
+    if [ $# -gt 0 ]; then
+	count=$1 ; shift
+    fi
+    local i=0 ; local j=0
+    local ret=""
+    while [ $i -le $_N ]; do
+	local oo=$(lc ${opts[$i]})
+	#echo "searching for: $opt in $oo" 1>&2
+	if [ "x$opt" == "x$oo" ]; then
+	    if [ $count -gt 0 ]; then
+		ret=""
+		for j in `seq 1 $count` ; do
+		    let i++
+		    ret="$ret ${opts[$i]}"
+		done
+	    else
+		ret=1
+	    fi
+	    break
+	fi
+	let i++
+    done
+    [ ${#ret} -eq 0 ] && [ $count -eq 0 ] && ret=0
+    # trim white-space
+    ret=${ret%% }
+    ret=${ret## }
+    printf "%b" "$ret"
+}
+
+# Expands the "short-hand" keywords to their full equivalents
+function expand_key {
+    local nm="$1" ; shift
+    local i=1
+    local reg_opt=
+    local opt=
+    local val=
+    while : ; do
+	reg_opt=$(get_opt -$nm$i 1)
+	# If the option does not exist, we simply return
+	[ ${#reg_opt} -eq 0 ] && [ $i -gt 6 ] && break
+	# Delete, so that we only keep the wanted options (not necessary, but)
+	rem_opt -$nm$i
+	reg_opt="${reg_opt//,/ }"
+	reg_opt="${reg_opt//;/ }"
+        # Loop over comma/semi-colon separated entries
+	for reg in $reg_opt ; do
+            # Add to the option list
+            # Retrive the option and value
+	    opt=${reg%=*}
+	    val=${reg#*=}
+	    add_opt -r "-$nm$i-$opt" "$val"
+	done
+	let i++
+    done
+}
+
+function reduce {
+    local -a tmp
+    local i=0
+    local j=0
+    while [ $j -le $_N ]; do
+	if [ ${#opts[$j]} -gt 0 ]; then
+	    tmp[$i]=${opts[$j]}
+	    let i++
+	fi
+	let j++
+    done
+    opts=( "${tmp[@]}" )
+    _N=$i
+}
+
+for i in `seq 1 $def` ; do
+    # The user has requested a default setting for some electrodes
+    # We add them to the option stack
+    add_opt "-mu$i-name" "mu-$i"
+    add_opt "-mu$i-mu" "<mu-mu$i>"
+    add_opt "-el$i-name" "el-$i"
+    add_opt "-el$i-mu" "mu-$i"
+    add_opt "-el$i-pos" "<input-value>"
+    add_opt "-el$i-semi-inf" "+a3"
+done
+
+# In case of only two electrodes we have
+# the particular case of the traditional transiesta code
+case $def in
+    1)
+	add_opt -r "-mu1-name" "single"
+	add_opt -r "-mu1-mu" "0."
+	add_opt -r "-el1-name" "single"
+	add_opt -r "-el1-mu" "single"
+	rem_opt -el1-pos -el1-end-pos -el1-semi-inf
+	add_opt "-el1-pos" "1"
+	add_opt "-el1-semi-inf" "-a3"
+	;;
+    2)
+	add_opt -r "-mu1-name" "Left"
+	add_opt -r "-mu1-mu" "V/2"
+	add_opt -r "-el1-name" "Left"
+	add_opt -r "-el1-mu" "Left"
+	rem_opt -el1-pos -el1-end-pos -el1-semi-inf
+	add_opt "-el1-pos" "1"
+	add_opt "-el1-semi-inf" "-a3"
+	add_opt -r "-mu2-name" "Right"
+	add_opt -r "-mu2-mu" "-V/2"
+	add_opt -r "-el2-name" "Right"
+	add_opt -r "-el2-mu" "Right"
+	rem_opt -el2-pos -el2-end-pos -el2-semi-inf
+	add_opt "-el2-end-pos" "-1"
+	add_opt "-el2-semi-inf" "+a3"
+	;;
+esac    
+
+
+# Add the user options
+while [ $# -gt 0 ]; do
+    opt=$1 ; shift
+    case $opt in
+	--*)
+	    opt=${1:1}
+	    rem_opt $opt ;;
+	-el*|-mu*)
+	    rem_opt $opt ;;
+    esac
+    let _N++
+    opts[$_N]="$opt"
+done
+
+# Expand again after the user things
+expand_key mu
+expand_key el
+reduce
+
+
+function count {
+    local nm=$1 ; shift
+    local n=0
+    local i=
+    local oldn=-1
+    local opt=
+    local j=
+    while [ $oldn -ne $n ]; do
+	oldn=$n
+	i=$n
+	let i++
+	j=0
+	while [ $j -lt $_N ]; do
+	    opt=${opts[$j]}
+	    case $opt in 
+		-${nm}${i}*)
+		    n=$i ;;
+	    esac
+	    let j++
+	done
+    done
+    printf "%b" "$n"
+}
+
+# Save the options 
+old_opts=( "${opts[@]}" )
+i=0
+while [ $i -lt $_N ]; do
+    case ${opts[$i]} in
+	-el*)
+	    unset opts[$i]
+	    let i++
+	    unset opts[$i]
+	    ;;
+    esac
+    let i++
+done
+# Remove empty array elements
+reduce
+
+# Energy for the contour
+emin=$(get_opt -emin 1)
+[ ${#emin} -eq 0 ] && emin=-40.
+emin="$emin eV"
+de=$(get_opt -de 1)
+[ ${#de} -eq 0 ] && de=0.01
+
+
+
+# The chemical potentials should be denoted like this option:
+# -mu1 <options for the chemical potential in comma separated list>
+# Another way of supplying information regarding the chemical potentials
+#  -mu1-name <name>
+#  -mu1-mu   <the value of the chemical potential>
+_mus=$(count mu)
+# First we need to figure out the number of chemical potentials:
+# I don't suspect ANY user will go above 10 chemical potentials
+[ $_mus -eq 0 ] && \
+    error "You have zero chemical potentials." "Please supply at least one..."
+
+# Create the chemical potential array
+_mu_names=()
+j=0
+for i in `seq 1 $_mus` ; do
+    _mu_names+=( "$(get_opt -mu$i-name 1)" )
+    let j++
+done
+
+function mu_e_correct {
+    local mu=$1 ; shift
+    if [ "x${mu:0:1}" == "x-" ]; then
+	# we have a negative sign
+	if [ "x${mu:1:2}" != "x|V" ] && [ "x${mu:1:1}" != "xV" ]; then
+	    mu="$mu eV"
+	fi
+    else
+	if [ "x${mu:0:2}" != "x|V" ] && [ "x${mu:0:1}" != "xV" ]; then
+	    mu="$mu eV"
+	fi
+    fi
+    printf "%b" "$mu"
+}
+	
+if [ $print_mu -eq 1 ]; then
+
+# Print out chemical potential block:
+echo "$prefix.Voltage $volt eV"
+echo "%block $prefix.ChemPots"
+for mu in ${_mu_names[@]} ; do
+    echo "  $mu"
+done
+echo "%endblock $prefix.ChemPots"
+echo ""
+
+
+function create_mu {
+    local mu=$1 ; shift
+    local name=$(get_opt -mu$mu-name 1)
+    if [ ${#name} -eq 0 ]; then
+	error "Chemical potential: $mu" "Has not received a name (-mu$mu-name)"
+    fi
+    local chem=$(get_opt -mu$mu-mu 1)
+    if [ ${#chem} -eq 0 ]; then
+	error "Chemical potential: $mu" "Has not received a chemical potential value (-mu$mu-mu)"
+    fi
+    chem=$(mu_e_correct $chem)
+    echo "%block $prefix.ChemPot.$name"
+    echo "  mu $chem"
+    if [ $tbt -eq 0 ]; then
+	echo "  contour.eq"
+	echo "    begin"
+	echo "      C-$name"
+	echo "      T-$name"
+	echo "    end"
+    fi
+    echo "%endblock $prefix.ChemPot.$name"
+}
+for i in `seq 1 $_mus` ; do
+    create_mu $i
+done
+fi
+
+# Create all the contours
+
+if [ $print_c -eq 1 ]; then
+
+echo ""
+echo "TS.Contours.Eq.Pole 2.5 eV"
+
+for i in `seq 1 $_mus` ; do
+    mu=$(get_opt -mu$i-name 1)
+    # Create the contours
+    echo "%block TS.Contour.C-$mu"
+    echo "  part circle"
+    e=$(get_opt -mu$i-mu 1)
+    e=$(mu_e_correct $e)
+    [ ${e:0:1} != "-" ] && e="+ $e"
+    echo "   from $emin $e to -10 kT $e"
+    echo "     points 25"
+    echo "      method g-legendre"
+    echo "%endblock TS.Contour.C-$mu"
+
+    # Create the contours
+    echo "%block TS.Contour.T-$mu"
+    echo "  part tail"
+    echo "   from prev to inf"
+    echo "     points 10"
+    echo "      method g-fermi"
+    echo "%endblock TS.Contour.T-$mu"
+
+done
+
+echo ""
+echo "MSG: You may need to correct the TS.Contours.nEq for signs of bias" >&2
+echo "MSG: The energies must be in increasing order (you may use |V|/2 to designate absolute values)" >&2
+
+if [ $_mus -gt 1 ]; then
+    # Create the non-equilbrium contour
+    echo "%block TS.Contours.nEq"
+    for i in `seq 1 $_mus` ; do
+	[ $i -eq $_mus ] && continue
+	echo "  neq-$i"
+    done
+    echo "%endblock TS.Contours.nEq"
+fi
+
+# Create array of chemical potentials sorted
+mus=()
+j=1
+for i in `seq $_mus -1 1` ; do
+    tmp="$(get_opt -mu$i-mu 1)"
+    tmp=${tmp//\|V\|/V}
+    tmp=${tmp//V/\|V\|}
+    mus[$j]=$(mu_e_correct $tmp)
+    let j++
+done
+
+
+# Create the contours
+if [ $_mus -gt 2 ]; then
+    echo "MSG: You can do with a single contour line in the bias window" >&2
+    echo "MSG: Use the lowest bias to the highest bias." >&2
+    echo "MSG: Here we supply as many different parts as there are chemical potentials" >&2
+fi
+
+if [ $_mus -gt 1 ]; then
+
+    i=1
+    echo "%block TS.Contour.nEq.neq-$i"
+    echo "  part line"
+    j=$((i+1))
+    if [ $j -eq $_mus ]; then
+	echo "   from ${mus[$i]} - 5 kT to ${mus[$j]} + 5 kT"
+    else
+	echo "   from ${mus[$i]} - 5 kT to ${mus[$j]}"
+    fi
+    echo "     delta $de eV"
+    echo "      method mid-rule"
+    echo "%endblock TS.Contour.nEq.neq-$i"
+
+    # Sort all chemical potentials
+    for i in `seq 3 $((_mus))` ; do
+	# Create the contours
+	echo "%block TS.Contour.nEq.neq-$((i-1))"
+	echo "  part line"
+	if [ $i -eq $_mus ]; then
+	    echo "   from prev to ${mus[$i]} + 5 kT"
+	else
+	    echo "   from prev to ${mus[$i]}"
+	fi
+	echo "     delta $de eV"
+	echo "      method mid-rule"
+	echo "%endblock TS.Contour.nEq.neq-$((i-1))"
+    done
+fi
+
+fi
+
+
+
+# Recreate the old 
+opts=( "${old_opts[@]}" )
+_N=${#opts[@]}
+# remove all mu
+i=0
+while [ $i -lt $_N ]; do
+    case ${opts[$i]} in
+	-mu*)
+	    unset opts[$i]
+	    let i++
+	    unset opts[$i] ;;
+    esac
+    let i++
+done
+# Remove empty things
+reduce
+
+
+# Now we will concentrate on the electrodes
+# First we need to figure out the number of electrode
+# I don't suspect ANY user will go above 10 electrodes
+_els=$(count el)
+[ $_els -eq 0 ] && \
+    error "You have zero electrodes." "Please supply at least one..."
+[ $_els -lt $_mus ] && \
+    error "You have fewer electrodes than chemical potentials." "Please correct your input..."
+
+
+if [ $print_el -eq 1 ]; then
+
+# Print out electrodes block:
+echo ""
+echo "%block $prefix.Elecs"
+for i in `seq 1 $_els` ; do
+    echo "  $(get_opt -el$i-name 1)"
+done
+echo "%endblock $prefix.Elecs"
+echo ""
+
+# Get whether or not we should print everything    
+_all=$(get_opt -print-all)
+[ $_all -eq 0 ] && _all=$(get_opt -el-all)
+
+
+function create_el {
+    local el=$1 ; shift
+    local name=$(get_opt -el$el-name 1)
+    [ ${#name} -eq 0 ] && error "Electrode: $el" "Has not received a name (-el$el-name)"
+    local chem=$(get_opt -el$el-mu 1)
+    [ ${#chem} -eq 0 ] && error "Electrode: $el" "Has not received a chemical potential reference (-el$el-mu)"
+    # Currently we don't check for the existance of
+    # the chemical potential.
+    local found=0
+    local i=$_mus
+    let i--
+    for j in `seq 0 $i` ; do
+	if [ "x${_mu_names[$j]}" == "x$chem" ]; then
+	    found=1
+	fi
+    done
+    [ $found -ne 1 ] && error "Electrode: $el could not be associated with a matching chemical potential: $chem"
+    echo "%block $prefix.Elec.$name"
+    local tshs=$(get_opt -el$el-tshs 1)
+    
+    if [ ${#tshs} -ne 0 ]; then
+	echo "  HS $tshs"
+    else
+	echo "  HS <input file>"
+    fi
+    echo "  chemical-potential $chem"
+    echo "  semi-inf-direction $(get_opt -el$el-semi-inf 1)"
+    local pos=$(get_opt -el$el-pos 1)
+    if [ ${#pos} -eq 0 ]; then
+	pos=$(get_opt -el$el-end-pos 1)
+	[ ${#pos} -ne 0 ] && pos="end $pos"
+    fi
+    [ ${#pos} -eq 0 ] && pos="<set value>"
+    echo "  electrode-position $pos"
+    # Now we can set all the other things (non-important)
+    function print_if {
+	local opt=
+	for opt in ${1//:/ } ; do
+	    local val=$(get_opt -el$el-$opt 1)
+	    if [ ${#val} -gt 0 ]; then
+		echo "  $2 $val"
+		break
+	    fi
+	done
+    }
+    print_if used-atoms:ua "used-atoms"
+    print_if bloch-a:bloch-a1 "bloch-a1"
+    print_if bloch-b:bloch-a2 "bloch-a2"
+    print_if bloch-c:bloch-a3 "bloch-a3"
+    print_if bulk "bulk"
+    print_if gf "GF"
+    echo "%endblock $prefix.Elec.$name"
+}
+for i in `seq 1 $_els` ; do
+    create_el $i
+done
+
+fi
+
+
+# Recreate the options
+opts=( "${old_opts[@]}" )
+_N=${#opts[@]}
+