read_sphere()

int read_sphere	(	int	igeom,
		const char *	geoname,
		struct sphere_geom *	sphere
	)

                                                                          {
  int ltmp;
  double dtmp;
  char varname[MAX_PATH]; // temp var to read dynamic geometry &configuration
  int rchk, mchk;

  logger(LOG_DEBUG, "Sphere\n");

  sprintf(varname, "%s.radius", geoname); CONFIG_GET_FLOAT(config, varname, sphere->radius);
  sprintf(varname, "%s.segments", geoname); CONFIG_GET_INT(config, varname, sphere->segments);
  sprintf(varname, "%s.rings", geoname); CONFIG_GET_INT(config, varname, sphere->rings);

  geometries[igeom].type = SPHERE;
  geometries[igeom].ptrgeom = (void*) sphere;
  geometries[igeom].nelmb = (sphere->segments) * (sphere->rings);
  geometries[igeom].nnodb = (sphere->segments + 1) * (sphere->rings + 1);
  sprintf(varname, "%s.ncntr", geoname);
  CONFIG_GET_INT(config, varname, geometries[igeom].ncntr);
  if (0 == geometries[igeom].ncntr) {
    if (0 == runinfo->ksymmi) {
      geometries[igeom].ncntr = geometries[igeom].nelmb;
    } else {
      geometries[igeom].ncntr = (int) (geometries[igeom].nelmb / (runinfo->nsymm));
    }
  } else {
    rchk = geometries[igeom].nelmb / geometries[igeom].ncntr;
    mchk = geometries[igeom].nelmb % geometries[igeom].ncntr;
    if (rchk != runinfo->nsymm)
      logger(LOG_ERROR, "Number of elements and centers not consistent with symmetry level required: nelmb=%d, ncntr=%d, nsymm=%d \n", geometries[igeom].nelmb, geometries[igeom].ncntr, runinfo->nsymm);
    if (mchk != 0)
      logger(LOG_ERROR, "Number of elements is not a multiple integer of ncntr: nelmb(mod)ncntr=%d \n", mchk);
    if ((rchk != runinfo->nsymm) || (mchk != 0)) {
      return -1;
    }
  }
  logger(LOG_DEBUG, "------------------------\n");
  logger(LOG_DEBUG, " sphere\n");
  logger(LOG_DEBUG, "------------------------\n");
  logger(LOG_DEBUG, "  - radius      = %f\n", sphere->radius);
  logger(LOG_DEBUG, "  - rings       = %d\n", sphere->rings);
  logger(LOG_DEBUG, "  - segments    = %d\n", sphere->segments);

  return 0;

}