359   PetscInt    ilim = grid->ni - 1, jlim = grid->nj - 1;  in CalcSecInv()  local
364   if (i == ilim) i--;  in CalcSecInv()
521   PetscInt    ilim = grid->ni - 1;  in ZMomentumResidual()  local
531     if (i == ilim) TE = TW;  in ZMomentumResidual()
548   if (i == ilim) dwdxE = etaE * (x[j][i].u - x[j + 1][i].u) / dz;  in ZMomentumResidual()
588   PetscInt    ilim = grid->ni - 1, jlim = grid->nj - 1, jlid = grid->jlid;  in EnergyResidual()  local
617   if (param->adv_scheme == ADVECT_FV || i == ilim - 1 || j == jlim - 1 || i == 1 || j == 1) {  in EnergyResidual()
646   PetscInt    ilim = grid->ni - 1, jlim = grid->nj - 1;  in ShearStress()  local
649   if (j <= grid->jlid || i < j || i == ilim || j == jlim) return EPS_ZERO;  in ShearStress()
681   PetscInt    ilim = grid->ni - 1, jlim = grid->nj - 1, ivisc;  in XNormalStress()  local
699     if (i == ilim || j == jlim) return EPS_ZERO;  in XNormalStress()
720   PetscInt    ilim = grid->ni - 1, jlim = grid->nj - 1, ivisc;  in ZNormalStress()  local
738     if ((i == ilim) || (j == jlim)) return EPS_ZERO;  in ZNormalStress()
1100   PetscInt   i, j, is, js, im, jm, ilim, jlim, ivt;  in ViscosityField()  local
1118   ilim = grid->ni - 1;  in ViscosityField()
1126       if (i < ilim && j < jlim) {  in ViscosityField()
1294   PetscInt    i, j, mx, mz, ilim, jlim;  in FormFunctionLocal()  local
1301   ilim = mx - 1;  in FormFunctionLocal()
1318       } else if (i == ilim) {  in FormFunctionLocal()
1357       } else if (i == ilim) {  in FormFunctionLocal()
1377       } else if ((i == ilim || j == jlim) && ibound == BC_ANALYTIC) {  in FormFunctionLocal()
1395       } else if (i == ilim) {  in FormFunctionLocal()