[alanlove]

Hi,everyone
I am writing a marco of DEFINE_SOURCE as followings:
DEFINE_SOURCE(mass_source,c, t, dS, eqn)
{
real source_mass;
if (C_T(c,t)>800.0)
{
source_mass=Density*PI*area*root/C_VOLUME(c,t);
dS[eqn]=0.0;
}
else
dS[eqn]=source_mass=0.0;
return source_mass;
}
My problem is: In order to validate my codes, How can I get the value of source_mass?Any ideas,thank u very much.

[nonyear]

you can do that
(1) Use Messsage function (if the mesh file is small)
Add
Message(" source_mass = %g\n", source_mass);
before

return source_mass;

then you can get the message

(2) Through a UDM(if the mesh file is big)
define a UDM
then add
C_UDMI(c,t, 0)=source_mass;
before
return source_mass;

then you can get the message through display the coutour

[alanlove]

Big Thanks

[nonyear]

you are welcome.