#include "global.h"
void init()
{
int x,y;/* initialise distributions for u=v=0 */	
  for(y=0;y<max_y;y++)	
   {
     for(x=0;x<max_x;x++)		
      {			
        f[0][y][x]=rho_0*4.0/9.0;			
        f[1][y][x]=f[2][y][x]=f[3][y][x]=f[4][y][x]=rho_0/9.0;
        f[5][y][x]=f[6][y][x]=f[7][y][x]=f[8][y][x]=rho_0/36.0;
       }	
    }
}