function [U,V,W]=read_ml

% This script reads in data from homogeneous isotropic turbulence DNS
% in Matlab. 
%
% Running 
% [u,v,w] = read_box reads in the three components of velocity,
%
% Make sure you run the script from the same folder where you saved 
% the data file, or specify the correct path to the data file below
%

% filename = '/path/to/data/file/ML_field_sp.dat';
filename = 'ML_field_sp.dat';

fid = fopen(filename,'r');

dims = fread(fid,4,'integer*4','ieee-le');

nx = dims(1);
ny = dims(2);
nz = dims(3);
nvar = dims(4);

nsize = nx*ny*nz;

dummy = fread(fid,nsize,'real*4','ieee-le');
U = reshape(dummy,nx,ny,nz);
dummy = fread(fid,nsize,'real*4','ieee-le');
V = reshape(dummy,nx,ny,nz);
dummy = fread(fid,nsize,'real*4','ieee-le');
W = reshape(dummy,nx,ny,nz);

fclose(fid);