It also looks like in many cases the code runs into the divide-by-0 issue in calculating factor in the _gain_harmonic_n function in cyclotronphysics.py file. This seems to be due to beta_ortho being calculated to 0 in _get_analytic_power function when the pitcjh angle is 0 and the issue arises since it’s used in the denominator in calculating 'factor'. Need to take a closer look.
It also looks like in many cases the code runs into the divide-by-0 issue in calculating factor in the _gain_harmonic_n function in cyclotronphysics.py file. This seems to be due to beta_ortho being calculated to 0 in _get_analytic_power function when the pitcjh angle is 0 and the issue arises since it’s used in the denominator in calculating 'factor'. Need to take a closer look.