/raid/lucy/Density/Cluster2/DW
/raid/lucy/Density/Cluster2/Old has cid files using the old 1.096E-5 density cutoff
/raid/lucy/Density/Cluster2/Streams has cid files for streams:
Galaxy 19, stream 2, snapshot 158
Galaxy 48, stream 1, snapshot 78
Galaxy 48, stream 2, snapshot 78
Galaxy 57, stream 1, snapshot 158
Galaxy 57, stream 2, snapshot 158

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

all_c2.sm
     # runall - runs clumping code on all snapshots
     # icl_all prints file with:
     # snapshot, numtot, numicl, num(streams>500)
     # output: c2_iclfrac.list

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

     # plotfracgals - plots the fraction of ICL in clumps > 500
     # for each of our test galaxies
     # output: c2_iclfrac_g$g.gif

     # plots the number of ICL particles in clumps > 500
     # for each test galaxy
     # output:  c2_iclnum_g$g.gif

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

     # icl_gal makes file with the following for each snapshot, for galaxies 19 and 48:
     # numtot(in galaxy), numicl(from galaxy), num(streams>500 from galaxy)
     # output: c2_gals_iclfrac.list

all_c2_old.sm
     # runall - runs clumping code on all snapshots
     # icl_all - prints files with number of ICL particles in clumps of
     # various sizes

     # c2_icl.list prints the number of particles in clumps of the
     # given size or smaller (cumulative)

     # c2_icl_diff.list prints the number of particles in clumps
     # in the given size range

     # c2_iclfrac.list prints the number of particles in clumps
     # >500, for the cluster and each test galaxy

     # plot - plots the number of particles in various clump sizes (counted cumulatively)

     # plotdiff - plots number of particles in clumps of various sizes,
     # counted separately

     # plotfrac - plots the fraction of ICL in clumps > 500

     # plotfracgals - plots the fraction of ICL in clumps > 500
     # for each of our test galaxies

     # plots the number of ICL particles in clumps > 500
     # for each test galaxy

     # plotnum - plots the number of ICL particles in clumps > 500

cluster2.sm # mkfiles makes ascii position and velocity files for the given snapshots # plot makes graph with XY plots at 8 snapshots # limits: -3000 to 3000 # output: all_3000.eps and all_3000.gif # plot5 makes graph with XY plots at 8 snapshots # limits: -5000 to 5000 # output: all_5000.eps and all_5000.gif

header.sm
     # reads snapshot headers and creates file with higher-precision redshifts and expansion factors

movie.sm
     # for each snapshot, movie outputs:
     #   -for the whole cluster: c2_xy_000.gif, c2_xz_000.gif, c2_yz_000.gif
     #   -for each galaxy: c2_xy_000.g??.gif, c2_xz_000.g??.gif, c2_cl_000.g??.gif

plots.sm
     # plots makes ICL plots for the cluster for each snapshot
     # output movies/c2_iclnum$(snapshot).gif

     # plots_der makes plots of the derivatives fo ICL fraction and stream fraction
     # output: movies/c2_der$(snapshot).gif

plots_gals.sm
     # plots makes ICL plots for each snapshot, for each galaxy
     # output: movies/c2_iclnum_g$(gal)_$(snapshot).gif

stream.sm
     # runstr snapshot galaxy stream
     # runs clumping code on given stream
     # run with:
     # runstr 158 57 1
     # runstr 158 57 2
     # runstr 158 19 2
     # runstr 78 48 1
     # runstr 78 48 2
     # output: c2_g$(gal)_str$(stream)_$(snapshot).cid
     # all output files moved to ./Streams

     # run1 snapshot_to_run
     # calc snapshot galaxy stream
     # calc creates file with: 
     # snapshot, numtot, numicl, num(streams<500), num(streams>500), num(all streams)
     # num(biggest stream)
     # run with:
     # calc 158 57 1
     # calc 158 57 2
     # calc 158 19 2
     # calc 78 48 1
     # calc 78 48 2
     # output: c2_g$(gal)_str$(stream)$(snapshot)_icl.list