#//u(i,j)        x方向的速度u;或者θ方向速度uθ
#//u(i,j)        y方向的速度v;或者径向速度ur'h\K
#//pc(i,j)       压力修正 p'[OS
#//p(i,j)        压力p-=5-+
#//p(i,j)        密度ρBk
#//p(i,j)        扩散系数γQ"G48E
#//t(i,j)        温度Tn7/bKr
#//ake(i,j)        湍流脉动能量k971E
#//dis(i,j)      动能的耗散率ε4
 
//amut(i,j)     湍动扩散系数+
//gen(i,j)      湍流能量的生成率S
//f(i,jnf)     不同的φ变量M
//lsolve(nf)   ＝1时,求解变量f(i,j,nf)  8
//lprint(nf)   ＝1时,打印变量f(i,j,nf)   ~f
//lblk(nf)     ＝1时,对变量f(i,j,nf)应用块修正b5|{@U
//mode         选择坐标系的变量. 目前只支持mode=1.
mode=1 直角坐标系(x~y)BMI
mode=2 圆柱坐标系(r~z)^pg>)
mode=3 极坐标(r~θ)   Z^r[LE
在solve子程序中求解变量f(i,j,nf)的重复扫描次数 b
变量f（i,j,nf）的字符性标题 |L;
 
//xl    计算区域在x方向上的宽度)?
//yl    计算区域在y方向上的宽度@ob
//l1     x方向上主控制体的网格数。也是x方向上压力节点位置的最后一个i值。
//m1   y方向上主控制体的网格数。也是y方向上压力节点位置的最后一个j值。
//dt           时间步长 △t
//第三部分网格设定的变量
//x(i)        节点位置的x值
//xu(i)       主控制容积在x方向的界面位置，即速度u(i,j)所在位置
//xdif(i)     差值 x(i)-x(i-1)
 //xcv(i)      主控制容积在x方向上的宽度
//xcvs(i)     速度u(i,j)的控制容积在x方向上的宽度
//y(j)        节点位置的y值
//yv(j)       主控制容积在y方向的界面位置，即速度v(i,j)所在位置
//ydif(j)     差值 y(j)-x(j-1)
//ycv(j)      主控制容积在y方向上的宽度
//ycvs(j)     速度v(i,j)的控制容积在y方向上的宽度
////r(i,j)      主网络节点的半径r
////rmn(j)      在速度v(i,j)所在处的半径r之值
////sx (j)      主网格节点位置y(j)处x方向上的标尺因子
//sxmn(j)     在界面位置yv(j)上x方向的尺度因子
       //上面四个变量应用于非直角坐标系的情况
//xcvi(i,j)   xcv(i) 中与u(i,j)的控制容积相覆盖的部分
//xcvi(i,j)   xcv(i) 中与u(i+1,j)的控制容积相覆盖的部分
//ycvr(j)     主控制容积垂直于x方向的面的面积
//ycvrs(j)    速度v(i,j)的控制容积垂直于x方向的面的面积
//arx(j)      与x方向相垂直的控制容积的面积
//arxj(j)     arx(j)中与速度v(i,j)的控制容积相覆盖的部分
//arxjp(j)    arx(j)中与速度v(i,j+1)的控制容积相覆盖的部分
//arxj与arjxp实际上对应于x方向上的ycvi与ycvip
 
//第四部分差分方程系数设定的变量
//con(i,j)    离散方程中的常数项b，在子程序gamsor中又作为存储sc的单元
//aip(i,j)    系数ae
//aim(i,j)    系数aw
//ajp(i,j)    系数an
//ajm(i,j)    系数as
//ap(i,j)     系数ap,在在子程序gamsor中又作为存储sp的单元
//flow         穿过控制容积界面的质量流率
//diff         扩散传导性 D
//acof         由DIFLOW子程序计算的量，它给出了对流与扩散作用的联合影响
 
//第五部分求解差分方程过程中的变量求解
//du(i,j)     影响u(i,j)的de
//dv(i,j)     影响v(i,j)的dn
//pv(j)       用于计算主网络节点i,j上的质量流率的ρvr的插值因子：计算式如下： //fvp(j)      说明同上
//fx(i)       用于计算主控制容积界面(即速度u(i,j)所在处)的密度
 
//            rhom的插值因子，计算式如下：
 
//fxm(i)      说明同上
//pt(i)或pt(j) tdma中的转换系数(消元过程中)l
//qt(i)或qt(j) tdma中的转换系数(消元过程中)
//第六部分 index变量
int nf;               //nf           标明不同φ变量的下标值
int nfmax;         //nfmax        设有存储单元的nf的最大值
int np;              //(nfmax)      p(i,j)实际为f(i,j,nfmax)
int nrho;           //(nfmax+1)    rho(i,j)实际为f(i,j,nfmax)
int ngam;         //(nfmax+2)    gam(i,j)实际为f(i,j,nfmax+2)n
int l2;              //l2           (l1-1)
int l3;              //l3           (l1-2)
int m2;           //m2           (m1-1)
int m3;           //m3           (m1-2) 2
int ist;                    //ist          i的第一个内节点值`
int jst;                    //jst          j的第一个内节点值
int iter;                   //iter         不稳态问题的步进计数
int last;                   //last         用户所规定的最大迭代计数
int iter1;             //iter1        一个时间点求解setup2( )的迭代次数
double time;              //time         不稳态问题中的时间t
int ipref;              //ipref        压力参考结点的i值
int jpref;                //jpref        压力参考结点的j值=/
 
//第七部分其它变量
double rhocon;           //rhocon       密度为常数的问题中的ρ
int lstop;                 //lstop        =1时，停止计算
double  smax;            //smax         p'方程中的“质源”的最大值
double  ssum;            //ssum         p'方程中的“质源”的代数和