[joshmccraney]

Hi All!

I am running OpenFOAM 6. When executing funkySetFields -time 0 I get the error:

/*---------------------------------------------------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  6
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
Build  : 6-fa1285188035
Exec   : funkySetFields -time 0
Date   : Sep 27 2020
Time   : 23:57:05
Host   : "josh-Super-Server"
PID    : 833
I/O    : uncollated
Case   : /home/josh/OpenFOAM/josh-6/run/cases/ICF4/tank_2/case
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
swakVersion: 0.4.3 (Release date: Next release)
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0


Time = 0
 Using funkySetFieldsDict 



Part: initFieldAlpha
 Modifying field alpha.water of type volScalarField

 Putting "0" into field alpha.water at t = "0" if condition "true" is true
 Keeping patches unaltered

 Setting 104904 of 104904 cells
 Writing to "alpha.water"


Part: initFieldU
 Modifying field U of type volVectorField

 Putting "vector(0.0,0.0,0.0)" into field U at t = "0" if condition "true" is true
 Keeping patches unaltered

 Setting 104904 of 104904 cells
 Writing to "U"


Part: pressureAir
 Modifying field p_rgh of type volScalarField

 Putting "0" into field p_rgh at t = "0" if condition "true" is true
 Keeping patches unaltered

 Setting 104904 of 104904 cells
 Writing to "p_rgh"


Part: alphaIC
 Modifying field alpha.water of type volScalarField

 Putting "1" into field alpha.water at t = "0" if condition "pos().y<0.0006+0.00468023*pow(459.772*pos().z-19944.1*pow(pos().z,8/3)+177465*pow(pos().z,13/3),1/3)" is true
 Keeping patches unaltered

#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3  ? in "/lib/x86_64-linux-gnu/libm.so.6"
#4  pow in "/lib/x86_64-linux-gnu/libm.so.6"
#5  Foam::pow(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
#6  Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::pow<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
#7  parserField::FieldValueExpressionParser::parse() at ??:?
#8  Foam::FieldValueExpressionDriver::parseInternal(int) at ??:?
#9  Foam::CommonValueExpressionDriver::parse(Foam::exprString const&, Foam::word const&) at ??:?
#10  ? in "/home/josh/OpenFOAM/josh-6/platforms/linux64GccDPInt32Opt/bin/funkySetFields"
#11  ? in "/home/josh/OpenFOAM/josh-6/platforms/linux64GccDPInt32Opt/bin/funkySetFields"
#12  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#13  ? in "/home/josh/OpenFOAM/josh-6/platforms/linux64GccDPInt32Opt/bin/funkySetFields"
Floating point exception (core dumped)

Do you have any idea what the issue is? Tough for me to see. The odd thing is, funkySetFields works on my computer for other cases, so I copied a working case into the new directory and ran it: it runs fine. Then, I copy the geometry files and blockMesh etc. from the case I want to run, and I get the above errors. It's odd because when I change to setFields (rather than funkySetFields) the case runs. I've tried self-diagnosing this issue, but I have no clue what this error implies. Help is greatly appreciated.

EDIT: FIXED IT. I THINK the error was a domain issue with the "pos().y<0.0006+0.00468023*pow(459.772*pos().z-19944.1*pow(pos().z,8/3)+177465*pow(pos().z,13/3),1/3)" condition (fractional exponents were ill-defined outside a domain).

Thanks though!