do i = 2, i2Tmp
do j = 2, j2Tmp


    ! ************************************************************************************************************
    ! Step 1a): Locate point that is farthest away from wall (at far field (FF)) and save its flow quantities
    ! ************************************************************************************************************


    ! loop over all points of the wall normal "pillar", locate local maximum of distx within the far field range
    maxDistxLocalFF = 0.
    idxMaxDistxLocalFF = 2

    do k = 2, k2Tmp

        if (distxTmp(i,j,k) > maxDistxLocalFF .and. distxTmp(i,j,k) <= maxDistxAllowed) then

            maxDistxLocalFF = distxTmp(i,j,k)
            idxMaxDistxLocalFF = j

        endif
    enddo

    ! depending on the orientation of the block, the index must be inversed,
    ! and if the far field range is reached, those cells are not modified
    if (wallNormDirIJK >0) then
        idxB = 2
        idxE = idxMaxDistxLocalFF
    else
        idxB = idxMaxDistxLocalFF
        idxE = k2Tmp
    endif

    ! maximum global wall distance of far field cell
    ! -> also consider dummy layers that contain the information of neighboring blocks
    idxMaxDistxGlobalFF = maxloc(vortShieldStep1Tmp(i,j,1:k2Tmp+1,1), 1)
    maxDistxGlobalFF = vortShieldStep1Tmp(i,j,idxMaxDistxGlobalFF,1)