139     for (PetscInt fieldA = 0; fieldA < ctx->num_species; fieldA++) {  in LandauFormJacobian_Internal()
1189   ctx->num_species = nt;  in ProcessOptions()
1190 …for (ii = 0; ii < ctx->num_species; ii++) ctx->thermal_temps[ii] *= 1.1604525e7; /* convert to Kel…  in ProcessOptions()
1193 … == ctx->num_species - 1, ctx->comm, PETSC_ERR_ARG_WRONG, "num ion masses %" PetscInt_FMT " != num…  in ProcessOptions()
1196 … nm == ctx->num_species, ctx->comm, PETSC_ERR_ARG_WRONG, "wrong num n: %" PetscInt_FMT " != num sp…  in ProcessOptions()
1201 …c == ctx->num_species - 1, ctx->comm, PETSC_ERR_ARG_WRONG, "num charges %" PetscInt_FMT " != num s…  in ProcessOptions()
1209 …num_species == nt, ctx->comm, PETSC_ERR_ARG_WRONG, "-dm_landau_num_species_grid: sum %" PetscInt_F…  in ProcessOptions()
1212     if (ctx->num_species > LANDAU_MAX_GRIDS) {  in ProcessOptions()
1214       num_species_grid[1] = ctx->num_species - 1;  in ProcessOptions()
1217       ctx->num_grids = ctx->num_species;  in ProcessOptions()
1222 …PetscCheck(ctx->species_offset[ctx->num_grids] == ctx->num_species, ctx->comm, PETSC_ERR_ARG_WRONG…  in ProcessOptions()
1223              ctx->num_species);  in ProcessOptions()
1240 …PetscCheck(non_dim_grid >= 0 && non_dim_grid < ctx->num_species, ctx->comm, PETSC_ERR_ARG_WRONG, "…  in ProcessOptions()
1343 …for (ii = ctx->num_species; ii < LANDAU_MAX_SPECIES; ii++) ctx->masses[ii] = ctx->thermal_temps[ii…  in ProcessOptions()
1347 …le)ctx->masses[0], (double)(ctx->masses[1] / pmassunit), (double)(ctx->num_species > 2 ? ctx->mass…  in ProcessOptions()
1348 …harges[0], (double)(-ctx->charges[1] / ctx->charges[0]), (double)(ctx->num_species > 2 ? -ctx->cha…  in ProcessOptions()
1349 … %10.3e %10.3e\n", (double)ctx->n[0], (double)ctx->n[1], (double)(ctx->num_species > 2 ? ctx->n[2]…  in ProcessOptions()
1351 …(double)ctx->thermal_temps[1], (double)((ctx->num_species > 2) ? ctx->thermal_temps[2] : 0), non_d…  in ProcessOptions()
1481     const PetscInt  *plex_batch = NULL, elMatSz = Nb * Nb * ctx->num_species * ctx->num_species;  in CreateStaticData()
1713                             ctx->num_species, Nb, dim, N, (PetscInt)info.nz_used));  in CreateStaticData()
1830 …PetscCall(PetscMalloc4(ctx->num_species, &nu_alpha_p, ctx->num_species, &nu_beta_p, ctx->num_speci…  in CreateStaticData()
1831       for (PetscInt ii = 0; ii < ctx->num_species; ii++) {  in CreateStaticData()
2277   for (PetscInt ii = 0; ii < ctx->num_species; ii++) PetscCall(PetscFEDestroy(&ctx->fe[ii]));  in DMPlexLandauDestroyVelocitySpace()
2521       if (ctx->num_species > 1) PetscCall(PetscPrintf(PETSC_COMM_WORLD, "\n"));  in DMPlexLandauPrintNorms()
2528   if (ctx->num_species > 1) {  in DMPlexLandauPrintNorms()