Skip to content

Commit

Permalink
Finalize example
Browse files Browse the repository at this point in the history
  • Loading branch information
dccowan committed Oct 6, 2020
1 parent 27036df commit f0b0d93
Show file tree
Hide file tree
Showing 19 changed files with 16 additions and 17 deletions.
Binary file added _ext/__pycache__/edit_on_github.cpython-38.pyc
Binary file not shown.
4 changes: 1 addition & 3 deletions conf.py
Original file line number Diff line number Diff line change
Expand Up @@ -35,9 +35,7 @@
extensions = [
'sphinx.ext.mathjax',
'sphinxcontrib.bibtex',
'edit_on_github',
'purpose',
'example'
'edit_on_github'
]

# Number figures
Expand Down
3 changes: 1 addition & 2 deletions content/elements.rst
Original file line number Diff line number Diff line change
Expand Up @@ -31,7 +31,6 @@ There are seven general files which are used in MAG3D. All are in ASCII text for
Mesh <files/meshfile>
Topography <files/topo>
Observation/Location <files/magfile>
Model <files/model>
Model <files/modelfile>
Weighting <files/weight>
Active model <files/model>

10 changes: 6 additions & 4 deletions content/example/fwd.rst
Original file line number Diff line number Diff line change
@@ -1,19 +1,21 @@
.. _example_fwd:


.. note:: This example has been developed to demonstrate functionality specific to v6.0. The example can be completed using v5.0, however some functionality may not exist in v5.0 and the format of certain input files may differ slightly.

Forward Modeling
================

Here the code **gzfor3d.exe** is used to forward model the gravity anomaly for the mesh and density contrast model provided. We consider a ground-based survey with a uniform station spacing of 40 m.
Here the code **magfor3d_v60.exe** is used to forward model the magnetic anomaly for the mesh and susceptibility model provided. We consider an airborne survey with a uniform station spacing of 40 m and a flight height of 30 m. The inclination, declination and background field intensity are 60 degrees, 25 degrees and 50,000 nT, respectively.

Files relevant to this part of the example are in the sub-folder *fwd*. We used the same model that was created in the :ref:`create model <example_model>` section. Before running this example, you may want to do the following:

- `Download and open the zip folder containing the entire grav3d example <https://github.com/ubcgif/grav3d/raw/master/assets/grav3d_example.zip>`__ (if not done already)
- Learn how to run :ref:`gzfor3d<gzfor3d>`
- `Download and open the zip folder containing the entire mag3d example <https://github.com/ubcgif/mag3d/raw/master/assets/mag3d_example.zip>`__ (if not done already)
- Learn how to run :ref:`magfor3d_v60 <magfor3d>`
- There is no input file


The resulting gravity anomaly (mgal) is shown below. According to the convention used by *grav3d*, a positive density contrast results in a positive gravity anomaly.
The total magnetic intensity anomaly (nT) is shown below.


.. figure:: images/fwd_anomaly.png
Expand Down
Binary file modified content/example/images/distance_weights.png
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Binary file modified content/example/images/fwd_anomaly.png
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Binary file modified content/example/images/inv_L2_input.png
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Binary file modified content/example/images/inv_sparse_input.png
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Binary file modified content/example/images/model_L2.png
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Binary file modified content/example/images/model_input.png
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Binary file modified content/example/images/model_sparse.png
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Binary file modified content/example/images/pfweight_input.png
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Binary file modified content/example/images/sensitivity_L2_input.png
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Binary file modified content/example/images/sensitivity_sparse_input.png
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Binary file modified content/example/images/true_model.png
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
8 changes: 4 additions & 4 deletions content/example/inv.rst
Original file line number Diff line number Diff line change
Expand Up @@ -5,16 +5,16 @@
Inversion
=========

