############################################################
# Building and running Diffusion Limited Aggregation
# program 'dla' on the SGIs.
############################################################
einstein% pwd; ls
/usr2/people/phy329/rand/ex2
Makefile   dla.f

einstein% make
        f77 -g -c dla.f
        f77 -g -L/usr/local/lib dla.o -lp329f -o dla

einstein% dla
 usage: dla <size> <npart> [<r0> <bias>]

        Current default <r0>:          0.75
        Current default <bias>:  0.0000E+00

        Program reads initial fixed-particle coordinates 
        (integers x,y; 1 <= x,y <= size)  from standard
        input, writes final fixed positions to standard
        output.

############################################################
# "Interactive" run: 10 particles on a 100 x 100 arena
############################################################
einstein% dla 100 10
 dla: size            100
 dla: npart            10
 dla: r0       0.7500000000000000    
 dla: bias     0.0000000000000000E+00

############################################################
# Put a single fixed cell at the center of the arena
############################################################
50 50
^D
 dla: Read            1 particle positions.
 dla: Particle            1 fixed after         3192 steps
 dla: Particle            2 fixed after         7100 steps
 dla: Particle            3 fixed after        11890 steps
 dla: Particle            4 fixed after         6564 steps
 dla: Particle            5 fixed after        18517 steps
 dla: Particle            6 fixed after        15199 steps
 dla: Particle            7 fixed after         1340 steps
 dla: Particle            8 fixed after         3049 steps
 dla: Particle            9 fixed after        19606 steps
 dla: Particle           10 fixed after         9202 steps
           51          48
           51          49
           55          49
           50          50
           51          50
           53          50
           54          50
           48          51
           49          51
           52          51
           52          52
 dla: Wrote           11 particle positions.

############################################################
# The same calculation but with initial positions read 
# from file 'state0' and final positions written to 
# 'state10'.
############################################################
einstein% more state0
50 50
einstein% dla 100 10 < state0 > state10
 dla: size            100
 dla: npart            10
 dla: r0       0.7500000000000000    
 dla: bias     0.0000000000000000E+00
 dla: Read            1 particle positions.
 dla: Particle            1 fixed after         3192 steps
 dla: Particle            2 fixed after         7100 steps
 dla: Particle            3 fixed after        11890 steps
 dla: Particle            4 fixed after         6564 steps
 dla: Particle            5 fixed after        18517 steps
 dla: Particle            6 fixed after        15199 steps
 dla: Particle            7 fixed after         1340 steps
 dla: Particle            8 fixed after         3049 steps
 dla: Particle            9 fixed after        19606 steps
 dla: Particle           10 fixed after         9202 steps
 dla: Wrote           11 particle positions.

############################################################
# Add ten more particles to the simulation and save 
# final results in file 'state20'.
############################################################
einstein% dla 100 10 < state10 > state20
 dla: size            100
 dla: npart            10
 dla: r0       0.7500000000000000    
 dla: bias     0.0000000000000000E+00
 dla: Read           11 particle positions.
 dla: Particle            1 fixed after         3129 steps
 dla: Particle            2 fixed after         7146 steps
 dla: Particle            3 fixed after         5198 steps
 dla: Particle            4 fixed after         7751 steps
 dla: Particle            5 fixed after         7628 steps
 dla: Particle            6 fixed after        12129 steps
 dla: Particle            7 fixed after          735 steps
 dla: Particle            8 fixed after         4638 steps
 dla: Particle            9 fixed after         1594 steps
 dla: Particle           10 fixed after         1611 steps
 dla: Wrote           21 particle positions.