Accuracy problems about empymod

[prisae]

Hi Zhengguan,

Thanks for getting in touch. I run your example, and I cannot see any issues with it. What is, in your opinion, incorrect in the result?

Just to confirm a few things, you have:

• a 100 m long, x-directed electric dipole source, centred at x=y=z=0;
• 41 x-directed electric point-receivers, from (x, y, z) = (-10 km, +5 km, 0) to (+10 km, +5 km, 0).

If you can tell me what you expect to be different and why I might be able to help more.

I can see two potential points of confusion.

strength

One thing that might be confusing is the source strength. From the docstring:

strength : float, default: 0.0
    Source strength (A):

    - If 0, output is normalized to source and receiver of 1 m length, and
      source strength of 1 A.
    - If != 0, output is returned for given source and receiver length, and
      source strength.

You provide strength=1, which means that the returned result corresponds to a 100 Am source (100 m long, 1 A strong).

srcpts

The number of source points is by default set to 1. That means that your 100 m long dipole source is approximated by a single point source at (0,0,0). If you set srcpts=3, for instance, it uses Gaussian Quadrature to approximate the finite length dipole with three points.

However, given your source-receiver distances this will hardly make any difference at all, and approximating it with a point source is absolutely fine.

[liuzg-23]

Hi Dieter,

Thanks very much for your detailed response. You are so kind. According to your answers, I have checked my code again. I made a particularly small mistake about my code.

Best Regards,
Zhengguang.

[prisae]

You're welcome Zhengguang. I assume it works now, great!