From f3bd31b25826ae25ba4cd9aea60ffd4bef9897ac Mon Sep 17 00:00:00 2001 From: rowanc1 Date: Sat, 12 Jul 2014 11:52:25 -0500 Subject: [PATCH] added current to the dipole transmitter --- simpegDC/BaseDC.py | 5 ++++- 1 file changed, 4 insertions(+), 1 deletion(-) diff --git a/simpegDC/BaseDC.py b/simpegDC/BaseDC.py index d9a7b7b6..929f226c 100644 --- a/simpegDC/BaseDC.py +++ b/simpegDC/BaseDC.py @@ -3,6 +3,9 @@ from SimPEG import * class DipoleTx(Survey.BaseTx): """A dipole transmitter, locA and locB are moved to the closest cell-centers""" + + current = 1 + def __init__(self, locA, locB, rxList, **kwargs): super(DipoleTx, self).__init__((locA, locB), 'dipole', rxList, **kwargs) self._rhsDict = {} @@ -12,7 +15,7 @@ class DipoleTx(Survey.BaseTx): pts = [self.loc[0], self.loc[1]] inds = Utils.closestPoints(mesh, pts) q = np.zeros(mesh.nC) - q[inds] = - ( np.r_[1., -1.] / mesh.vol[inds] ) + q[inds] = - self.current * ( np.r_[1., -1.] / mesh.vol[inds] ) self._rhsDict[mesh] = q return self._rhsDict[mesh]