e=itimes(i);

    <span class="comment">% SHOULDN'T HAVE TO CHANGE BELOW HERE</span>
    tic
    <span class="comment">%initialize dataset object</span>
    nc=<a href="../../njtbxhelp/njTools/njTBX-2.0/@mDataset/mDataset.html">mDataset</a>(uri);
    <span class="comment">%get variable object</span>
    nc_var = nc{var};
    <span class="comment">% read data at specified time step (all nodes)</span>
    zeta=nc_var(itime,:);

    <span class="comment">% Find the coordinate variables</span>
    <