[morourke]

I am trying to implement a boundary condition for a transported species named RP. As I understand it the boundary condition is a flux boundary condition j(x,t) and is simply a linear function of the species RP.

j(x,t) = K . RP(x,t)

where K is a constant.

Any advice on how I implement this?

Thanks in advance.