
Thank you, thank you of reading me again.
Firstly and regrettely, there is a bug in the problem set.
The code we supply is for a 2-D hexogonal model.

Secondly, ilb2 need the following files

binaryfluid.c     initlat.c         model.c           testcons.c
initrock.c        ran1.c            util.c
gravity.c         mech_surf_tens.c  stream.c          writevel.c
collision.c       initlat2.c        user.c

boltzmann.h

plus some script files
bubble.scrpt mech.scrpt rayleigh.scrpt
----------------------------------------------------------------

On the sun workstations of 
athena please type:
"make CC=suncc"
to link and create the "ilb2" executabled.
For any other workstations use the apropiate
C compiler.
If errors occur during compiling please
let me know.

For your convenience, several script files are included:
	mech.scrpt
	bubble.scrpt
	rayleigh.scrpt
The first two scripts measure the surface tension.
The third script simulates the Rayleigh-Taylor unstability.

As usual you can run ilb2 by supplying it with an script file:

ilb2< <scriptfile>

The format of each script file is:
- each line in the script file 
   schould begin with a prefix listed below, followed 
   by a single space.
- prefixs are:
    '.':           indicates the end of the scrpt file
    'u':           set universe: read NX, NY
    'f':           output file name
    'n':           seed random generator
    's':           read solid and set universe
    'i':           initial flow param:1) prob_of_red,
                                      2) density
                                      3) density fluct. < 1.0
                                      4) vx
                                      5) vy
    'p':           set parameters : LAMBDA, perturb. of 
                                            surf. ten.
    'g':           set gravity 1) gx
                               2) gy
    'e':           experiment + 1) followed by a char
                                   which can be: 
                                   m: mechanical surface 
	                                 tension
                                   b: bubble
                                   r: Rayleigh Taylor 
	                              instability
                                   s: single phase flow 
	                              thru porous media
                                2) # of frames
                                3) # of steps between each frame


Please find supplied include the following scriptfiles:
bubble.scrpt        -compute the surface tension of a bubble
mech.scrpt          -compute the mech surface tension of layer















