ss="keyword">end</span>
            <span class="keyword">if</span> (~isempty(drv.particles))
                ravg = mean([drv.particles.radius]);
                rdev = std([drv.particles.radius]);
                rmin = min([drv.particles.radius]);
                rmax = max([drv.particles.radius]);
                <span class="keyword">if</span> (rdev/rmin &lt; 1e-8)
                    rdev = 0;
                <span class="keyword">end</span>
                fprintf(<span class="string">'Average radius