Load coord and field into global variables
Copy pointers into constant memory for the device to access
Merge branch 'develop' into harm

Conflicts:
	src/main.cc
Stop integration if photons are too far; dump data at the end if GL=0
Add a simple IDL script to visualize the image output
Compute red shift k_alpha u^alpha from the HARM data
Read the spin parameter from the coordinate data file
Set the spin parameter to the value read from the coordinate data
Update operation count
Update bandwidth estimate
Print the maximum density and energy
Construct the four vectors u and b in modified KS coordinates
Transform vector u and b from KSP to KS coordinates
Transform vector u and b from KS to BL coordinates
Compute red shift and angle cosine between b and k
kt, kr, etc are components of vectors; k0, k1 etc are of covectors
Compute magnitude of magnetic field before coordinate transformation
Reformulate the torus source
Put explicit casting in front of constants to avoid warnings
Convert gcov[][] and dxdxp[][] float to avoid warning; skip fields