#//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'方程中的“质源”的代数和
Search - mode 3 S