read_nodes()

int read_nodes	(	int	igeom,
		const char *	geoname,
		struct nodes_geom *	nodesgeom
	)

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

  nodesgeom->filename= NULL;
  sprintf(varname, "%s.filename", geoname); CONFIG_GET_STRING(config, varname, nodesgeom->filename);
  sprintf(varname, "%s.ncntr", geoname); CONFIG_GET_INT(config, varname, geometries[igeom].ncntr);
  sprintf(varname, "%s.nnodb", geoname); CONFIG_GET_INT(config, varname, geometries[igeom].nnodb);
  sprintf(varname, "%s.nelmb", geoname); CONFIG_GET_INT(config, varname, geometries[igeom].nelmb);

  geometries[igeom].type = NODES;
  geometries[igeom].ptrgeom = (void*) nodesgeom;
  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, " file\n");
  logger(LOG_DEBUG, "------------------------\n");
  logger(LOG_DEBUG, "  - filename    = \"%s\"\n", nodesgeom->filename);


  return 0;


}