The MPI - Adaptive Mesh Refinement - Versatile Advection Code (development version)
Frequently Asked Questions

List of questions

Below you can find a list of frequently asked questions and answers.

Is there a mailinglist?

Yes, we have a mailinglist that you can send questions to:

You can join the mailing-list by subscribing, so that you will be informed about important changes or (bug)fixes. You can also search the mailing list archive.

Is there a way I can define my own parameters somewhere in mod_usr.t and configure them through <tt>amrvac.par</tt> ?

Indeed, there is a quick and time-saving way to read your own parameters without having to give an explicit value in the usr file and recompile each time. Instead, add this in your usr file :

  1. at the end of the usr_init subroutine, add this line : call params_read(par_files)
  2. and just after the usr_init subroutine, define the params_read subroutine you just called. For instance :
!> Read parameters from a file
subroutine params_read(files)
character(len=*), intent(in) :: files(:)
integer :: n
namelist /my_list/ my_parameter_1, my_parameter_2 ! where you tell the code to read your own parameters
do n = 1, size(files)
open(unitpar, file=trim(files(n)), status="old")
read(unitpar, my_list, end=111)
111 close(unitpar)
end do
end subroutine params_read
This module contains definitions of global parameters and variables and some generic functions/subrou...
integer, parameter unitpar
file handle for IO

with my_parameter_1 and my_parameter_2 to be defined at the very beginning of the user file, just before the “contains” statement. Doing so, you can use them anywhere in the usr file and they will have the value you defined in your par file adding the following lines :

my_parameter_1 = 0.2d0
my_parameter_2 = 1.d0

I'm looking for a way to translate some .blk files (using ‘convert_type= 'onegrid’<tt>) back into standard output</tt>.dat` that I can use to restart a simulation. My goal is to be able to change my data mid-simulation through another program.

The idea of the onegrid option for conversion is that a hierarchical block AMR grid (as stored in the .dat files) can be saved to an equivalent uniform grid representation at a user-chosen level (combine level_io with onegrid). You then can use any software you like to handle uniform grid data. Converting back to the .dat format is then impossible. However, you can write a user-routine to read in the uniform data, and then use it to restart a similar uniform-grid (no AMR, just domain decomposed) simulation. The .dat files are all that is needed to do restarts.

If you want to change the data during the simulation, in principle you do not need another program to do that. For that purpose, we provide the generic subroutines

  1. usr_process_grid
  2. usr_process_global

or immediately after the advance step:

  1. usr_process_adv_grid
  2. usr_process_adv_global

or just once, at restart (and after reading in a .dat file), then use

  1. usr_before_main_loop

which allows you to modify things during runtime (of course, it should make physical sense). See their interface in mod_usr_methods.t.

Is there a way to save (and visualize) the data in the boundary ghost cells?

The logical variable save_physical_boundary can be set to true, which enforces the .dat files to also contain the ghost cell information for physical boundaries (i.e., those beyond the left or right edge of the domain, in each direction). You can use this file (like any other .dat file, to restart, and this helps if you want to use saved boundary info in your boundary value handling. However, all our present conversion options (like e.g. to .vtu files) do not store this extra info, and you can therefore not use them for visualizing the ghost cell info. For that, you will need to handle the .dat files directly, e.g. using python.

I get an error: mpi.mod was created by a different version of GNU Fortran

This means that your MPI library was compiled with a different version of gfortran (and GCC) than you are using now. If you run on a cluster, contact their support and ask them to fix it.

If you encounter this problem on your own machine, you can try to:

  1. Change compilers: Use the version of gfortran that MPI was compiled with
  2. Reinstall your MPI library, or install a different one (e.g. MPICH or OpenMPI), which is then hopefully compiled with the right version of GCC/gfortran.
  3. Install an operating system with working package management, e.g. Debian ;)

Can the MHD module handle Hall-MHD?

The MPI-AMRVAC code has Hall-MHD included, as detailed in some of our publications, e.g. in the method paper

  • ‘MPI-AMRVAC for Solar and Astrophysics’, O. Porth, C. Xia, T. Hendrix, S.P. Moschou, & R. Keppens, 2014, ApJS 214,4 doi:10.1088/0067-0049/214/1/4

or also in the Kelvin-Helmholtz related application paper

  • ‘On the influence of environmental parameters on mixing and reconnection caused by the Kelvin-Helmholtz instability at the magnetopause’, M.H.J. Leroy & R. Keppens, 2017, PoP 24, 012906 (15pp), doi:10.1063/1.4974758

The implementation details are given in the first reference (Porth et al.), and although it works properly on several tests and applications, we note that the time step constraint of our explicit implementation may become prohibitive for particular applications. We just limit the CFL according to \(\Delta t < \Delta x/ c_w\), with time step and spatial step linked by the speed \(c_w\), but in that speed we set \(c_w= |v|+\max(c_{fast}, \eta_h B/\rho * k_{max})\) and \(k_{max}\) is the maximal wavenumber we can represent (i.e. linked to \(\Delta x\)). The dispersive nature of the Hall-MHD system may then make \(\Delta t\) going like \(\Delta x^2\), and this limits the current implementation.

In MPI-AMRVAC 3.0, the Hall effect is included when setting the following in the mhd_list namelist part

mhd_etah= .... (some positive number, quantifying the hall parameter)

In the same namelist, the optional logical mhd_4th_order=.true. implies a 4th order evaluation for currents, default it is only second order. In any case, you may also need to activate an additional ghost cell layer (or 2 for 4th-order evaluations), through setting appropriately the parameter nghostcells.