请帮忙看一下错误之处，这是为了得到处so2的份额

yuanwanli

#include "udf.h"

#define M 1.4e-3
#define N 0.5
#define A 1.1e-4
#define B 1.8369

DEFINE_PROFILE(Fraction_profile, cell, thread)
{
real s1;
real fraction=C_T(cell,thread);
real s1 = species_mf[3];
if (s1>5e-5)
fraction = M*(M-s1)/(N*(M+s1)-A*pow((M-s1),B));
else
fraction = 0.0;
return fraction;
}

gearboy78

为啥不用C_YI宏？

yuanwanli

什么意思？denfine C_YI？可应用于边界条件吗？我查一下