[mps]

I have seen some people looking for dpm_vof as a macro in udf. I provided the following udf, writing dpm-vof of each cell in an udm. I hope everybody finds it useful

DEFINE_ADJUST(dpm_vof, domain)
{
real Rop;
cell_t c;
Thread *t;
Particle *p;
Injection *Ilist = Get_dpm_injections();
Injection *I;
loop(I, Ilist)
{
loop(p, I->p)
{
Rop = P_RHO(p);
}
}

thread_loop_c(t, domain)
{
begin_c_loop(c, t)
{
C_UDMI(c, t, 0) = C_DPMS_CONCENTRATION(c, t) / Rop; //dpm_vof saved in a udm
//this macro works as well: C_STORAGE_R(c, t, SV_DPMS_CONCENTRATION)
}
end_c_loop(c, t);
}
}

[mrishi]

Sorry to pull up an old post.



But I am facing problems in implementing a very similar piece of code.
I store the value of C_DPMS_CONCENTRATION(c,t) in C_UDSI(c,t,0)
but when I plot value of this scalar I obtain garbage values (negative values and values which are not same as the variable DPM Concentration available in contour list). i cannot understand the reason.

[souza.emer]

Originally Posted by mps

I have seen some people looking for dpm_vof as a macro in udf. I provided the following udf, writing dpm-vof of each cell in an udm. I hope everybody finds it useful

DEFINE_ADJUST(dpm_vof, domain)
{
real Rop;
cell_t c;
Thread *t;
Particle *p;
Injection *Ilist = Get_dpm_injections();
Injection *I;
loop(I, Ilist)
{
loop(p, I->p)
{
Rop = P_RHO(p);
}
}

thread_loop_c(t, domain)
{
begin_c_loop(c, t)
{
C_UDMI(c, t, 0) = C_DPMS_CONCENTRATION(c, t) / Rop; //dpm_vof saved in a udm
//this macro works as well: C_STORAGE_R(c, t, SV_DPMS_CONCENTRATION)
}
end_c_loop(c, t);
}
}

Thank you very very much sr!