="n">Bmag</span><span class="p">(</span><span class="n">cc</span><span class="p">)</span><span class="o">*</span><span class="n">Bmag</span><span class="p">(</span><span class="n">cc</span><span class="p">))</span>
     
  <span class="k">end do</span>
  <span class="c">!$OMP END SIMD</span>
  
   
<span class="c">!  write(output_unit_write,&#39;(&quot;PSIp: &quot;,E17.10)&#39;) PSIp</span>
  
<span class="c">!  write(output_unit_write,&#39;(&quot;Y_R: &quot;,E17.10)&#39;) Y_R</span>
<span class="c">!  writ