Searching for Domain in Define VR

Arthur Valais

Has anyone used DEFINE_VR_RATE?

1. The following trivial example results in a non-zero reaction rate. What am I missing here?

2. Actually my problem is more complex. I want to set the reaction rate to be zero in all zones except ignitionZone (id: 9)

If I try...

Thread *ignitionThread = Lookup_Thread(domain,ignitionZone);

it should work, however I have no way of finding the value *domain since it didn't come along with the function call.

Any work-arounds?

Cheers,

Arthur Valais

Trivial function attached...

DEFINE_VR_RATE(zero_rate,cell,thread,r,mole_weight ,species_mf,rate,rr_t) {

cell_t c;

begin_c_loop_all (c,thread)

{

*rate = 0.;

*rr_t = 0.;

}

end_c_loop_all (c,thread) }