Skip to content

Arbitrarily shaped sources
Compare
Choose a tag to compare
@prisae prisae released this 13 Jan 16:30
v0.16.0
ece2baf
This commit was created on GitHub.com and signed with GitHub’s verified signature. The key has expired.
GPG key ID: 4AEE18F83AFDEB23
Expired
Learn about vigilant mode.

  • fields.get_source_field:

    • Arbitrarily shaped sources (and therefore also loops) can now be created by
      providing a src that consists of x-, y-, and z-coordinates of all
      endpoints of the individual segments.

    • Simple "magnetic dipole" sources can now be created by providing a point
      dipole ([x, y, z, azm, dip]) and set msrc=True. This will create a
      square loop of lengthxlength m perpendicular to the defined point
      dipole, hence simulating a magnetic source. Default length is 1 meter.

    • Point dipoles and finite length dipoles were before treated differently.
      Point dipoles are new converted into finite length dipoles of provided
      length (default is 1 meter), and treated as finite length dipoles. This is
      backwards incompatible and means that the source field for point dipoles
      might not be exactly the same as before. However, in any properly set-up
      simulation this should have no influence on the result.

    • Bugfix: Fix floating point issue when the smaller coordinate of a finite
      length dipole source was very close to a node, but not exactly. This is
      done by rounding the grid locations and source position, and the precision
      can be controlled via decimals; default is micrometer.

  • fields: Values outside the grid in get_receiver and
    get_receiver_response are new set to NaN's instead of zeroes.
    Additionally, the first and last values in each direction of the fields are
    ignored, to avoid effects form the boundary condition (receivers should not
    be placed that close to the boundary anyway).

  • simulations:

    • Within the automatic gridding the properties are estimated much more
      conservative now, if not provided: before the log10-average of the last
      slice in a given direction was used; now it uses the maximum resistivity.
      This is usually the air value for x/y and positive z. This is very
      conservative, but avoids that users use too small computational domains in
      the case of land and shallow marine surveys. The downside is that it
      heavily over-estimates the required domain in the deep marine case.
      However, slower but safe is better in this case.
    • New method print_grids, which prints the info of all created grids.
      This is also used for logging in the CLI interface.

  • maps: interp3d takes a new keyword cval, which is passed to
    map_coordinates.

Assets 2
Loading