8.2.2. Forward model magnetic amplitude data in GIFtools

In this recipe, we step through forward modeling magnetic amplitude data, starting from a TMI data file, a topography file, and a 3D mesh. At the bottom of this page, you will find the files used in this recipe as well as the complete project and output files.

This recipe requires the following steps:

  1. Import the necessary files

  2. Invert for an equivalent source model

  3. Forward model amplitude data

STEP 1: Import the necessary files

Start GIFtools or open an already existing project. Because we will be loading in a few files, it is easiest to set the working directory to the folder where the files are located and/or where you would like to do the forward model.

The first step is then to import the magnetic data, the topography file, and the 3D mesh. Importation of data can slightly differ, depending on the file type (GIF format, XYZ, etc). For this example, the magnetic data and the topography file are both in GIF format. Both are imported using the following menus:

  • For magnetic survey locations: ImportDataMagnetics → GIF format

  • For topography: ImportDataTopography → 3D GIF (3D XYZ)

  • For the 3D mesh: ImportMesh3D

../../../_images/importdata.png

Note

Note that the mesh does not have any padding in the vertical direction, as required for the equivalent source problem.

Once the data are loaded, GIFtools will have the three items shown in the project tree on the left. It’s a good idea to regularly save your project!

Note

If you have a susceptibility model, you can import the model and skip directly to forward modeling the amplitude data.

Before we can solve for an equivalent source model for this magnetic data, we need to set the I/O for the data and uncertainties. This means we need to assign uncertainties (link needed).

Use the menu structure to set the uncertainty for the magnetic data:

Data ManipulationAssign uncertainties

A dialog pops up to add a percentage and/or a floor to each header. We completed the dialog by clicking the Assign checkbox for the header “B” and using 0 for percentage and 0.5 nT for the floor. Click “OK” to apply the changes. Now use the menu structure to set the I/O headers for the data:

Data ManipulationSet I/O headers

For this example, make sure the total-field anomaly is set the “B” and the uncertainty to “B_uncert” by choosing these from the drop-down boxes. Click “OK” to close the dialog and apply the changes. The info panel for the magnetic data item will have updated and all minimum data requirements for the inversion are now set (indicated using ** in the info panel).

STEP 2: Invert for an equivalent source model

Now we can invert for an equivalent source model. To create the equivalent source item, use the menu structure in GIFtools as following:

CreateProcessingMagnetic equivalent source

../../../_images/step2.png

Choose a directory to write files to and run the equivalent source in. Click “Open”. You will now see the equivalent source item in the GIFtools project tree on the left. For this example, it’s called “MagEquivalentSource (MAGequivalentSource)”. You can always rename the first part of the name.

../../../_images/step21.png

The info panel shows you which data item, topography, and mesh are being used. In addition, it will show you the sensitivity filename.

../../../_images/step2a.png

The equivalent source item has its own menu which contains the following:

  • Working directory (to set or view)

  • Edit options

  • Copy options

  • Write

  • Run

  • Load results

It also has a visualization menu. These menus are identical to those for an inversion item so we refer you to the inversion recipe for more information on each one.

Use the menu and select Edit options. The dialog shown below pops up:

../../../_images/step2b.png

See the MAGES manual (which can be accessed via the GIF documentation page too) for detailed information about each parameter.

In this example, set the mesh, data, and topography. We will leave the remaining parameters set to the default values. Click “Apply and write files” to set the changes and write the files to the working directory. A message will appear to state that the files have been written. Click “OK”.

Next, use the menu to run the equivalent source. Use the “all” options to run both the sensitivity and the inversion. A command window will pop up and the code will run. Once finished, close the command window. In GIFtools, use the menu to load the last result:

Magnetic equivalent sourceLoad resultsLast result. The equivalent source item will become a folder, with two subfolders (one for predicted data and one for the models).

If you view the model (use the visualization menu for the recovered GIFmodel), you’ll see it is a single layer, just below the topography.

STEP 3: Forward model amplitude data

The final step is to forward model amplitude data from the equivalent source model. Using the menus, create a forward modeling item:

CreateForward modellingMagneticsMAG3D (v6.0)

../../../_images/step3.png

Select the magnetic data (MAGdata) item that we used for the equivalent source calculation. Convert it to an amplitude data item using the menus:

Data ManipulationChoose different data typeMagnetic amplitude data

../../../_images/step3a.png

The following dialog will appear to set the I/O headers:

../../../_images/step3b.png

Unset the data output using the drop-down menus.

../../../_images/step3c.png

Click “OK”. We now have a MAGAMPdata item in the GIFtools project tree.

../../../_images/step3d.png

Select the forward modelling item. Use the menus to edit its options:

Magnetics forward modellingEdit options

The dialog is the exact same as in the forward modeling recipe. Select the equivalent source model, the amplitude data item, and the topography. Then click “Apply”.

Use the menu to write the files:

Magnetics forward modellingWrite files

Then use the menu to run the forward model:

Magnetics forward modellingRun MAGFOR3D

A command window will open up and the code will start running. Once the code has finished, close the command window. Finally, use the menu to import the predicted data into GIFtools:

Magnetics forward modellingLoadPredicted data

The forward modelling item will change into a folder, with the predicted data as an item within that folder. Select the predicted data and view it to see the amplitude data.

Example

The magnetics data file, topography file, and 3D mesh along with a GIFtools project to repeat the above steps to forward model amplitude data: download