ignotur / magnetar_spectrum Goto Github PK

The way how k-vector is initialised (lines 141-143 of photon.cpp) seems to be wrong

Now I always emit single photon.

It is worth to integrate the optical depth only in the region where the scattering is allowed

The code should reproduce a benchmark from Fernandez & Thompson (2007), so the makefile should include an instruction according to which the code is compiled with proper f (beta gamma) distribution

The photons are probably not emitted correctly because mu value is nearly constant for fixed location at the NS surface

Code should be in a folder src

The code solve_twisted_magnetosphere.py should also print \Delta\phi N-S angle

For testing plot photon trajectory in mu-omega_c/omega plane