[CIG-MC] Reflective boundary condition in CitcomCU

[ting yang]

Hi All,

In my understanding, reflective boundary condition means that Vn=0 at the
boundary layer.
Here Vn means velocity normal to the boundary.

I noticed that in subroutine velocity_refl_vert_bc in Boundary_condition.c
in CitcomCU, the reflective boundary condition is set by making Vx=0 and
Szx=0 (taking x direction as an example):

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;
                }

I wonder why it sets the reflective boundary condition like this.

How about commenting the lines about E->VB[3] or aading other lines like
setting Syx=0, just like below?

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

OR

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;
                }

Thank you very much!
-- 
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/20130122/989f2e5c/attachment.htm