[T.D.]

hi foamers,

I have a volVectorField Wrel, and a volScalarField gammadot (of same dimensions)

I need to get the volScalarField rhok,
where rhok= (2*magnitude(Wrel) ) / (gammadot+magnitude(Wrel))
so i tried the following:

volScalarField rhok=2*mag(Wrel)/(gammadot+mag(Wrel));

but it didn't work !!!! it compiles well, but when i run it it gives this error:

.
.
.
nProcs : 1
SigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0

Reading transportProperties

Reading field p

Reading field U

Reading field c

Reading field gammadot

Reading/calculating face flux field phi

Reading field SigmaL

Reading field SigmaR

Reading field SigmaP

Reading field Sigma

Reading field SigmaF

Reading field J


Starting time loop

Time = 0.1

Courant Number mean: 0 max: 1.71165
#0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam171/lib/linuxGccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/opt/openfoam171/lib/linuxGccDPOpt/libOpenFOAM.so"
#2 Uninterpreted:
#3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/opt/openfoam171/lib/linuxGccDPOpt/libOpenFOAM.so"
#4 void Foam::divide<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/home/talib/OpenFOAM/talib-1.7.1/applications/bin/linuxGccDPOpt/TalibicoFoam2DframeInvariant"
#5 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam:perator/<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<doub le, Foam::fvPatchField, Foam::volMesh> > const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) in "/home/talib/OpenFOAM/talib-1.7.1/applications/bin/linuxGccDPOpt/TalibicoFoam2DframeInvariant"
#6
in "/home/talib/OpenFOAM/talib-1.7.1/applications/bin/linuxGccDPOpt/TalibicoFoam2DframeInvariant"
#7 __libc_start_main in "/lib/tls/i686/cmov/libc.so.6"
#8
in "/home/talib/OpenFOAM/talib-1.7.1/applications/bin/linuxGccDPOpt/TalibicoFoam2DframeInvariant"
Floating point exception
talib@jml-port:~/Desktop/talib/OpenFOAM/talib-1.7.1/run/tutorials/incompressible/icoFoam/ChannelFlow2Dfi30percent_constantInletVelocityTali bicoFoam2DframeInvariant$

any ideas ???

thanks

T.D.