[dmoroian]

I wrote a udf which reads the temperature field using macro C_T(c,t). However, if the energy equation is not enabled the udf generates a SEGMENTATION FAULT.
Is it possible to check from within the udf if there is a temperature field?

Dragos

[gearboy]

Originally Posted by dmoroian

I wrote a udf which reads the temperature field using macro C_T(c,t). However, if the energy equation is not enabled the udf generates a SEGMENTATION FAULT.
Is it possible to check from within the udf if there is a temperature field?

Dragos

#include "udf.h"
DEFINE_ON_DEMAND(test_temp_field)
{
Domain *domain;
Thread *c_thread;
domain=Get_Domain(1);
thread_loop_c(c_thread, domain) /*loops over all cell threads in domain*/
{

if(NNULLP(THREAD_STORAGE(c_thread,SV_T)))
Message0("Temp. field exists\n");
else
Message0("Temp. field doesn't exist\n");

}

[dmoroian]

Thanks Ji,
Do you have a similar check for species?

Dragos

[gearboy]

Originally Posted by dmoroian

Thanks Ji,
Do you have a similar check for species?

Dragos

Use the following for species transport model.
if(NNULLP(THREAD_STORAGE(c_thread,SV_Y_I(i))))

i is the number of your specified species.

[dmoroian]

SV_Y_I is not working (sintax error), instead I think SV_Y_0 should be used and then check for particular SV_Y_I(i)

[dmoroian]

Originally Posted by dmoroian

SV_Y_I is not working (sintax error), instead I think SV_Y_0 should be used and then check for particular SV_Y_I(i)

The above test gives false positives, the following check seems to work well (12.0.16):

if(!rp_spe)
   Error("Species transport is disabled!\n");
if(2 > n_spe)
   Error("Not enough species, minim 3, current %d.\n",n_spe);