有两个错误，都在第16行
提示parse error和invalid lvalue in assignment
#include "udf.h"
#define rho 1.225
#define pi 3.14159
#define omega 1166.924043
#define theta0 11
#define dr 0.034275
#define dophy 9.0
#define r0 0.6855
#define c 0.1
#define h 0.025
DEFINE_SOURCE(xmom_source,cell,thread,dS,eqn)
{
real x[ND_ND], x1,y1,z1,r,v1,v2,v3,vx,vy,vz,v;
real phi,theta,alpha,ophy,cy,cx,l,d,Fn,Ft;
real Fx,Fy,Fz,sx,sy,sz,source;
begin_c_loop(c,c_thread)
{
C_CENTROID(x,cell,thread);
x1 = x[0];
y1 = x[1];
r=sqrt(x1*x1+y1*y1);
if (x1>0.&&y1>0.)ophy=asin(fabs(y1/r));
else if (x1<0.&&y1>0.)ophy=asin(fabs(x1/r))+pi/2.;
else if (x1<0.&&y1<0.)ophy=asin(fabs(y1/r))+pi;
else if (x1>0.&&y1<0.)ophy=asin(fabs(x1/r))+3.*pi/2.;
else if (x1==0.&&y1<0.)ophy=3./2.*pi;
else if (x1==0.&&y1>0.)ophy=0.5*pi;
else if (x1>0.&&y1==0.)ophy=0.;
else if (x1<0.&&y1==0.)ophy=pi;
v1=C_U(cell,thread);
v2=C_V(cell,thread);
v3=C_W(cell,thread);
vx=v1+omega*r;
vy=v2-omega*r;
vz=v3;
v=sqrt(vx*vx+vy*vy+vz*vz);
phi=atan(-vz/(sqrt(vx*vx+vy*vy)));
alpha=theta0-phi;
cy=0.11*alpha*180.0/pi;
cx=0.001*alpha*180.0/pi;
l=0.5*rho*v*v*cy*c*dr;
d=0.5*rho*v*v*cx*c*dr;
Fn=l*cos(phi)-d*sin(phi);
Ft=l*sin(phi)-d*cos(phi);
Fx=Ft*sin(ophy);
Fy=-Ft*cos(ophy);
Fz=Fn;
sx=-2*dophy*Fx/(2*pi);
sy=-2*dophy*Fy/(2*pi);
sz=-2*dophy*Fz/(2*pi);
source=sx;
}
end_c_loop(c,c_thread)
dS[eqn]=-dophy*rho*c*dr*(2*vx+vy*vy+vz*vz)*(cy*sin(phi)-cx*cos(phi))*sin(ophy)/(2*pi);
return source;
}
#include <stdio.h>
int jc(int m){
if(m!=1) return m*jc(m-1);
else return 1;
}
int c(int m,int n){
if(m>=n) return jc(m)/(jc(n)*jc(m-n));
}
int main(void){
int m,n;
scanf("%d%d",&m,&n);
printf("%d\n",c(m,n));
return 0;
}本回答被网友采纳
