init_geo.c

一个关于LB的经典程序。希望大家有好的东东就那出来

#include "global.h"
void init_geo()
{
int x,y;/* Driven Cavity */
/*   *BC* y=max_y-1  |-------- *|***| solid ** solid *|*** y=0 |*********|solid*|______________________*||* x=0  x=max_x-1     */
for(y=0;y<max_y-1;y++) 
wall[y][0]=wall[y][max_x-1]=WALL;	
for(x=0;x<max_x;x++)
wall[0][x]=WALL;
for(x=0;x<max_x;x++) 
wall[max_y-1][x]=SET_U;/* just for fun: inner obstacle */
/*    * for(y=max_y/4;y<3*max_y/4;y++) * 	
{* for(x=max_x/4;x<3*max_x/4;x++) wall[y][x]=WALL;* }*/
}