Rayleigh is a 3-D convection code designed for the study of dynamo behavior in spherical shell geometry. It evolves the incompressible and anelastic MHD equations in spherical geometry using a pseudo-spectral approach.
Rayleigh employs spherical harmonics in the horizontal direction and Chebyshev polynomials in the radial direction. The code has undergone extensive accuracy testing using the Christensen et al. (2001) Boussinesq benchmarks and the Jones et al. (2011) anelastic benchmarks. Rayleigh has been developed with NSF support through the Computational Infrastructure for Geodynamics.
Bug-fix release: Version 1.0.1
Radially varying magnetic diffusivity is now implemented correctly using the Crank-Nicholson scheme.
We released Rayleigh v1.0.0 on November 12 with various new features and bugfixes. Among others, this release includes the following significant changes:
main_input, are now stored within each numbered checkpoint directory.
Nicholas Featherstone, Philipp Edelmann, Rene Gassmoeller, Loren Matilsky, Ryan Orvedahl, Cian Wilson, and many other contributors.
Full documentation is available at https://rayleigh-documentation.readthedocs.io/en/latest/index.html
NSF EAR-0949446 and EAR-1550901.