/raid/lucy/Density/Cluster1/DW
/raid/lucy/Density/Cluster1/movies has XY, XZ, and YZ plots, ICL and derivative plots
note: the position plots have not been re-run with the new density definition
/raid/lucy/Density/Cluster1/Old has cid files using the old 1.096E-5 density cutoff
/raid/lucy/Density/Cluster1/Streams has cid files for streams:
Galaxy 59, stream 2, defined at snapshot 158
Galaxy 84, stream 1, defined at snapshot 76

c1_???.pden files link to /raid/craig/data/cluster1/pden/c1_???.pden

all_c1.sm -- similar and better formatted plots are made by
     allicl_plots.sm in /raid/lucy/Density
     # runall - runs clumping code on all snapshots

     # icl_all creates c1_iclfrac.list file with snapshot, numtot, numicl, num in streams>500

     # plot - plots the number of particles in various clump sizes (counted cumulatively)
     # output: c1_iclnum_all.gif

     # plotdiff - plots number of particles in clumps of various sizes,
     # counted separately
     # output: c1_iclnum_diff.gif

     # plotfrac - plots the fraction of ICL in clumps > 500
     # output: c1_iclfrac.gif

     # plotnum - plots the number of ICL particles in clumps > 500
     # output: c1_iclnum.gif

movie.sm
     # for each snapshot, movie outputs:
     #   -the whole cluster: c1_xy_000.gif, c1_xz_000.gif, c1_yz_000.g
     # all plots are in ./movies/ directory

plots.sm
     # plots creates ICL plots for each snapshot, with a dot indicating the current snapshot
     # output: movies/c1_iclnum$(snapshot).gif
     # blue - total ICL fraction in cluster
     # red - fraction of ICL in streams with > 500 particles

     # plotd makes plots of the derivative of ICL fraction and stream fraction 
     # as a function of time, with dots indicating the current snapshot
     # output: movies/c1_der$(snapshot).gif
     # note: limits are the same for all three clusters, so they seem large

stream.sm
     # runstr snapshot galaxy stream
     # runstr 76 84 1
     # runstr 158 60 2
     # runs clump code on given stream
     # output: c1_g$2_str$3_$(snapshot).cid(.stats and .lum)
     # all output files moved to directory ./Streams

     # I belatedly realized that streams defined at z=0 are not interesting, 
     # so subsequent analysis was only run on those defined earlier
     # calc creates .list file with 
     # num numtot numicl num(streams<500) num(streams>500) num(allstreams) num(biggest stream)
     # calc snapshot galaxy stream
     # calc 76 84 1