Rayleigh version:0.9.1 date:: 2018-04-28


Rayleigh-0.9.1.tar.gz [2018-04-28]

Primary Developers

Nicholas Featherstone
University of Colorado, Boulder


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 (CIG).

Cite as

Featherstone, N. (2018), Rayleigh 0.9.1, Available from:, doi: 10.5281/zenodo.1236565, url:

Primary References

Featherstone, N.A.; Hindman, B.W. (2016), The Spectral Amplitude Of Stellar Convection And Its Scaling In The High-Rayleigh-Number Regime, The Astrophysical Journal, 818 (1) , 32, doi: 10.3847/0004-637X/818/1/32, url:

Secondary References

Please also cite:

  • Matsui, H.; Heien, E.; Aubert, J.; Aurnou, J.M.; Avery, M.; Brown, B.; Buffett, B.A.; Busse, F.; Christensen, U.R.; Davies, C.J.; Featherstone, N.; Gastine, T.; Glatzmaier, G.A.; Gubbins, D.; Guermond, J.-L.; Hayashi, Y.-Y.; Hollerbach, R.; Hwang, L.J.; Jackson, A.; Jones, C.A.; Jiang, W.; Kellogg, L.H.; Kuang, W.; Landeau, M.; Marti, P.H.; Olson, P.; Ribeiro, A.; Sasaki, Y.; Schaeffer, N.; Simitev, R.D.; Sheyko, A.; Silva, L.; Stanley, S.; Takahashi, F.; Takehiro, S.-ichi; Wicht, J.; Willis, A.P. (2016), Performance benchmarks for a next generation numerical dynamo model, Geochemistry, Geophysics, Geosystems, 17 (5) , 1586-1607, doi: 10.1002/2015GC006159, url:


Primary Developer
Nicholas Featherstone
Primary Manual
Nicholas Featherstone
National Science Foundation 0949446
National Science Foundation 1550901

If you see errors or have further questions, please email

Sign In