read_cylinder()

int read_cylinder	(	int	igeom,
		const char *	geoname,
		struct cylinder_geom *	cylinder
	)

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

  logger(LOG_DEBUG, "Cylinder\n");
  cylinder = malloc(sizeof(struct cylinder_geom));
  sprintf(varname, "%s.start_radius", geoname); CONFIG_GET_FLOAT(config, varname, cylinder->start_radius);              // gp mod
  sprintf(varname, "%s.end_radius", geoname); CONFIG_GET_FLOAT(config, varname, cylinder->end_radius);                  // gp mod
  sprintf(varname, "%s.rings", geoname); CONFIG_GET_INT(config, varname, cylinder->rings);
  sprintf(varname, "%s.segments", geoname); CONFIG_GET_INT(config, varname, cylinder->segments);
  sprintf(varname, "%s.length", geoname); CONFIG_GET_FLOAT(config, varname, cylinder->length);

  geometries[igeom].type = CYLINDER;
  geometries[igeom].ptrgeom = (void*) cylinder;
  geometries[igeom].nelmb = (cylinder->segments) * (cylinder->rings);
  geometries[igeom].nnodb = (cylinder->segments + 1) * (cylinder->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, " cylinder\n");
  logger(LOG_DEBUG, "------------------------\n");
  logger(LOG_DEBUG, "  - first radius      = %f\n", cylinder->start_radius);            // gp mod
  logger(LOG_DEBUG, "  - last radius      = %f\n", cylinder->end_radius);                       // gp mod
  logger(LOG_DEBUG, "  - length      = %f\n", cylinder->length);
  logger(LOG_DEBUG, "  - rings       = %d\n", cylinder->rings);
  logger(LOG_DEBUG, "  - segments    = %d\n", cylinder->segments);


  return 0;


}