gs_b function

Default boundary generation