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

Fortran source codes for Boussinesq systems

• bsnsq_rot-3d-shell_ut_mpi_module_sypack_abcn_hdiff_Omega.f90
  • A source code for followup calculations of Heimpel and Aurnou (2007) (Takehiro et al. 2024).
  • Designed for the whole spherical domain calculations.
  • Optional setting for the benchmark calculation (Case 0 of Chistensen et al. 2001).
  • Spmodel library is needed to compile and execute.
  • A sample script for exectuion
• Scripts for dynamo benchmark case0
  • Analysis : dynamobench_case0_analysis.rb
  • Drawings : 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.
• Heimpel, M., Aurnou, J., 2007 : Turbulent convection in rapidly rotating spherical shells: A model for equatorial and high latitude jets on Jupiter and Saturn. Icarus, 187, 540--557.
• Takehiro, S., Sasaki, Y., Ishioka, K., Enomoto, T., Nakajima, K., Hayashi, Y.-Y., 2024: Asymptotic profiles of mean zonal flows generated by thermal convection of Boussinesq fluid in a rapidly rotating thin spherical shell, submitted to Icarus.