#include?“udf.h“
#include?“math.h“
#define?T?10.0
#define?PIE?3.14159265
#define?miu?1.7894e-5
#define?density?1.225


DEFINE_SOURCE（xmom_sourcectdSeqn）
{
	real?x[ND_ND];
	real?sourcer;
	C_CENTROID（xct）;
	
	r=x[0]*x[0]+x[1]*x[1]+x[2]*x[2];/*?r=r*r?*/
	source=0;
	dS[eqn]=0;
	return?source;
}

DEFINE_SOURCE（ymom_sourcectdSeqn）
{
	real?x[ND_ND];
	real?sourcertime;
	C_CENTROID（xct）;
	
	r=x[0]*x[0]+x[1]*x[1]+x[2]*x[2];/*?r=r*r?*/
	time=0.001;
	if（r<