Hi.
I need to find the derivative of a variale (reactant.molconc) at particular location(x=0.002) which i m able to do it and then specify two sources, NEGATIVE SOURCE at left of the location and POSITIVE SOURCE at the right of the location, where source is some simple function of the gradient of the reactant.molconc .
The geometry is simple reactangle with a boundary at x=0.002
I thougt of a way but don't know(90% sure that it won't work) if it will work
the part of the code i written is :
DO ILOC = 1, NLOC
IF (X(ILOC).GE.0.002 .AND. X(ILOC).LE.0.002 .AND.
& Y(ILOC).GE.0.00 .AND. Y(ILOC).LE.0.001) THEN
SOURCE(ILOC-1) = -(1e-9)*GRAD_PHI(1,ILOC)
SOURCE(ILOC+1) = (1e-9)*GRAD_PHI(1,ILOC)
END IF
END DO
but i m not sure if ILOC-1 will be just left of the location and ILOC+1 at the right of the boundary. The problem is how the vertices are assigned to NLOC..
Plz help on this...Suggest some other and better ways to do the same.

Hi.
but i m not sure if ILOC-1 will be just left of the location and ILOC+1 at the right of the boundary. The problem is how the vertices are assigned to NLOC..

No, it's not going to work. ILOC has nothing to do with the physical location of a vertex. More importantly, you cannot define a point or line source in the finite volume method. You should add your source like the "User Defined Momentum Source" example in the documentation.