[cico0815]

Hello everyone,

yesterday I downloaded the updated version of SOWFA from github. The solver
suite can be found here:

https://github.com/NREL/SOWFA

After several small tweaks everything compiles fine on my XUbuntu 14.04. I tried
the precursorABL/neutral tutorial case and after some updates (fixedFluxPressure
instead of buoyantPressure) I started the case. I used blockMesh, setFieldsABL
(which already failed because there was the following error)

Build  : 2.3.0-f5222ca19ce6
Exec   : setFieldsABL
Date   : Dec 05 2014
Time   : 14:08:06
Host   : "cfd"
PID    : 2650
Case   : /home/cfd/OpenFOAM/cfd-2.3.0/run/SOWFA/tutorials/precursorABL/neutral
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

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

Create mesh for time = 0

Reading field U
Reading field T
Reading field p_rgh
Creating/Calculating face flux field, phi...
Reading gravitational acceleration...
Selecting incompressible transport model Newtonian
Creating the kinematic density field, rhok...

u* = 0.35863751503 m/s

<U_1>  = (6.92820323026 3.99999999999 0)   <U_s> = (6.92820323026  3.99999999999 0)   <dU/dn> = (-1.16193474595e-18  -3.25610680393e-18 0)


--> FOAM FATAL ERROR: 
updateCoeffs(const scalarField& snGradp) MUST be called before updateCoeffs() or evaluate() to set the boundary gradient.

    From function fixedFluxPressureFvPatchScalarField::updateCoeffs()
    in file fields/fvPatchFields/derived/fixedFluxPressure/fixedFluxPressureFvPatchScalarField.C at line 151.

FOAM exiting

I anyways tried to let the ABLSolver run but an error alike the above occured. I
then commented out the line 325 in setFieldsABL.C (p_rgh.correctBoundaryConditions() and the setFieldABL utility ran afterwards.
I also commented the line 92 in ABLSolver.C out (see before) and the solver ran
up to this point

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.3.0                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 2.3.0-f5222ca19ce6
Exec   : ABLSolver
Date   : Dec 05 2014
Time   : 12:28:26
Host   : "cfd"
PID    : 10092
Case   : /home/cfd/OpenFOAM/cfd-2.3.0/run/SOWFA/tutorials/precursorABL/neutral
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

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

Create mesh for time = 0

Reading the gravitational acceleration, g...
Reading planetary rotation rate...
Reading the latitude...
Creating and reading potential temperature field, T...
Creating and reading modified pressure field, p_rgh...
Creating and reading velocity field, U...
Creating and calculating velocity flux field, phi...
Reading/calculating face flux field phi

Reading transport properties...
Selecting incompressible transport model Newtonian
Reading the atmospheric boundary layer properties...
     Specified wind at 80 m:
              Time 0,   wind from 240 degrees at 8 m/s
              Time 50000,   wind from 240 degrees at 8 m/s

     +x is east and +y is north
                   N
                   0
                   |

        W 270 --       --  90 E

                   |
                  180
                   S
Omega [0 0 -1 0 0 0 0] (0 5.40430167656e-05 4.86605508618e-05)
Creating turbulence model...
Selecting turbulence model type LESModel
Selecting LES turbulence model dynLagrangianCsBound
Selecting LES delta type smooth
Selecting LES delta type cubeRootVol
dynLagrangianCsBoundCoeffs
{
    filter          simple;
    ce              1.048;
    theta           1.5;
}

Creating the Coriolis force vector, fCoriolis...
Creating kinematic (Boussinesq) density field, rhok...
Creating the kinematic thermal conductivity field, kappat...
Reading and creating the wall shear stress field, Rwall...
Reading and creating the wall temperature flux field, qwall...
Creating and calculating the gravity potential field, gh and ghf...
Creating and calculating the static pressure field, p...
Setting up the pressure reference cell information...
Creating mean temperature field, Tmean...
Creating fluctuating temperature field, Tprime...
Creating mean velocity field, Umean...
Creating fluctuating velocity field, Uprime...
Creating mean SGS stress field, Rmean...
Creating mean SGS temperature flux field, qmean...
Creating mean SGS viscosity field, nuSGSmean...
Initializing with zero pressure gradient...

Courant Number mean: 0.373333333393 max: 0.373333333333

Total number of cell center height levels: 100
78400   90000000
...
78400   90000000
Local Flux Continuity Error:  Min 0   Max 3.39686165615e-14   Weighted Mean 8.87516162301e-17
Total Boundary Flux: 3.92910790963e-24

PIMPLE: Operating solver in PISO mode

Starting time loop

<U_1> = (8 0 0)   <U_s> = (8 0 0)   <dU/dn> = (0 0 0)
DILUPBiCG:  Solving for flm, Initial residual = 1, Final residual = 0.00180675762046, No Iterations 1
DILUPBiCG:  Solving for fmm, Initial residual = 1, Final residual = 0.00180675761966, No Iterations 1
TRef = 290
uStarMean = 0.817991099632   LMean = 9.99999999996e+29   phiMMean = 1
UParallelMeanMag = 8   UParallelPMeanMag = 8
RwMagMean  = 0.669109439077   RwMeanMag = 0.669109439077   sqrt(RwMagMean) =  0.817991099632   sqrt(RwMeanMag) = 0.817991099632   uStarMean =  0.817991099632
z1Mean = 4.99999999999
TRef = 290
TSurfaceMean = 300   TAdjacentMean = 300   deltaTMean = 0
#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2   in "/lib/x86_64-linux-gnu/libc.so.6"
#3   Foam::fixedHeatingRateFvPatchField::qwEvaluate(double&,  double&, double&, double&, double&, double, double,  double, double, double, double, double, double, double, double, double,  double, double, double, int) at ??:?
#4  Foam::fixedHeatingRateFvPatchField::evaluate(Foam::UPstream::commsTypes) at ??:?
#5   Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField,  Foam::volMesh>::GeometricBoundaryField::evaluate() at ??:?
#6  
 at ??:?
#7  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#8  
 at ??:?

Anyone out there who could help me resolve this issue? Is there a function that

is used to replace the p_rgh.correctBoundaryConditions()???

Any help is greatly appreciated!

Thanks a lot!

Best, Thomas