[Nishant_Sharma]

Hi all,

I have written a udf to calculate the current in a MFC but I am getting some error:
Error: / (devide) invalid argument (wrong type): not a number.
Error Object: nan

I think the error is in udf (attached with this thread). Can anyone please check report error if any.

Thank you

[zobekenobe]

Error 1:

begine_thread_loop(t,d)

should read

begin_thread_loop(t,d)


Error 2:

begine_c_loop(c,t)

should read

begin_c_loop(c,t)



The 'e' after the begin

[Nishant_Sharma]

Originally Posted by zobekenobe

Error 1:

begine_thread_loop(t,d)

should read

begin_thread_loop(t,d)


Error 2:

begine_c_loop(c,t)

should read

begin_c_loop(c,t)



The 'e' after the begin

Yaah, that I figured out after, and even after correcting it I am getting an error / (divide): not an argument

[AlexanderZ]

Try this code

#include "udf.h"
real Voc = 0.92;
real Rext = 100.0;
real k2 = 9.13e-9;
real k1 = 5.57e-5;
real r1 = 1.62e-6;
real r2 = 4.32e-6;
real F = 96485;
real A = 0.6;
real B = 0.4;
real Kac = 0.592;
real Kox = 0.004;
real R = 8.314;
real T = 300.0;
real den = 998.0;

DEFINE_ADJUST(current, d)
{
    Thread *t;
    cell_t c;
    face_t f;
    real Vac, Vox, I, Cac, Cox, YI_AC, YI_OX, mw[6], V;
    thread_loop_c (t,d)
    {
        begin_c_loop (c,t)
        {
            YI_AC = C_YI(c,t,0);
            YI_OX = C_YI(c,t,1);
            Cac = YI_AC*den/60.0;
            Cox = YI_OX*den/32.0;
            Vac = R*T*log(r1*(Kac+Cac)/(Cac*k1))/(A*F);
            Vox = R*T*log(r2*(Kox+Cox)/(Cox*k2))/((B-1)*F);
            V = Voc - Vac - Vox;
            I = V/Rext;
            C_UDMI(c,t,0) = Cac;
            C_UDMI(c,t,1) = Cox;
            C_UDMI(c,t,2) = Vac;
            C_UDMI(c,t,3) = Vox;
            C_UDMI(c,t,4) = V;
            C_UDMI(c,t,5) = I;
        }
        end_c_loop(c,t);
    }
}

best regards