96   PetscInt        corner, ni, nj, jlid, jfault, inose;  member
313   x  = (i - grid->jlid) * grid->dx;  in HorizVelocity()
314   z  = (j - grid->jlid - 0.5) * grid->dz;  in HorizVelocity()
330   x  = (i - grid->jlid - 0.5) * grid->dx;  in VertVelocity()
331   z  = (j - grid->jlid) * grid->dz;  in VertVelocity()
346   x  = (i - grid->jlid - 0.5) * grid->dx;  in Pressure()
347   z  = (j - grid->jlid - 0.5) * grid->dz;  in Pressure()
368     if (j <= grid->jlid) return EPS_ZERO;  in CalcSecInv()
383     if (j == grid->jlid + 1) uS = 0.0;  in CalcSecInv()
387     if (j < grid->jlid) return EPS_ZERO;  in CalcSecInv()
401     if (j == grid->jlid) {  in CalcSecInv()
588   PetscInt    ilim = grid->ni - 1, jlim = grid->nj - 1, jlid = grid->jlid;  in EnergyResidual()  local
602   if (j <= jlid && i >= j) {  in EnergyResidual()
649   if (j <= grid->jlid || i < j || i == ilim || j == jlim) return EPS_ZERO;  in ShearStress()
661     if (j == grid->jlid + 1) uS = 0.0;  in ShearStress()
683   if (i < j || j <= grid->jlid) return EPS_ZERO;  in XNormalStress()
723   if (i < j || j <= grid->jlid) return EPS_ZERO;  in ZNormalStress()
877   grid->jlid          = (PetscInt)(param->lid_depth / grid->dz);   /* gridcells */  in SetParams()
879   param->lid_depth    = grid->jlid * grid->dz;                     /* dim'less */  in SetParams()
881   grid->corner        = grid->jlid + 1;                            /* gridcells */  in SetParams()
908 …jlid = %3" PetscInt_FMT "              jfault = %3" PetscInt_FMT " \n", grid->jlid, grid->jfault));  in ReportParams()
982       else if (j <= grid->jlid) x[j][i].u = 0.0;  in Initialize()
986       else if (j <= grid->jlid) x[j][i].w = 0.0;  in Initialize()
989       if (i < j || j <= grid->jlid) x[j][i].p = 0.0;  in Initialize()
1044       PetscCall(VecSetValue(pars, 8, (PetscScalar)grid->jlid, INSERT_VALUES));  in DoOutput()
1315 …  else if (j <= grid->jlid || (j < grid->corner + grid->inose && i < grid->corner + grid->inose)) {  in FormFunctionLocal()
1345 …} else if (j <= grid->jlid || (j < grid->corner + grid->inose && i < grid->corner + grid->inose)) {  in FormFunctionLocal()
1373 …if (i < j || j <= grid->jlid || (j < grid->corner + grid->inose && i < grid->corner + grid->inose)…  in FormFunctionLocal()