Add new flags 'p' and 'P' to set size and density scale at the poles
Make sure that we pick the correct GPU before configuring parameters
Performing additional checks in j_synchr() to avoid NaNs
Automatically change the data dimension according to the harm data
Use the black hole spin provided by a harm data file if it is provided
Cast constants to real to avoid the "Demoting to float" warning
Copy a_spin to GPU even if HARM is disabled
Rename j_nu to L_j_nu to indicate length scale is taken into account
Fix bug: make sure the cut-out region is computed correctly
Switch back to the original configuration mechanism for a_spin
Automatically load the locations of the radial boundaries
The values we read in are actually ln(r)
Compute the photon grid by input data such as lnrmin and lnrmax
Add a parameter 'O' to let the user change the output file names
cuda_gl_interop.h may cause errors; skip it when OpenGL is disable
Assume "usgdump2d" and the fieldline files are in the same directory
Fix bug: grid path is now correct even at the current directory
Output intensity and optical depth in single precision
Switch between a wind Tp/Te and a disk Tp/Te by looking at beta
Fix bug: correct the equation for magnetic field