[CIG-MC] Replicating in/out flow boundary conditions from CitcomCU in CitcomS

[scott tarlow]

Hello all,

I have been trying to find a way to to open up the the two xz planes in my regional mesh, having a flow in/out boundary condition. This seems to be done in CitcomCU, but when looking at the code in Boundary_conditions.c from the CitcomCU-1.0.2-br-inflow package, the changes from the reflective boundary condition did not seem obvious, though my programming experience is limited. Is there a way to implement this in CitcomS ? In CitcomS's Regional_boundary_conditions.c, I opened up the bottom of the mesh by adding this code:

 else if(E->mesh.botvbc == 2) {
    /* This extra BC is for a open bottom */
        horizontal_bc(E,E->sphere.cap[j].VB,1,1,0.0,VBX,0,lv,j);
        horizontal_bc(E,E->sphere.cap[j].VB,1,3,0.0,VBZ,0,lv,j);
        horizontal_bc(E,E->sphere.cap[j].VB,1,2,0.0,VBY,0,lv,j);
        horizontal_bc(E,E->sphere.cap[j].VB,1,1,E->control.VBXbotval,SBX,1,lv,j);
        horizontal_bc(E,E->sphere.cap[j].VB,1,3,0.0,SBZ,1,lv,j);
        horizontal_bc(E,E->sphere.cap[j].VB,1,2,E->control.VBYbotval,SBY,1,lv,j);
        }


and am wondering if there is a way to modify the existing reflective boundary condition code to open up the sides?

Thank you for your time,

Scott Tarlow
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://geodynamics.org/pipermail/cig-mc/attachments/20130520/3c6d7878/attachment.html>