I'm trying to set wall BC as rough solid face in uscpl file (coal combustion) with roughness height 1cm:
Code: Select all
CALL GETFBR('wall',NLELT,LSTELT)
!==========
do ilelt = 1, nlelt
ifac = lstelt(ilelt)
! Kind of boundary conditions for standard variables
itypfb(ifac,iphas) = iparug
rcodcl(ifac, iu(iphas), 3) = 0.01d0
! zone's number (from 1 to n) izone = 5
! - Allocation of the zone's number to the face
izfppp(ifac) = izone
enddoUnfortunately I receive the following error:
@
@@ WARNING: ABORT DURING THE BOUNDARY CONDITIONS VERIF.
@ ========
@
@ Uninitialized boundary conditions : 0
@ Unexpected boundary conditions:
@ on the scalars : 0
@ on the scalars representing
@ a variance : 0
@ Incoherencies:
@ between velocity and scalars : 240
@
@ The calculation will not be run.
@
@ Verify the parameters given via the interface or
@ usclim.
@
Call stack:
1: 0x7f8bb9d6a844 <vericl_+0x5d94> (libsaturne.so.0)
2: 0x7f8bb9b852a2 <condli_+0x15c2> (libsaturne.so.0)
3: 0x7f8bb9d3c821 <tridim_+0x6931> (libsaturne.so.0)
4: 0x7f8bb9b74855 <caltri_+0x5085> (libsaturne.so.0)
5: 0x7f8bb9b4fc63 <cs_run+0x8a3> (libsaturne.so.0)
6: 0x7f8bb9b4ff45 <main+0x1f5> (libsaturne.so.0)
7: 0x7f8bb6b80abd <__libc_start_main+0xfd> (libc.so.6)
8: 0x402c69 ? (?)
End of stack
As far as I'm concerned it arises from incoherency between velocity and some scalars. Should I add specific conditions for scalars when using iparug BC?
Thank you in advance
Norbert