This readme file is for the experimental dataset associated with the publication "Pollutant dispersion around a single tall building"

File names with P1, P2, P3, describe the location of the field-of-view, with P1 being closest to the source. The source concentrations of all datasets have been normalised during PLIF calibration.
File names ending with a number then H, eg 2H, 4H, refer to stereo fields and the distance behind the tall building that each was taken. The X component of these fields denotes the distance from the dye source.
In all mean datasets the coordinate system is that x is downstream, y is cross-stream, and z is vertical. The dye release location is taken as the origin.
Instantaneous PLIF results are included in the instantaneous sub folder for the centreline data that was used to generate figure 8a.
Most data is stored as doubles, apart from the stereo fields between 2H and 12H. This is due to data size limitations.

All files apart from those in the instantaneous sub folder contain mean results. These slices contain results from both PIV and PLIF data and have the following fields:

x coordinate (streamwise): X
y coordinate (spanwise)  : Y
Z coordinate (wall-norm) : Z
Mean streamwise velocity : U
Mean span wise velocity  : V
Mean vertical velocity   : W
Mean concentration	 : C
Mean advective flux in x : CU
Mean advective flux in y : CV
Mean advective flux in z : CW
Mean turb flux in x      : cu
Mean turb flux in y    	 : cv
Mean turb flux in z    	 : cw
Mean u'u' (variance) 	 : uu
Mean u'v'		 : uv
Mean v'v'		 : vv
Mean u'w'		 : uw
Mean v'w'		 : vw
Mean w'w'		 : ww
Mean c'c' (C variance)	 : cc

All variables involving V are not present in P1, P2, and P3. This is due to these datasets being planar PIV results with no component in the y direction.

An example way to print concentration maps from the P1 mean centreline structure would be as follows:

imagesc((ACMP1.X(:,1)),fliplr(ACMP1.Z(1,:)),rot90(ACMP1.C));
set(gca,'YDir','normal')
axis equal

And an example for the 2H stereo structure:

imagesc(ACM2H.Y(:,1),fliplr(ACM2H.Z(1,:)),rot90(ACM2H.C));
set(gca,'YDir','normal')
axis equal

And an example for the P1 instantaneous slice:

imagesc(ACM1.AC_frame.X(1,:),(ACM1.AC_frame.Y(1,:)),(ACM1.AC_frame.C));
set(gca,'YDir','normal')
axis equal