make CMDs
power series
emission profile weighted mean intensity
plot schechter function:
remove SDSS sources with photometric errors above some threshold
join individual sources and compile database of src with the same lines
generate velocity field spider diagrams