Call face value of a parameter
 

Hi FOAMERS
I am trying to solve a set of equations in the following form
--------
(Density)*fvm::ddt(U)
+ (Density/gL)*fvm::div(phi, U)
- (Visc)*(fvm::laplacian(U))
- GravityVector*g*Density*((BetaT*(T-TNot))+BetaC*(C-CNot))
+ fvm::Sp(gL*Visc/Perm,U)
--------
Perm is the permeability which is calculated by the following relation

Perm[celli] = (pow(SDAS.value(),2))*pow((gL[celli]),3)/(180*(pow((1-gL[celli]),2)));

Is there any way to make sure that the permeability at the cell face is calculated by mean harmonic interpolation?