span><span class="p">))</span>     <span class="c">! MRC</span>
    <span class="c">!	feta = 0.5_rp*A*EXP(A*COS(deg2rad(eta)))/SINH(A)                            ! MRC</span>

    <span class="n">fRE_H_3D</span> <span class="o">=</span> <span class="n">fRE_H_3D</span><span class="o">*</span><span class="n">feta</span>

    <span class="c">!if (fRE_H_3D.eq.0._rp) then</span>
    <span class="c">!   write(6,*) &#39;f_RE_H&#39;,f0 + m*(g - g0)</span>
    <span class="c">!   write(6,*) &#39;feta&#39;,feta</span