From 1240d65eaf3de114438185302a98be19c9ff22c7 Mon Sep 17 00:00:00 2001 From: Lindsey Date: Thu, 4 Jun 2015 16:12:29 -0700 Subject: [PATCH] B fields for casing mag dipole --- simpegEM/Analytics/FDEMcasing.py | 10 ++++++++-- 1 file changed, 8 insertions(+), 2 deletions(-) diff --git a/simpegEM/Analytics/FDEMcasing.py b/simpegEM/Analytics/FDEMcasing.py index 447f4986..547d40ac 100644 --- a/simpegEM/Analytics/FDEMcasing.py +++ b/simpegEM/Analytics/FDEMcasing.py @@ -87,5 +87,11 @@ def getCasingHrMagDipole(srcloc,obsloc,freq,sigma,a,b,mu=mu_0,eps=epsilon_0,mome def getCasingHzMagDipole(srcloc,obsloc,freq,sigma,a,b,mu=mu_0,eps=epsilon_0,moment=1): d2HertzZdz2 = _getCasingHertzMagDipole2Deriv_z_z(srcloc,obsloc,freq,sigma,a,b,mu,eps,moment) k2 = k(freq,sigma[2],mu,eps) - HertzZ(srcloc,obsloc,freq,sigma,a,b,mu,eps,moment) - return d2HertzZdz2 + k2**2 * HertzZ \ No newline at end of file + HertzZ = _getCasingHertzMagDipole(srcloc,obsloc,freq,sigma,a,b,mu,eps,moment) + return d2HertzZdz2 + k2**2 * HertzZ + +def getCasingBrMagDipole(srcloc,obsloc,freq,sigma,a,b,mu=mu_0,eps=epsilon_0,moment=1): + return mu_0 * getCasingHrMagDipole(srcloc,obsloc,freq,sigma,a,b,mu,eps,moment) + +def getCasingBzMagDipole(srcloc,obsloc,freq,sigma,a,b,mu=mu_0,eps=epsilon_0,moment=1): + return mu_0 * getCasingHzMagDipole(srcloc,obsloc,freq,sigma,a,b,mu,eps,moment) \ No newline at end of file