双椭球热源与高斯热源合并的udf_尴尬做一次伸手党谁有高斯热源模型的UDF编程的程序啊，跪求...

1，

# include"udf.h"

DEFINE_PROFILE(heat_flux,t,i)

{

real xx[ND_ND];

real a, pi, r, ef,p, heat, x,y;

face_t f;

ef= 0.8;

pi = 3.141596;

r = 0.003;

p=1400;

heat = 3 *ef*p / (pi*r*r);

begin_f_loop(f,t)

{

x = xx[0];

y = xx[1];

if (sqrt(x*x + y*y) <= 0.003)

F_PROFILE(f, t, i) = heat*exp(-3 * (x*x + y*y) / (r*r));

else F_PROFILE(f,t,i) = 0;

}

end_f_loop(f,t)

}

2，

# include "udf.h"

# define v 10

# define delt .05

# define delx .5

# define I0 2*10000000

# define n 39

DEFINE_PROFILE(g_profile,t,i)

{

int x[ND_ND];

int j,k;

double xreg;

double xl[39],I[39],y[3