   int s2b = dom->Gfy._s2b;
 
   if((ti < dom->Gfy._inb) && (tj < dom->Gfy._jnb)) {
-    //v[ti + tj*s1b + dom->Gfy._ksb*s2b] = 2. * bc[ti + tj*dom->Gfy.inb]
-    //  - v[ti + tj*s1b + dom->Gfy._ks*s2b];
 	v[ti + tj*s1b + dom->Gfy._ksb*s2b] = bc_b[ti + tj*dom->Gfy.inb];
 	v[ti + tj*s1b + dom->Gfy._ks*s2b] = bc_t[ti + tj*dom->Gfy.inb];
  }

  int s2b = dom->Gfz._s2b;

  if((ti < dom->Gfz._inb) && (tk < dom->Gfz._knb)) {
    w[ti + (dom->Gfz._jeb-1)*s1b + tk*s2b] = bc_s[tk + ti*dom->Gfz.knb];
    w[ti + dom->Gfz._jeb*s1b + tk*s2b] = bc_n[tk + ti*dom->Gfz.knb];
  }
}


  int s2b = dom->Gfy._s2b;

  if((ti < dom->Gfy._inb) && (tj < dom->Gfy._jnb)) {
    v[ti + tj*s1b + (dom->Gfy._ke-1)*s2b] = bc_s[ti + tj*dom->Gfy.inb];
    // velocity on ghost cell
    v[ti + tj*s1b + (dom->Gfy._keb-1)*s2b] = bc_n[ti + tj*dom->Gfy.inb];
  }
}

