[Mascarius]

Hi, i have a problem with the syntax of the following code. But I do not know where the problem is. I will appreciate if you can help me with this.

#include "udf.h"
static const real Arrhenius = 1.4e08;
static const real E_Activation = 1.612e08;
#define SMALL_S 1.e-29
DEFINE_HET_RXN_RATE(arrh,c,t,hr,mw,yi,rr,rr_t)
{
Domain **domain_reactant = hr->domain_reactant;
real *stoich_reactant = hr->stoich_reactant;
Thread **pt=THREAD_SUB_THREADS(t);
/* for each phase */
Thread *tp = pt[0];
Thread *ts = pt[1];
int *reactant = hr->reactant;
int i;
int sp_id;
int dindex;
Thread *t_reactant;
real ci;
/* obtain temperature from cell for each phase */
real Tp=C_T(c,tp);
real Ts=C_T(c,ts);
real T=(Tp+Ts)/2.;
/* instead of compute rr directly, compute log(rr) and then
take exp */
*rr = 0;
for (i=0; i < hr->n_reactants; i++)
{
sp_id = reactant[i]; /* species ID to access mw and yi */
if (sp_id == -1) sp_id = 0; /* if phase does not have species,
mw, etc. will be stored at index 0 */
dindex = DOMAIN_INDEX(domain_reactant[i]);
/* domain index to access mw & yi */
t_reactant = THREAD_SUB_THREAD(t,dindex);
/* get conc. */
ci = yi[dindex][sp_id]*C_R(c,t_reactant)/mw[dindex][sp_id];
ci = MAX(ci,SMALL_S);
*rr += stoich_reactant[i]*log(ci);
}
*rr += log(Arrhenius + SMALL_S) -
E_Activation/(UNIVERSAL_GAS_CONSTANT*T);
/* 1.e-40 < rr < 1.e40 */
*rr = MAX(*rr,-40);
*rr = MIN(*rr,40);
*rr = exp(*rr);

[gearboy]

Try using VC++ Udf Studio. It can tell you where to revise and it's free for less than two macros.

Originally Posted by Mascarius

Hi, i have a problem with the syntax of the following code. But I do not know where the problem is. I will appreciate if you can help me with this.

#include "udf.h"
static const real Arrhenius = 1.4e08;
static const real E_Activation = 1.612e08;
#define SMALL_S 1.e-29
DEFINE_HET_RXN_RATE(arrh,c,t,hr,mw,yi,rr,rr_t)
{
Domain **domain_reactant = hr->domain_reactant;
real *stoich_reactant = hr->stoich_reactant;
Thread **pt=THREAD_SUB_THREADS(t);
/* for each phase */
Thread *tp = pt[0];
Thread *ts = pt[1];
int *reactant = hr->reactant;
int i;
int sp_id;
int dindex;
Thread *t_reactant;
real ci;
/* obtain temperature from cell for each phase */
real Tp=C_T(c,tp);
real Ts=C_T(c,ts);
real T=(Tp+Ts)/2.;
/* instead of compute rr directly, compute log(rr) and then
take exp */
*rr = 0;
for (i=0; i < hr->n_reactants; i++)
{
sp_id = reactant[i]; /* species ID to access mw and yi */
if (sp_id == -1) sp_id = 0; /* if phase does not have species,
mw, etc. will be stored at index 0 */
dindex = DOMAIN_INDEX(domain_reactant[i]);
/* domain index to access mw & yi */
t_reactant = THREAD_SUB_THREAD(t,dindex);
/* get conc. */
ci = yi[dindex][sp_id]*C_R(c,t_reactant)/mw[dindex][sp_id];
ci = MAX(ci,SMALL_S);
*rr += stoich_reactant[i]*log(ci);
}
*rr += log(Arrhenius + SMALL_S) -
E_Activation/(UNIVERSAL_GAS_CONSTANT*T);
/* 1.e-40 < rr < 1.e40 */
*rr = MAX(*rr,-40);
*rr = MIN(*rr,40);
*rr = exp(*rr);

[AlexanderZ]

how do you know there is a problem?
put here all logs you have

best regards