
Copyright (C) 1997 D.H. Rothman and S. Zaleski. 

Covered by the GNU General Public License, see the file ../README.


             Table look up, 6 velocity code. 
                  Author S. Zaleski

This code is a simple FHP-I code. 

Changelog of small changes:
Olav van Genabeek march 1997:modified the forcing 
	                     and cast the table to 
	                     unsigned short

COMPILING THE CODE
Copy the contents of directory to 
your local directory.
Before you can compile you have to make sure that
you have an CC compiler.


RUNNING THE CODE 
To run the code:
Edit the parameters
in the file parameters.h
Specifically, edit the following lines

#define NX 16
#define NY 20
#define TPRINT 1000
#define TMAX 7000
#define DENSITY 0.5
#define FORCING_RATE .005 
#define TAVG 2000


Make sure NX and NY are the same as in Solid. 
The code will output results every TPRINT. 
Right now it will only 
output x and y-momentum and mass, but you can change this if
you wish. DENSITY is the initial density per link
and must have a value between 0 an 1. 
FORCING_RATE is the forcing in the x direction. 
TMAX is the maximum number of experiments before averaging. TAVG
is number of averages.

The following lines in the parameters file you can use for the
channel flow 

#define PARABOLAS 1
#if PARABOLAS
#define WALL_THICKNESS 1
#endif

When PARABOLAS is to to 1, the code will simulate flow
through a channel at y=1 and y=NY.

Compile it using 

make fhp6_simp

and run it using 

 fhp6_simp


- Changes of some of the parameters 
  in parameters.h will only become 
  effective after recompiling that is, type:
  make ( if does not work do first a " make clean").
  	

- This program will output two momentum files: x_mom and y_mom.
  These files contains values for the x-momentum times 2
  and y-momentum times 2/sqrt(3) (see file macros_6.h), summed
  over TAVG time steps.
  Then by including the mass at each 
  point (listed in file "mass")
  the velocity can be computed.
- make sure not to take the forcing rate too large. 


You may want to specify the location of solid grains. You may
do this using the program Solid in directory "processing". 
Read the README file in  directory "processing" for more information. 

SOLID_SITES 

If you do not wish solid sites, it may accelerate the code
then comment out the line:

#define SOLID_SITES












