converting to fits

started the conversion to using fits files
done: starformation.py
todo: analysis.py and plot.py