 class="n">A</span> <span class="o">=</span> <span class="n">A</span><span class="o">*</span><span class="n">h_params</span><span class="p">%</span><span class="n">A_fact</span>

    <span class="c">!write(output_unit_write,*) &#39;A&#39;,A</span>
    <span class="c">!flush(output_unit_write)    </span>

    <span class="n">feta</span> <span class="o">=</span> <span class="n">A</span><span class="o">*</span><span class="nb">EXP</span><span class="p">(</span><span class="o">-</span><span class="n">A</span><s