radB</span><span class="p">,</span><span class="n">vars</span><span class="p">%</span><span class="n">curlb</span><span class="p">,</span> <span class="n">vars</span><span class="p">%</span><span class="n">flagCon</span><span class="p">,</span><span class="n">vars</span><span class="p">%</span><span class="n">PSI_P</span><span class="p">)</span>

       <span class="k">end if</span>

<span class="k">    </span><span class="n">endif</span>

  <span class="k">end subroutine </span><span class="n">get_analytic