medium.alpha_power;

<span class="comment">% extract the pressure amplitude at each position</span>
p = extractAmpPhase(sensor_data.p, 1/kgrid.dt, freq);

<span class="comment">% reshape the data, and calculate the volume rate of heat deposition</span>
p = reshape(p, Nx, Ny);
Q = alpha_np .* p.^2 ./ (medium.density .* medium.sound_speed);
</pre>

<p>The volume rate of heat deposition is assigned to <code>source.Q</code>, and <code><a href="kWaveDiffusion.html">kWaveDiffusion</a></code> is used to simulate t