The Llama code is a 3-dimensional multiblock infrastructure for Cactus based on Carpet. It provides different patch systems that cover the simulation domain by a set of overlapping patches. Each of these patches has local cooordinates with a well-defined relation to global Cartesian coordinates. Information between the different patches is communicated via interpolation in the overlap zones.
Here we show an example evolution of a wave equation on a Kerr background using the "Thornburg2004nc" patch-system, consisting of 90 degree by 90 degree wedges centred on each of the +x, -x, +y, -y, +z, -z coordinate axes. This provides a complete covering of a region between rmin and rmax with smooth inner and outer boundaries. The initial data is an l=2, m=2 spherical harmonic in the angular directions, and a Gaussian in the radial direction. An alternative patch system, "Thornburg2004", adds a cubical patch enclosing r < rmin in which standard Carpet box-in-box mesh refinement can be employed. This is used in the GW150914 gravitational wave gallery example.
We ask that if you make use of the Llama code, then please cite Llama, the Einstein Toolkit, the Carpet mesh-refinement driver and Cactus.
Parameter File | Kerr-Schild_Multipole.par |
---|---|
Submission command | simfactory/bin/sim create-submit Kerr-Schild_Multipole --parfile arrangements/Llama/LlamaWaveToy/par/Kerr-Schild_Multipole.par |
Total memory | 800 MB |
Run time | A little over 1 minute using 56 MPI tasks and 56 cores on Frontera; or about 10 minutes on two laptop cores |
Results (gzip 226MB, uncompressed 460MB) | Kerr-Schild_Multipole-20240620.tar.gz |
This example was last tested on 20-June-2024.
visit
executable.File > Restore session with sources ...
.SOURCE00
is selected in the Source Identifiers
pane, click on the ...
button of the Source
pane.HDF5
files are located.Files
pane, locate and click on the u.file_* database
icon.OK
button.OK
button.If the file structure is different than expected in the saved *.session files above, the following VisIt directions might help:
visit
executable.Open
Directories
panel, navigate to the location of the simulation's HDF5
outputFiles
panel, navigate to the u.file_* database
entry and click it. The database will turn blue, indicating that it is selected.Open file as type
dropdown menu, select CarpetHDF5
. Click the OK
button.Plots
panel, click Add > Pseudocolor > LLAMAWAVETOY--u_lp_MP_rp_
.Plots
panel, click Operators > Slicing > Slice
.Plots
panel, open the dropdown menu of the Pseudocolor - Slice(LLAMAWAVETOY--u_lp_MP_rp_)
Pseudocolor - Slice(LLAMAWAVETOY--u_lp_MP_rp_)
, double-click the Slice
component.Slice operator attributes
window, under the Normal
panel, select the Z Axis
radio button. Click apply.Slice operator attributes
window.Plots
panel, click Draw
.Pseudocolor - Slice(LLAMAWAVETOY--u_lp_MP_rp_)
double-click on Pseudocolor
.Pseudocolor plot attributes
window, on the Data > Scale > Limits
panel, check the Minimum
and Maximum
checkboxes. Change the Minimum
field to -0.5
and Maximum
field to 0.5
.Color
panel, click the Color Table
dropdown menu and select the orangehot
item. Check the Invert
checkbox.Apply
and close the Pseudocolor plot attributes
window.Swap background and foreground colors
button, represented by a black triangle with two green arrows.Time
panel, type 56
on the time field and press the ENTER
key on your keyboard. The plot obtained reproduces exactly the first image on this gallery page. Users are encouraged to try and experiment with other settings, particularly other color maps and simulation times or 2D slices.File > Set save options ...
.Filename
panel, set the Output directory
field to the desired location. Click Apply
and close this window.File > Save window
to save the image.The Python script needs to be executed within VisIt's Python parser. Download the script and navigate to the directory where the visit binary is located (this is not necessary if VisIt is in your PATH
environment variable). In your command line, issue
visit -cli -nowin -s plot.py <HDF5 dir> <output dir> <output name>
where
<HDF5 dir>
is the directory containing the output HDF5
files.<output dir>
is the directory where the image will be saved.<output name>
is the name of the output file (excluding the file extension).