ij * SMAX(0.0, face_wall_external_acceleration);</div><div class="line"><a name="l00170"></a><span class="lineno">  170</span>&#160;                        Real rho_in_wall = fluid_k-&gt;DensityFromPressure(p_in_wall);</div><div class="line"><a name="l00171"></a><span class="lineno">  171</span>&#160;                        Vecd vel_in_wall = 2.0 * vel_ave_i - vel_n_k[index_j];</div><div class="line"><a name="l00172"></a><span class="lineno">  172</span>&#160;</div><div class="line"><a name="l00173"></a><sp