Out of bound in Simulation::removeAccretedGas

In Simulation::removeAccretedGas

case 3:
      h_min = 0;
      h_max = ni;
      while((phi(h_min) < vAngle - 2*RRoche/planetDist) && (h_min < ni)) h_min++;
      while((phi(h_max) > vAngle + 2*RRoche/planetDist) && (h_max > 0)) h_max--;
      break;
.....
Will fail potentially (and in some of our cases) there:
              real vpcell = flat ? real(0) : 0.5*(vphi(i,h,j)+vphi(i,h+1,j));