Question

fluent 编译udf错误parse error 请高手帮帮忙啊

Error: E:\shisuan\source.c: line 23: parse error.
Error: E:\shisuan\source.c: line 26: parse error.
Error: E:\shisuan\source.c: line 28: parse error.
Error: E:\shisuan\source.c: line 33: parse error.
这是原程序#include "udf.h"
#define C_DPMS_K(c,t)C_STORAGE_R(c,t,SV_K)
#define C_DPMS_D(c,t)C_STORAGE_R(c,t,SV_D)
#define P_POS(p)
real C1 = 17.5; /* main constant */
real height = 0.42; /* vessel height */
real bottom_fraction = 0.2; /* fraction below which k is damped */
real top_fraction = 0.85; /* fraction above which e is augmented */
real bottom;
real top;
real Ce = 1.29; /* epsilon constant */
real Ct1 = 4; /* top epsilon constant */
real rho_U = 980; /* upper density limit */
real rho_L = 3; /* lower density limit */
real Vslip; /* slip velocity */
DEFINE_DPM_SOURCE(dpm_source,c,t,S,strength,p)
{
bottom = bottom_fraction*height;
top = top_fraction*height;
Vslip = P_VEL(p)[1]-C_V(c,t);
if (C_R(c,t)>=rho_L)
{
if (P_POS(p)[1]<bottom)
{
C_DPMS_K(c,t) += P_POS(p)[1]/bottom*C1*9.81*C_R(c,t)*P_MASS(p)/P_RHO(p)*Vslip*strength*CURRENT_TIMESTEP;
C_DPMS_D(c,t) += P_POS(p)[1]/bottom*Ce*C1*9.81*C_R(c,t)*P_MASS(p)/P_RHO(p)*Vslip*strength*CURRENT_TIMESTEP*C_(c,t)/C_K(c,t);
}
if (P_POS(p)[1]>bottom)
{
C_DPMS_K(c,t) += C1*9.81*C_R(c,t)*P_MASS(p)/P_RHO(p)*Vslip*strength*CURRENT_TIMESTEP;
C_DPMS_D(c,t) += Ce*C1*9.81*C_R(c,t)*P_MASS(p)/P_RHO(p)*Vslip*strength*CURRENT_TIMESTEP*C_D(c,t)/C_K(c,t);
}
if (P_POS(p)[1]>=top)
{
C_DPMS_D(c,t) += Ct1*strength*(P_POS(p)[1]/height-top_fraction)*CURRENT_TIMESTEP;
}
}
return source;
}原程序

Answer 1

P_POS(p) 这个宏是什么，按照上面的代码展开自然是错误的