User Tools

Site Tools


software:mag:start

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revision Previous revision
Next revision
Previous revision
software:mag:start [2013/12/31 00:15]
emheien
software:mag:start [2020/12/29 01:00] (current)
ljhwang
Line 1: Line 1:
 ====== MAG ====== ====== MAG ======
 +
 +//This content has been moved to hubzero. LJ Hwang 2020.12.28//​
  
 [[http://​geodynamics.org/​cig/​software/​mag|MAG Software Page]] [[http://​geodynamics.org/​cig/​software/​mag|MAG Software Page]]
Line 12: Line 14:
  
 A high-resolution version of the above phenomenon. ({{:​software:​mag:​magmv500.mov|original hi-res}}, {{:​software:​mag:​magmv315l.mov|original low-res}}) A high-resolution version of the above phenomenon. ({{:​software:​mag:​magmv500.mov|original hi-res}}, {{:​software:​mag:​magmv315l.mov|original low-res}})
 +
 +===== Tutorials =====
 +
 +The following tutorials demonstrate (1) how to run the prepackaged benchmark data, and (2) how to produce a magnetic field reversal.
 +
 +==== Prerequisites and Run Notes ====
 +
 +  * **Install MAG** - The following tutorials assume you have already installed MAG following the instructions in Chapter 2 in the [[http://​geodynamics.org/​cig/​software/​mag/​mag_book.pdf|MAG User Manual]].
 +
 +  * **Delete, Move, or Rename Output Files** - After each run, you must delete, move, or rename all of the output files in the current directory before re-running with the same "​output"​ filename. Retaining same-named output files in the current directory causes MAG to crash. ​
 +
 +  * **Changing Parameters** - Physical and time-step parameters can be changed in the par-file namelist without re-compiling MAG. See Appendix B in the [[http://​geodynamics.org/​cig/​software/​mag/​mag_book.pdf|MAG User Manual]] for a list of the input parameter names and definitions. Grid parameters must be changed in param.f and MAG must be then re-compiled. There are some numerical restrictions on the grid parameter combinations,​ which are given in Appendix A. 
 +
 +==== Tutorial 1: Running Benchmarks ====
 +Uncompress all files, and create a path (see Chapter 2 of [[http://​geodynamics.org/​cig/​software/​mag/​mag_book.pdf|MAG User Manual]]) ​
 +Link the grid parameter file to param.f, which enters into most subroutines through "​include"​ statements. For example, a grid parameter file named '​param32f4.f'​ (32 is the spherical harmonics truncation degree; 4 is the longitude symmetry) is linked using:
 +
 +    $ ln -sf param32s4.f param.f ​
 +
 +Compile the program as follows (renaming the executable in line 2 is optional):
 +
 +    $ make
 +    $ mv magx magx32s4  ​
 +
 +MAG uses a standard input file. Background execute using par.XXX as the input file and .YYY as the output file's extension, e.g., '$ magx32s4 <par.XXX >p.YYY &',​ so to run with the benchmark input files (par.bnch0 or par.bnch1), the execution statement should be:
 + 
 +    $ magx32s4 <​par.bnch0 >​p.bench0 &
 +
 +If there is a problem with the input file list, it is often the final three lines; with some systems, a "​$"​ may be required at the end.
 +
 +MAG produces a series of output files. For example, when using input file par.bnch0 (the example in step 4) MAG generates: '​l.bench0,​ ls.bench0, g[i].bench0 and d[i].bench0',​ where i=0,​1,​2...9. See Appendix B in the [[http://​geodynamics.org/​cig/​software/​mag/​mag_book.pdf|MAG User Manual]] for details on MAG's output files. Compare your output files with the data provided in the directory ~/​bench-data/​data_bench0.
 +
 +==== Tutorial 2: Reversal Dynamo Case ====
 +
 +In this example, we produce a magnetic field reversal using MAG. The input parameter in the source directory for this case is ~/​src/​par.Rev. There is no longitudinal symmetry in this case, so when you compile MAG, use param32s1.f linking to param.f. The Ekman number is E=0.02, the Prandtl number is Pr=1 and the magnetic Prandtl number is Pm=10. The Rayleigh number is Ra=12000.
 +
 +Link param32s1.f ​ to param.f:
 +
 +    $ ln -sf param32s1.f ​ param.f ​
 +
 +Compile the program with:
 +
 +    $ make
 +    $ mv magx magx32s4 ​
 +
 +Background execute using par.Rev as the input file and p.rev as the output file (Note: The output file's extension is specified in the par-file, e.g., here par.Rev specifies that the outfile should be named *filename*.rev):​
 +
 +    $ magx32s4 <par.Rev >p.rev &
 +
 +MAG will then execute and produce datafiles in the same directory. See Appendix B in the [[http://​geodynamics.org/​cig/​software/​mag/​mag_book.pdf|MAG User Manual]] for details on MAG's output files. ​
 +
 +=== Results and Discussion ===
 +
 +This case has run on 32-bit and 64-bit Intel processors. Figure 1 below shows a plot of mean velocity Vrms, mean magnetic field Brms, the axial dipole and the dipole tilt on the outer boundary. It indicated a magnetic field reversal between time steps 25 and 30. 
 +
 +{{:​software:​mag:​tutorial:​field-64.gif?​300|Figure 1}}
 +
 +Figure 2 shows a longer run of MAG, where we see the magnetic field reversed again. At this time, the magnetic field had weakened substantially. ​
 +
 +{{:​software:​mag:​tutorial:​field-64-revr.gif?​300|Figure 2}}
 +
 +In Figures 3 and 4 (below), the top figure is the pole plot before the second field reversal and the bottom is the pole plot after the second field reversal. ​
 +
 +{{:​software:​mag:​tutorial:​g1revr.gif?​300|Figure 3}}
 +{{:​software:​mag:​tutorial:​g7revr.gif?​300|Figure 4}}
 +
 +An additional tutorial on creating a Reversal Dynamo movie is found in Chapter 5 of the [[http://​geodynamics.org/​cig/​software/​mag/​mag_book.pdf|MAG User Manual]].
software/mag/start.1388448945.txt.gz ยท Last modified: 2013/12/31 00:15 by emheien