Sets the density, temperature, gravity and magnetic field attributes of the respective fields to a realistic solar atmosphere profile. This routine first interpolates the temperature and numberdensity table at n_interp resolution, then solves the following ODE for the density: using a fifth order Runge-Kutta method. If the optional argument save_to is provided then the density profiles are saved to that file, which can be loaded back in on subsequent runs through the optional argument load_from. The integration is done over the entire table, the curve is sampled on the Gaussian grid, meaning that grid variations can all use the same result.
Throws an error if the geometry is not Cartesian.
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| type(settings_t), | intent(inout) | :: | settings | |||
| type(background_t), | intent(inout) | :: | background | |||
| type(physics_t), | intent(inout) | :: | physics | |||
| integer, | intent(in), | optional | :: | n_interp | points used for interpolation, defaults to 4000 if not present |