From 31da380804f5ddf499277557738205f9781fbfbc Mon Sep 17 00:00:00 2001 From: rowanc1 Date: Mon, 10 Mar 2014 23:01:49 -0700 Subject: [PATCH] deleted example in fdem, use test for example for now. --- simpegEM/FDEM/FDEM.py | 46 ------------------------------------------- 1 file changed, 46 deletions(-) diff --git a/simpegEM/FDEM/FDEM.py b/simpegEM/FDEM/FDEM.py index 435990cf..af27bf09 100644 --- a/simpegEM/FDEM/FDEM.py +++ b/simpegEM/FDEM/FDEM.py @@ -123,52 +123,6 @@ class ProblemFDEM_e(Problem.BaseProblem): raise NotImplementedError('Jtvec todo!') -if __name__ == '__main__': - from SimPEG import * - import simpegEM as EM - from simpegEM.Utils.Ana import hzAnalyticDipoleT - from scipy.constants import mu_0 - import matplotlib.pyplot as plt - - cs = 5. - ncx = 6 - ncy = 6 - ncz = 6 - npad = 3 - hx = Utils.meshTensors(((npad,cs), (ncx,cs), (npad,cs))) - hy = Utils.meshTensors(((npad,cs), (ncy,cs), (npad,cs))) - hz = Utils.meshTensors(((npad,cs), (ncz,cs), (npad,cs))) - mesh = Mesh.TensorMesh([hx,hy,hz]) - - XY = Utils.ndgrid(np.linspace(20,50,3), np.linspace(20,50,3)) - rxLoc = np.c_[XY, np.ones(XY.shape[0])*40] - - model = Model.LogModel(mesh) - - opts = {'txLoc':0., - 'txType':'VMD_MVP', - 'rxLoc': rxLoc, - 'rxType':'bz', - 'freq': np.logspace(0,3,4), - } - survey = EM.FDEM.SurveyFDEM(**opts) - - prb = EM.FDEM.ProblemFDEM_e(mesh, model) - prb.pair(survey) - - sigma = np.log(np.ones(mesh.nC)*1e-3) - - j_sx = np.zeros(mesh.vnEx) - j_sx[6,6,6] = 1 - j_s = np.r_[Utils.mkvc(j_sx),np.zeros(mesh.nEy+mesh.nEz)] - - prb.j_s = j_s - f = prb.fields(sigma) - - plt.colorbar(mesh.plotSlice((f.get_e(3)), 'E', ind=11, normal='Z', view='real')[0]) - plt.show() - -