UDF中UDS调用问题

wqlian1990

看过一篇关于水气变化反应的文献上面关于反应速率的UDF。其编写如下：
real mole_total,Keq;
real CO_mole_frac,H2_mole_frac,C02_mole_frac,H20_mole_frac;
real P_C0,P_H2,P_C02,P_H20,P_t0tal;
mole_total=C_YI(c,t,0)/wk_0+C_YI(c,t,l)/wk_l+C_YI(c,t,2)/wk_2+C_YI(c,t,3)/wk_3;
CO_mole_frac=C_YI(c,t,0)/(wk_0*mole_total);
H2_mole_frac=C_YI(c,t,l)/(wk_l*mole_total);
C02_mole_frac=C_YI(c,t,2)/(wk_2*mole_total);
H20_mole_frac=l-CO_mole_frac-H2_mole_frac-C02_mole_frac;
P_total=AB S_P(C_P (c ,t) ,op_pres);
PC 0=P_total * C O mole_frac ;
P_H2=P_total * H2_mole_frac ;
P_C02=P—total*C02—mole—frac;
P_H20=P_t0tal-P_C0-P_H2-P_C02;
Keq=exp(5.0992511+4764.1 l/C_T(c,t)+5.6273903*0.001*C_T(c,t)-2.189647*0.000001*pow(C_T(c,t),2)+4.17603*0.0000000001*pow(C_T(c,t),3)-1.9439412*log(C_T(c,t)));
C_UDSI(c,t,CO_mole)=1.74*pow(10,13)*0.512*pow(C_T(c,t),-8.5)*exp(-35/(8.314*C_T(c,t)))*(P_C0*P_H20-P_C02*P_H2/Keq);

关于这篇文献上，反应速率为何不用宏DEFINE_VR_RATE(user_rate,c,t,r,mole_weight,species_mf,rate,rr_t)引出。文献提及采用了两个UDS，但该程序C_UDSI(c,t,CO_mole)是何意？CO_mole指的又是什么呢？上面C_UDSI(c,t,CO_mole)=1.74*pow(10,13)*0.512*pow(C_T(c,t),-8.5)*exp(-35/(8.314*C_T(c,t)))*(P_C0*P_H20-P_C02*P_H2/Keq);语句是给它赋值吗？