Source term in the mass balance equation

Mehdi BEN HAJ

Hi,

I want to introduce an internal source term (q) in the mass balance equation for generating a water wave.

d(rho)/dt + div(rho*V) = q (1)

I think that the Plant command is:

PLANTBEGIN

**Instruction for q calculation (not given)

PATCH (SP1,PHASEM,IX1,IX2,1,NY,1,NZ,1,LSTEP)

<SORC01> VAL=q

COVAL (SP1,P1 , FIXFLU , GRND )

PLANTEND

Could some one confirm me that this Plant command is right to introduce a source term in equation (1)?

Thanks.