求助 程序出错,实在不知道什么原因？源程序以贴出

hdlishuhua

#include "udf.h"
#include "sg.h"
DEFINE_SOURCE(source,c,t,dS,eqn)
{real NA;
real P_SO2;
real C_SO2;
real kg;//
real kL;//
real H_SO2;  //mol/m3*pa//
real E;
//  液相传值系数的确定//
real  d; //m//
real  D_SO2;
real  Sc;
real  nu_g;   //Pa*s//
real  ru_g;                 //g/mol//
//real  Re_d;  //
real M_air;                //g/mol//
real M_SO2;                //g/mol//
real V_air;               //g/mol//
real V_SO2;              //g/mol//
real Re_d;
//气相传值系数的确定//
real b;   //工作温度下水的动力粘度//
real a;    //气体的分子体积  cm3/mol//
real DL;
real tp;
real pi;
real A;
real source;
//数据//
P_SO2=C_P(c,t);
C_SO2=C_YI(c,t);
H_SO2=5.541e-3;  //mol/m3*pa//33//
E=0.0001;
NA=(H_SO2*P_SO2-C_SO2)/(1/kg+H_SO2/E/kL);//35//
//  液相传值系数的确定//
//d=0.001；   //m//
Sc=nu_g/ru_g/D_SO2;
nu_g=1.92e-5;   //Pa*s//
ru_g=1.11;                 //g/mol//
real  Re_d;
M_air=29.9;                //g/mol//42//
M_SO2=64;                //g/mol//
V_air =29.9;               //g/mol//
V_SO2=44.8;              //g/mol//
//Re_d=d*fabs(C_U(c,t))*ru_g/C_U(c,t);//
D_SO2=0.00986*C_T(c,t)*sqrt(1/M_air+1/M_SO2)/P_SO2/((V_air^(1.0/3)+V_SO2^(1.0/3))^2);//48//
//kg=(2+0.55*sqrt(Re_d)*(Sc^(1/3)))*D_SO2/d;//
kg=(2+0.55*sqrt(Re_d)*pow(Sc,(1.0/3)))*D_SO2/d;
//气相传值系数的确定//
b=1.005*10^(-3);   //工作温度下水的动力粘度//
a=22400;    //气体的分子体积  cm3/mol//
DL=14*0.0000000001/((a^0.6)*(b^1.1));
tp=3.0;
pi=3.14;
kL=sqrt(4*DL/pi/tp);
//源方程//
A=3.0/2.0/d;
source=-NA*A;
//参数变量//

//source=-NA*A;//
      dS[eqn]=0;
      return source;
   }   
报错

hdlishuhua

cpp -IC:\Fluent.Inc\fluent6.2.16/src -IC:\Fluent.Inc\fluent6.2.16/cortex/src -IC:\Fluent.Inc\fluent6.2.16/client/src -IC:\Fluent.Inc\fluent6.2.16/multiport/src -I. -DUDFCONFIG_H="<udfconfig.h>" D:\fluent.inc\modle3.c.txt
source definition shadows previous definition
Error: D:\fluent.inc\modle3.c.txt: line 33: parse error.
Error: D:\fluent.inc\modle3.c.txt: line 42: parse error.
Error: D:\fluent.inc\modle3.c.txt: line 48: invalid type for integral binary expression: float ^ float.

hdlishuhua

怎么没有高手指导一下？

tday

不是有提示行么？按提示的那几行去检查下，应该能查出来的吧

mingming999

这个是模拟脱硫的吧？