大涡模拟程序代码

大涡模拟的代码程序，对自己编程的有很好的指导作用

/* duct.c *************************************\
*  - Turbulent Flow in Duct modeled by LES -   *
*        two-stage time differencing           *
*  written by Andrei Chernousov Mar 16 2001   *
*  E-mail:                 *
*  http://www.geocities.com/andrei_chernousov  *
\**********************************************/

#include 
#include 
#include  

#include 

/* continuation flag （to be changed by SIGINT handler） */
int continFlag = 1;

/* SIGINT handler */
#pragma argsused
void termHandler（int sigType）
{
  sigType = sigType;
  if （continFlag--）; /* 1 -> 0 */
  fprintf（stderr “caught [Ctrl+C] terminating...“）;
}

/*--- Combinations of the ratio of specific heats ---*/
#define K        1.4
#define K_1      0.4
#define _K_1     （-0.4）
#define _05K     0.7
#define K_K_1    3.5
#define K_1_2    0.2
#define _1_2_K_1 1.25

#define Cv      717.75  /* specific constant volume heat capacity J/（kg*K） */
#define R_VOZD  287.10  /* gas constant for air */
#define KR_VOZD 401.94  /* R_VOZD * K */

#define Fx1 xU1 /* aliases for x-fluxes */
#define Fx2 xU2
#define Fx3 xU3
#define Fx4 xU4
#define Fx5 xU5

#define Fy1 yU1 /* aliases for y-fluxes */
#define Fy2 yU2
#define Fy3 yU3
#define Fy4 yU4
#define Fy5 yU5

#define Fz1 zU1 /* aliases for z-fluxes */
#define Fz2 zU2
#define Fz3 zU3
#define Fz4 zU4
#define Fz5 zU5

/*=== Global variables ===*/
float
  ***U1_ ***U2_ ***U3_ ***U4_ ***U5_ /* substitutional pointers */
  /* Cons. variables at the cells‘ centroids */
  ***U1  ***U2  ***U3  ***U4  ***U5 
  ***U1p ***U2p ***U3p ***U4p ***U5p
  /* Cons. variables at cell boundaries in x-direction */
  ***xU1 ***xU2 ***xU3 ***xU4 ***xU5 
  ***U1x ***U2x ***U3x ***U4x ***U5x
  /* Cons. variables at cell boundaries in y-direction */
  ***yU1 ***yU2 ***yU3 ***yU4 ***yU5 
  ***U1y ***U2y ***U3y ***U4y ***U5y
  /* Cons. variables at cell boundaries in z-direction */
  ***zU1 ***zU2 ***zU3 ***zU4 ***zU5 
  ***U1z ***U2z ***U3z ***U4z ***U5z
   
  R0 /* initial density */
  R U V W P C       /* vector of primitive flow parameters */
  u1 u2 u3 u4 u5     /* conservative flow variables         */
  deltaX deltaY deltaZ /* cell spacings                       */
  deltaT                 /* step in time                        */
  deltaT_X deltaT_Y deltaT_Z /* convenient ratios             */
  _deltaX _deltaY _deltaZ
  _4deltaX _4deltaY _4deltaZ
  _2deltaX _2deltaY _2deltaZ;

char 
  Stage /* indicator of the current stage */
  str[80]; /* buffer string */

unsigned
  LEN HIG DEP /* numbers o cells in x- y- z-directions for 3d box */
  LENN HIGG DEPP
  Answer  /* flag of continuing solution */
  i  j  k i__  j__  k__ /* indices */
  l
   
  mem     /* total amount of dynamic memory required （in bytes） */
  step    /* current time step */
  f_step  /* frame taking step */
  numstep; /* overall prescribed number of time steps

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件      66752  2009-08-06 23:20  les_code\example\duct\100

     文件      71110  2009-08-06 23:20  les_code\example\duct\1000

     文件      71310  2009-08-06 23:20  les_code\example\duct\1100

     文件      64865  2009-08-06 23:20  les_code\example\duct\12

     文件      71485  2009-08-06 23:20  les_code\example\duct\1200

     文件      71601  2009-08-06 23:20  les_code\example\duct\1300

     文件      71675  2009-08-06 23:20  les_code\example\duct\1400

     文件      71779  2009-08-06 23:20  les_code\example\duct\1500

     文件      71788  2009-08-06 23:20  les_code\example\duct\1600

     文件      71813  2009-08-06 23:20  les_code\example\duct\1700

     文件      71894  2009-08-06 23:20  les_code\example\duct\1800

     文件      71937  2009-08-06 23:20  les_code\example\duct\1900

     文件      67592  2009-08-06 23:20  les_code\example\duct\200

     文件      71979  2009-08-06 23:20  les_code\example\duct\2000

     文件      71955  2009-08-06 23:20  les_code\example\duct\2100

     文件      71956  2009-08-06 23:20  les_code\example\duct\2200

     文件      72009  2009-08-06 23:20  les_code\example\duct\2300

     文件      72002  2009-08-06 23:20  les_code\example\duct\2400

     文件      72119  2009-08-06 23:20  les_code\example\duct\2500

     文件      72182  2009-08-06 23:20  les_code\example\duct\2600

     文件      72224  2009-08-06 23:20  les_code\example\duct\2700

     文件      72332  2009-08-06 23:20  les_code\example\duct\2800

     文件      72185  2009-08-06 23:20  les_code\example\duct\2900

     文件      68570  2009-08-06 23:20  les_code\example\duct\300

     文件      72277  2009-08-06 23:20  les_code\example\duct\3000

     文件      72185  2009-08-06 23:20  les_code\example\duct\3100

     文件      72273  2009-08-06 23:20  les_code\example\duct\3200

     文件      72302  2009-08-06 23:20  les_code\example\duct\3300

     文件      72359  2009-08-06 23:20  les_code\example\duct\3400

     文件      72348  2009-08-06 23:20  les_code\example\duct\3500

............此处省略100个文件信息