Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Initial conditions using checkpoint files added to user supplied files #488

Merged
merged 3 commits into from
Dec 5, 2023

Conversation

bhishekmanek
Copy link
Contributor

The default initial conditions type in Rayleigh allow reading from a checkpoint file (init_type = -1 for thermal and velocity fields, magnetic_init_type = -2 for magnetic fields). Modifications made in this request allow for new IC's:

  1. init_type = -2 reads checkpoint file as well as a user-supplied entropy profile. The sum total of these two is the new IC.
  2. magnetic_init_type = -2 reads checkpoint file and user-supplied profiles of toroidal and poloidal field components. The sum total of checkpoint file and user-supplied profiles is the new IC.

An additional subroutine Get_RHS is added to src/Math_Layer/Linear_Solve.F90.

bhishekmanek and others added 3 commits October 3, 2023 11:43
This file includes init_type = -2 for adding entropy perturbations from checkpoint and an user input file. Also included is a magnetic_init_type = -2 that adds checkpoints to user input files of poloidal and toroidal component of magnetic fields.
Copy link
Contributor

@feathern feathern left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Looks good. I'm going to approve and then submit a small PR of my own removing the !Nick comment.

@feathern feathern merged commit 7e4a084 into geodynamics:main Dec 5, 2023
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

2 participants