[seasoul]

I am trying to define a udf for surface injection and I think it would be using

DEFINE_DPM_INJECTION_INIT

, and a minimal test as below works, by simply defining the particle diameter and flow rate.

DEFINE_DPM_INJECTION_INIT(SURFINJ, I)
{
Particle *p;
face_t f;
Thread *cthread;

loop(p,I->p)  
    {
      P_DIAM(p) = 1e-5;
      P_FLOW_RATE(p) = 1.5e-3;
    }

}

I have found that the number of streams is the number of faces of the boundary surface, exactly the same as when we define the injection by surface without udf. But if I want to define the number of the streams, for example, twice of the number of the faces and then define the location for each stream, which function should I call?

Thank you in advance.