Searching for Domain in Define VR
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) } |
All times are GMT -4. The time now is 07:07. |