I would like to make a species source at a Wall boundary in COMET. I am trying to use spesou.f, but getting nowhere. Can someone write a brief example where the faces of region 1 generate 0.001 kg/sec of species 1 please? Thanks in advance.

EITHER you use spesou.f and in that case you have to set the source term in the cell next to the wall, since spesou.f is called for the fluid domain (in which case you need to switch on the flag -- see the SPECIE command -- to be able to use this subroutine)

OR you use the subroutine userbc.f (in which case you need to set the user subroutine flag on -- see the BDEF command in Cometpp) to set boundary values

Thank you for your help. I figureed out last night that spesou.f would only specify for cells. Do I use the LK(K) variable to figure out the cell indices?

I would prefer using the wall boundaries themselves as the source, how do I do this using userbc.f? (Which variable?)

Thanks again.

I think you meant LB(K) for getting the cell next to the wall, where K is the boundary index. For imposing the value at the wall within userbc.f, give the value into cb(ibofs(is)+K) where is is the species index (for the rest, look at the file ListOfVariables in the userprog directory).