[CIG-MC] Reflective boundary condition in CitcomCU

[ting yang]

Hi Zhong,

Thank you for all your  assistance.

What you said just confused me.
How could these three different kinds of settings all work well as the
reflective boundary condition?

Do you mean that the difference in settings below has little effect on the
convection field?
If this is what you mean, I will do some experiments to test this idea.

Reflective boundary condition for the first nodes in x direction:
(1) seted in CitcomCU

if(E->parallel.me_loc[1] == 0)
                {
                    E->VB[1][node1] = 0.0;
                    if((ii != 1) && (ii != E->mesh.noz))
                        E->VB[3][node1] = 0.0;
                }

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

OR (3)

if(E->parallel.me_loc[1] == 0)
                {
                    E->VB[1][node1] = 0.0;
                    if((ii != 1) && (ii != E->mesh.noz))
                        E->VB[3][node1] = 0.0;
                    if((XXXXXXXX) && (XXXXXXXX))
                        E->VB[2][node1] = 0.0;
                }


Thanks,

Best regards!

-- 
wish you happiness !
yours sincerely,
Yang Ting
School of Earth and Space Science
University of Science and Technology of China
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://geodynamics.org/pipermail/cig-mc/attachments/20130123/f7f2ca47/attachment.htm