Boundary condition in `PLUTO`

init.c

// Following must be put in `definitions.h` for this to work
/* #define  INTERNAL_BOUNDARY              YES */
if (side == 0) { // -- select active cells --
  RBox dom_box;
  int convert_to_cons = 0;
  DOM_LOOP(k,j,i) {
    int num_cells = 1;
    // IDIR --> x-dirction
    if ( (grid->xgc[IDIR][i] >= (grid->xl_glob[IDIR][grid->gend[IDIR]-num_cells+1])) 
      && (grid->xgc[IDIR][i] <= (grid->xr_glob[IDIR][grid->gend[IDIR]])) ){
      /* Put your conditions in here */
      // d-Vc[RHO][k][j][i] = 
      convert_to_cons = 1;
    }
    if (convert_to_cons) {
      RBoxDefine (i, i, j, j, k, k, CENTER, &dom_box);
      PrimToCons3D(d->Vc, d->Uc, &dom_box);
    }
  }
}