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:
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: Import → Data → Magnetics → GIF format
For topography: Import → Data → Topography → 3D GIF (3D XYZ)
For the 3D mesh: Import → Mesh → 3D
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 Manipulation → Assign 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 Manipulation → Set 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:
Create → Processing → Magnetic equivalent source
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.
The info panel shows you which data item, topography, and mesh are being used. In addition, it will show you the sensitivity filename.
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:
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 source → Load results → Last 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:
Create → Forward modelling → Magnetics → MAG3D (v6.0)
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 Manipulation → Choose different data type → Magnetic amplitude data
The following dialog will appear to set the I/O headers:
Unset the data output using the drop-down menus.
Click “OK”. We now have a MAGAMPdata item in the GIFtools project tree.
Select the forward modelling item. Use the menus to edit its options:
Magnetics forward modelling → Edit 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 modelling → Write files
Then use the menu to run the forward model:
Magnetics forward modelling → Run 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 modelling → Load → Predicted 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