From f020db31466e770a6ca9b237e5291277e25b837f Mon Sep 17 00:00:00 2001
From: craig-warren <craig_warren@me.com>
Date: Thu, 7 Jan 2016 22:37:37 +0000
Subject: [PATCH] Changed name of source position variable in transmission
 line. Experimenting with different lengths of line.

---
 gprMax/sources.py | 9 +++++----
 1 file changed, 5 insertions(+), 4 deletions(-)

diff --git a/gprMax/sources.py b/gprMax/sources.py
index 0a1316f7..191b9e38 100644
--- a/gprMax/sources.py
+++ b/gprMax/sources.py
@@ -184,10 +184,11 @@ class TransmissionLine:
         self.dl = np.sqrt(3) * c * G.dt
         
         # Nodal position of one-way injector excitation in the transmission line
-        self.source = 10
+        self.srcpos = 10
         
         # Number of nodes in the transmission line; add nodes to the length to account for position of one-way injector
-        self.nl = rvalue(self.length/self.dl) + self.source
+#        self.nl = rvalue(self.length/self.dl) + self.srcpos
+        self.nl = 20
 
         self.voltage = np.zeros(self.nl, dtype=floattype)
         self.current = np.zeros(self.nl, dtype=floattype)
@@ -220,7 +221,7 @@ class TransmissionLine:
             self.voltage[i] -= self.resistance * (c * G.dt / self.dl) * (self.current[i] - self.current[i - 1])
         
         # Update the voltage at the position of the one-way injector excitation
-        self.voltage[self.source] += (c * G.dt / self.dl) * waveform.amp * waveform.calculate_value(time - 0.5 * G.dt, G.dt)
+        self.voltage[self.srcpos] += (c * G.dt / self.dl) * waveform.amp * waveform.calculate_value(time - 0.5 * G.dt, G.dt)
 
         # Update ABC before updating current
         self.update_abc(G)
@@ -240,7 +241,7 @@ class TransmissionLine:
             self.current[i] -= (1 / self.resistance) * (c * G.dt / self.dl) * (self.voltage[i + 1] - self.voltage[i])
         
         # Update the current one node before the position of the one-way injector excitation
-        self.current[self.source - 1] += (c * G.dt / self.dl) * waveform.amp * waveform.calculate_value(time - 0.5 * G.dt, G.dt) * (1 / self.resistance)
+        self.current[self.srcpos - 1] += (c * G.dt / self.dl) * waveform.amp * waveform.calculate_value(time - 0.5 * G.dt, G.dt) * (1 / self.resistance)
 
     def update_E(self, abstime, Ex, Ey, Ez, G):
         """Updates electric field value in the main grid from voltage value in the transmission line.