cd RegionFunctionality
dmlist time_cube.fits cols | grep "Axis#" --after-context=4

ds9 time_cube.fits -scale log &

xpaset -p ds9 regions command "circle(3811,3456,100)"

xpaset -p ds9 cube 80

import -window "SAOImage ds9" ds9_04_regions_07.png 
display < ds9_04_regions_07.png

# Double click region, Plot 3D
# Close region info box
sleep 5

import -window "Circle" ds9_04_regions_08.png 
display < ds9_04_regions_08.png

xpaset -p ds9 quit