This PR will allow users to pass custom functions for plasma conditions and reaction rates as keyword arguments to either Stellarator or SourceMesh. The default behavior is to use the same plasma condition/reaction rate that it always has.
The plasma conditions function takes only one input, the plasma parameter s, and must return whatever inputs are required by the reaction rate function. The reaction rate function should return a single value.
Here's an example where the source strengths follow a sine profile:
Should I extend the parastell_example.py script to include an example of creating a custom source definition as well?
This PR will allow users to pass custom functions for plasma conditions and reaction rates as keyword arguments to either
StellaratororSourceMesh. The default behavior is to use the same plasma condition/reaction rate that it always has.The plasma conditions function takes only one input, the plasma parameter s, and must return whatever inputs are required by the reaction rate function. The reaction rate function should return a single value.
Here's an example where the source strengths follow a sine profile:
Should I extend the parastell_example.py script to include an example of creating a custom source definition as well?
Addresses #51