Two things do jump out as slightly odd in this script. Not that they will make any difference to the result.
On Jan 22, 2008 4:43 PM, asreeve <[EMAIL PROTECTED]> wrote: > discharge=dot(head.getFaceGrad(),kond.getHarmonicFaceValue()) kond.getHarmonicFaceValue() is a scalar and head.getFaceGrad() is a vector. Does it make any sense to dot them? Do you not just mean multiply? > transportEq2=(TransientTerm()== > -PowerLawConvectionTerm(coeff=vel()) > -DiffusionTerm(coeff=Disp)) I think that vel() evaluates the variable and returns a numpy array. This means that the coefficient will not be updated if kond or head change. I realize that vel never changes in this calculation. Cheers -- Daniel Wheeler