Here the code **gzinv3d_v60.exe** invert synthetic gravity anomaly data to recover a density contrast model. We perform both a smooth and a sparse inversion. We use the same data that were computed in :ref:`forward modeling example <example_fwd>`. Gaussian noise were added to these data with a standard deviation of 5e-4 mgal. Uncertainties of 5e-4 mgal were assigned to all data
Here the code **maginv3d_v60.exe** invert synthetic magntic data to recover a susceptibility model. We perform both a smooth and a sparse inversion. We use the same data that were computed in :ref:`forward modeling example <example_fwd>`. Gaussian noise were added to these data with a standard deviation of 0.5 nT. Uncertainties of 0.5 nT were assigned to all data

Before running this example, you may want to do the following:

- `Download and open the zip folder containing the entire grav3d example <https://github.com/ubcgif/grav3d/raw/master/assets/grav3d_example.zip>`__ (if not done already)
- Learn how to :ref:`run gzinv3d_v60 and learn the format of the input files <gzinv3d>`
- `Download and open the zip folder containing the entire mag3d example <https://github.com/ubcgif/mag3d/raw/master/assets/mag3d_example.zip>`__ (if not done already)
- Learn how to :ref:`run maginv3d_v60 and learn the format of the input files <maginv3d>`



.. important:: Since the sensitivities output by **gzsen3d_v60.exe** produce a large file, we have not provided them in the zip file. You need to complete the :ref:`compute sensitivities <example_sensitivity>` example to run the inversion.
.. important:: Since the sensitivities output by **magsen3d_v60.exe** produce a large file, we have not provided them in the zip file. You need to complete the :ref:`compute sensitivities <example_sensitivity>` example to run the inversion.


Smooth Inversion
Expand Down
6 changes: 3 additions & 3 deletions content/example/sensitivity.rst
Original file line number Diff line number Diff line change
Expand Up @@ -6,12 +6,12 @@
Computing Sensitivities
=======================

Here the code **gzsen3d_v60.exe** is used to compute and output the sensitivities for use in the inversion. Here, we generate we use the input file **sensitivity_L2.inp** to generate sensitivities for a standard smooth inversion. We then use the input file **sensitivity_sparse.inp** to generate sensitivities that can be used for sparse inversion.
Here the code **magsen3d_v60.exe** is used to compute and output the sensitivities for use in the inversion. Here, we generate we use the input file **sensitivity_L2.inp** to generate sensitivities for a standard smooth inversion. We then use the input file **sensitivity_sparse.inp** to generate sensitivities that can be used for sparse inversion.

Before running this example, you may want to do the following:

- `Download and open the zip folder containing the entire grav3d example <https://github.com/ubcgif/grav3d/raw/master/assets/grav3d_example.zip>`__ (if not done already)
- Learn how to :ref:`compute sensitivities and learn the format of the input files <gzsen3d>`
- `Download and open the zip folder containing the entire mag3d example <https://github.com/ubcgif/mag3d/raw/master/assets/mag3d_example.zip>`__ (if not done already)
- Learn how to :ref:`compute sensitivities and learn the format of the input files <magsen3d>`


Smooth Inversion
Expand Down
2 changes: 1 addition & 1 deletion content/example/weights.rst
Original file line number Diff line number Diff line change
Expand Up @@ -6,7 +6,7 @@ Depth/Distance Weighting

Here the code **pfweight.exe** and the input file **pfweight.inp** is used to generate distance weighting for the inversion. Files relevant to this part of the example are in the sub-folder *pfweights*. Before running this example, you may want to do the following:

- `Download and open the zip folder containing the entire grav3d example <https://github.com/ubcgif/grav3d/raw/master/assets/grav3d_example.zip>`__ (if not done already)
- `Download and open the zip folder containing the entire mag3d example <https://github.com/ubcgif/mag3d/raw/master/assets/mag3d_example.zip>`__ (if not done already)
- Learn how to :ref:`run pfweight and learn the format of the input files <pfweight>`


Expand Down
File renamed without changes.

0 comments on commit f0b0d93

Please sign in to comment.