************************************************************************
************************************************************************
            DATA FROM DNS OF IMPINGING SHOCK WAVE BOUNDARY
                      LAYER INTERACTION (M = 2.28)
************************************************************************
************************************************************************

This file contains some information to read and use some data from the DNS
database reported in Pirozzoli and Bernardini, AIAA J. 2010 (in press)

************************************************************************
NOTE THAT ONLY ONE FLOW SAMPLE IS UPLOADED FOR NOW, FOR ILLUSTRATIVE PURPOSES. 
IF YOU ARE INTERESTED THE FULL DATASET, PLEASE CONTACT THE AUTHORS 
sergio.pirozzoli@uniroma1.it
matteo.bernardini@uniroma1.it
************************************************************************

There are two-subdirectories

"grid"         ====> contains grid files (ascii format)
"flow_samples" ====> contains a 3D flow sample (binary format, little-endian)

The parameters of the simulation are

Free-stream Mach number ====> M = 2.28
Reynolds number based on the inlet BL thickness ====> Re = 16750.
Incidence angle of the shock generator ====> phi = 8 degree
gamma = 1.4 

The grid is stored in formatted files named "grid_ib.dat".
The fields are stored in unformatted files named "field_ib.bin".
The index "ib" refers to a specific block of the grid.
The computational domain is subdivided into 64 blocks along the
stream-wise direction (x coordinate). The wall-normal direction
is indicated with y and the spanwise with z.
The overall size of the domain is 80.58 in the stream-wise direction (x),
12.89 in the wall-normal direction (y), and 6.47 in the span-wise
direction (z). The reference length is the boundary-layer thickness
at the inlet section (x=0).
Each block consists of nx = 61 X ny = 344 X nz = 261 points uniformly
distributed in the x and z direction and clustered towards the wall.
The last point of each block coincides with the first one of the following
block.

In the "flow_samples" directory you will find 64 files of a single solution
whereas the corresponding grid files are in the directory "grid".

The grid file corresponding to block ib is written in "grid_ib.dat".
You can read it as

**************************************************************
                        FORTRAN CODE
**************************************************************
       open(unit=11,file='grid_'//ib//'.dat',form='formatted')
       read(11,*) nx, ny, nz
       do k=1,nz
        read(11,*) z(k)
       enddo
       do j=1,ny(ib)
        read(11,*) y(j)
       enddo
       do i=1,nx(ib)
        read(11,*) x(i)
       enddo
       close(11)
**************************************************************

The arrays x,y,z contain the grid coordinates.
The files containing the conservative variables can be read as

**************************************************************
                        FORTRAN CODE
**************************************************************
       open(11,FILE='field_'//ib//'.bin',form='unformatted')
       rewind (11)
       read  (11) telaps
       read  (11) ((((w(i,j,k,m),i=1,nx),
      .                          j=1,ny),
      .                          k=1,nz),
      .                          m=1,5)
       close  (11)
**************************************************************

where telaps is the elapsed time.
The order of variables in the array w is the following

 m  variable
------------
 1  rho
 2  rho*u
 3  rho*v
 4  rho*w
 5  p

where rho is the density, u, v, w are the velocity components and
p is the thermodynamic pressure.

All variables are nondimensional. The reference density, pressure,
temperature and velocity are taken to be the free-stream values.