#160;          vv(1:3,j) = rp%cx(i)%p(1:3,j)/rp%cx(i)%mass(j)</div>
<div class="line"><a name="l00265"></a><span class="lineno">  265</span>&#160;<span class="keywordflow">        enddo</span></div>
<div class="line"><a name="l00266"></a><span class="lineno">  266</span>&#160;        v(i,:) = reshape(vv,(/rp%na*3/))</div>
<div class="line"><a name="l00267"></a><span class="lineno">  267</span>&#160;        <span class="comment">!print*, &#39;INITIAL VALUES = &#39;, norm2(v(i,:)), norm2(reshape(rp%cx(i)%forc