<div dir="ltr"><div><div><div><div><div><div>Hi Zhong,<br><br></div>Thank you for all your 
assistance.<br></div><br>What you said just confused me.<br></div>How could these three different kinds of settings all work well as the reflective boundary condition?<br><br></div>Do you mean that the difference in settings below has little effect on the convection field?<br>
</div><div>If this is what you mean, I will do some experiments to test this idea.<br></div><div><br></div>Reflective boundary condition for the first nodes in x direction:<br></div>(1) seted in CitcomCU<br><div><div><pre>
if(E-&gt;parallel.me_loc[1] == 0)
                {
                    E-&gt;VB[1][node1] = 0.0;
                    if((ii != 1) &amp;&amp; (ii != E-&gt;mesh.noz))
                        E-&gt;VB[3][node1] = 0.0;
                }

(2) 
if(E-&gt;parallel.me_loc[1] == 0)
                {
                    E-&gt;VB[1][node1] = 0.0;
                }

OR (3)

if(E-&gt;parallel.me_loc[1] == 0)
                {
                    E-&gt;VB[1][node1] = 0.0;
                    if((ii != 1) &amp;&amp; (ii != E-&gt;mesh.noz))
                        E-&gt;VB[3][node1] = 0.0;
                    if((XXXXXXXX) &amp;&amp; (XXXXXXXX))
                        E-&gt;VB[2][node1] = 0.0;
                }</pre><br></div><div>Thanks,<br><br></div><div>Best regards!<br></div><div><div><div><div><div><div><br>-- <br>wish  you  happiness !<br>yours  sincerely,<br>Yang Ting<br>School of Earth and Space Science <br>
University of Science and Technology of China<br>
</div></div></div></div></div></div></div></div>