search

geodynamics / Rayleigh

Rayleigh: Pseudo-spectral MHD
GNU General Public License v3.0
61 stars 48 forks source link

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

Closed bhishekmanek closed 10 months ago

bhishekmanek commented 1 year ago

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.