span class="p">(</span><span class="o">-</span><span class="p">(</span><span class="n">rm</span><span class="o">/</span><span class="n">h_params</span><span class="p">%</span><span class="n">sigma_Z</span><span class="p">)</span><span class="o">**</span><span class="mi">2</span><span class="o">/</span><span class="mi">2</span><span class="p">)</span>

<span class="c">!    write(output_unit_write,&#39;(&quot;rm: &quot;,E17.10)&#39;) rm</span>
    
    <span class="n">A</span> <span class="o">=</span> <span c