UPPER branch at ReP=1517.073 (Pressure-gradient based Reynolds number: ReP= |dP/dx| h3/(2 nu2), where h is half-channel height).
LOWER branch at ReP=1513.24
This is only half of the data for -1< y <0. The other half can be rebuilt by symmetry i.e.
U(y)=U(-y), Vy(x, z, y)= -Vy(x, z, -y) and Oy(x, z, y)=Oy(x, z, -y).
The data has been reconstructed at the half
Gauss-Lobatto points for resolution nx=32, nz=32, ny=33
(33 for full channel, 17 to 33 for bottom half channel).
To read the full field in FORTRAN save the above data to "fort.9",
"fort.10" and "fort.11" and read with
do jy=17,33 read(9,110) y(jy), U(jy) do jz=1,32 do jx=1,32 write(10,110) Vy(jx,jz,jy) write(11,110) Oy(jx,jz,jy) enddo enddo enddo 110 format((G23.16))