_RL iflx(+1:+2)
if ( QUADROOT(aval,bval,cval,iflx) ) then
if ( ( iflx(1) .gt. -1. _d 0 )
& .and. ( iflx(1) .lt. +1. _d 0 ) ) then
& + iflx(1) * fhat(3) * 2. _d 0
& +(iflx(1) ** 2) * fhat(4) * 3. _d 0
& +(iflx(1) ** 3) * fhat(5) * 4. _d 0
if ( ( iflx(2) .gt. -1. _d 0 )
& .and. ( iflx(2) .lt. +1. _d 0 ) ) then
& + iflx(2) * fhat(3) * 2. _d 0
& +(iflx(2) ** 2) * fhat(4) * 3. _d 0
& +(iflx(2) ** 3) * fhat(5) * 4. _d 0
if (iflx.eq.1) then
if (iflx.eq.1.and.omegaC(i,j,k,bi,bj) .gt. 1.d0) then
INTEGER iflx
iflx=1