556 …PetscReal       lo[3] = {-ctx->radius[grid], -ctx->radius[grid], -ctx->radius[grid]}, hi[3] = {ctx…  in go()  local
563             if (dim == 2) lo[0] = 0; // Landau coordinate (r,z)  in go()
570             hp[0] = (hi[0] - lo[0]) / Npi;  in go()
571             hp[1] = (hi[1] - lo[1]) / Npj;  in go()
572             hp[2] = (hi[2] - lo[2]) / Npk;  in go()
574 …, " lo = %14.7e, hi = %14.7e; hp = %14.7e, %14.7e; kT_m = %g; \n", (double)lo[1], (double)hi[1], (…  in go()
583                     xx_t[grid][tid][pp] = lo[0] + hp[0] / 2.0 + pi * hp[0];  in go()
584                     yy_t[grid][tid][pp] = lo[1] + hp[1] / 2.0 + pj * hp[1];  in go()
585                     if (dim == 3) zz_t[grid][tid][pp] = lo[2] + hp[2] / 2.0 + pk * hp[2];  in go()
607 …           if (dim == 2 && xx_t[grid][tid][pp] < lo[0]) xx_t[grid][tid][pp] = -xx_t[grid][tid][pp];  in go()
608                       if (dim == 3) zz_t[grid][tid][pp] = lo[2] + hp[2] / 2.0 + pk * hp[2];  in go()
612 …    while (zz_t[grid][tid][pp] >= hi[2] || zz_t[grid][tid][pp] <= lo[2]) zz_t[grid][tid][pp] *= .9;  in go()
614 …    while (xx_t[grid][tid][pp] >= hi[0] || xx_t[grid][tid][pp] <= lo[0]) xx_t[grid][tid][pp] *= .9;  in go()
615 …    while (yy_t[grid][tid][pp] >= hi[1] || yy_t[grid][tid][pp] <= lo[1]) yy_t[grid][tid][pp] *= .9;  in go()
655               yy_t[grid][tid][pp++] = lo[1] + 5.e-7;  in go()
659               xx_t[grid][tid][p0 + 0] = lo[0];  in go()
661               yy_t[grid][tid][p0 + 2] = lo[1];  in go()
663               zz_t[grid][tid][p0 + 4] = lo[2];  in go()
859 …PetscReal lo[3] = {-ctx->radius[grid], -ctx->radius[grid], -ctx->radius[grid]}, hi[3] = {ctx->radi…  in go()  local
866               if (dim == 2) lo[0] = 0;  in go()
877               hp[0] = (hi[0] - lo[0]) / Npi;  in go()
878               hp[1] = (hi[1] - lo[1]) / Npj;  in go()
879               hp[2] = (hi[2] - lo[2]) / Npk;  in go()
886                     xx_t[grid][tid][pp] = lo[0] + hp[0] / 2.0 + pi * hp[0];  in go()
887                     yy_t[grid][tid][pp] = lo[1] + hp[1] / 2.0 + pj * hp[1];  in go()
888                     if (dim == 3) zz_t[grid][tid][pp] = lo[2] + hp[2] / 2.0 + pk * hp[2];  in go()
902                 yy_t[grid][tid][pp++] = lo[1] + 5.e-7;  in go()
906                 xx_t[grid][tid][p0 + 0] = lo[0];  in go()
908                 yy_t[grid][tid][p0 + 2] = lo[1];  in go()
910                 zz_t[grid][tid][p0 + 4] = lo[2];  in go()