This patch does initial setup of an 1/r**2 density profile with source at the x=y=z=0 corner of the box, 
w. constant temperature of 100K throughout the box. This is specifically done to run test 6 in Ilievs 
RT comparison project paper.
