TABLE OF CONTENTS

1. setups/twogamma/user_bnd

setups/twogamma/user_bnd

    tot_bnd

         call user_bnd(solnData,      block_no,neigh  , bndcnd)
         call user_bnd(real(:,:,:,:), integer, integer, integer)
 

              Called for a block if the boundary condition is bnd_user.
              The routine must fill the appropriate guardcells (as determined
              by `neigh', which is 1 for -X, 2 for +X, 3 for -Y, 4 for +Y,
              5 for -Z, 6 for +Z) for the given block number (block_no).
              For completeness, the boundary condition is passed as well;
              at the moment, we expect this to be bnd_user.

              This is a version with inflow boundaries on the -x side of the 
              domain borrowed from the wind tunnel.