Project of convection in rotating spheres and spherical shells 〜 Fortran source codes for anelastic systems

Fortran source codes

• A source code for followup calculations of anelasctic benchmark and execution script (Jones et al., 2011)
  • spmodel library ver.2 is needed to compile and execute.
  • Non-dimensional formulation using rotatiing angular velocity.
  • Hyper viscosity and diffusion are implemented.
  • Exectution script for Boussinesq benchmark (Case 0 of Chistensen et al. 2001)

Scripts for data analysis and visualization

• For anelastic benchmark calculatoin (Jones et al., 2011)
  • Visualization : anelastic_dynamobench_hydro_plot.rb
• For Boussinesq benchmark calculation (Case 0 of Chistensen et al. 2001)
  • Analysis : anelastic_dynamobench_case0_analysis.rb
  • Visualization : anelastic_dynamobench_case0_plot.rb

References

• Christensen, U. R., Aubert, J., Cardin, P., Dormy, E., Gibbons, S., Glatzmaier, G. A., Grote, E., Honkura, Y., Jones, C., Kono, M., Matsushima, M., Sakuraba, A., Takahashi, F., Tilgner, A., Wicht, J., Zhang, K., 2001: A numerical dynamo benchmark. Phys. Earth Planet. Inter., 128, 25--34.
• Jones, C. A., Boronski, P., Brun, A. S., Glatzmaier, G. A.,Gastine, T., Miesch, M. S., Wicht, J., 2011 : Anelastic convection-driven dynamo benchmarks. Icarus , 216, 120--135.