dw.sm
     # streams snapshot galaxy clumpnum
     # runs with: snapshot 58, gal 48, stream 1; snapshot 90, gal 19, stream 1
     # prints a file with a histogram of dw values for each snapshot
     # dw's are calculated for the given clump from the given galaxy,
     # defined at the given snapshot.
     # output: c2_g$1_c$2_$(snapshot).dw

     # meds galaxy clumpnum
     # reads in the dw file for each snapshot, finds the median and
     # upper quartile - lower quartile, prints to file
     # output: med$(gal).list

     # medplot galaxy clump#
     # plots median and quartile values over time, making karaoke plots
     # output: c2_g$1_c$2_med_$(snapshot).gif

     # mkhist galaxy clump#
     # makes a histogram of the dw values for each snapshot
     # output: c2_g$1_c$2_hist_$(snapshot).gif



plot_stream.sm
     # plot_stream
     # plots cluster and marks particles in the given clump defined at given timestep
     # runs for:
     # galaxy 48, stream 1, snapshot 58
     # galaxy 19, stream 1, snapshot 90
     # output: c2_dw_xy_g$(gal)_c$(stream)_$(snapshot).gif
     #         c2_dw_xz_g$(gal)_c$(stream)_$(snapshot).gif