359   PetscInt    ilim = grid->ni - 1, jlim = grid->nj - 1;  in CalcSecInv()  local
365   if (j == jlim) j--;  in CalcSecInv()
470   PetscInt    jlim = grid->nj - 1;  in XMomentumResidual()  local
476     if (j == jlim) TN = TS;  in XMomentumResidual()
493   if (j == jlim) dudzN = etaN * (x[j][i].w - x[j][i + 1].w) / dx;  in XMomentumResidual()
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
1119   jlim = grid->nj - 1;  in ViscosityField()
1126       if (i < ilim && j < jlim) {  in ViscosityField()
1294   PetscInt    i, j, mx, mz, ilim, jlim;  in FormFunctionLocal()  local
1302   jlim = mz - 1;  in FormFunctionLocal()
1326       } else if (j == jlim) {  in FormFunctionLocal()
1349       } else if (j == jlim) {  in FormFunctionLocal()
1377       } else if ((i == ilim || j == jlim) && ibound == BC_ANALYTIC) {  in FormFunctionLocal()
1400       } else if (j == jlim) {  in FormFunctionLocal()