Start rewriting all visualization/plotting scripts: a raw loader
Compose two images into one
Convert a loaded structure into flux (in Jy) or luminosity (in ergs/s)
Pass para to Data::output() for more flexible outputs; save nu0
Update the raw loader to take into account nu0
Initialize c.nu0[] because we will need to output it
Since we have access to the constants; we can easily skip unused nu0s
Make file matching easy
Skip the outermost cells
The integration is still accurate for a larger step cap
Further speed up the code by using larger step cap
Bug fix: missed a factor of 1/4pi in L_j_ff()
Add references to the emissivities
Setup and cleanup IDL environment for plotting to EPS files
Add the IDL scripts for interstellar scattering
Merge branch 'develop' into master