read_plate()

int read_plate	(	int	igeom,
		const char *	geoname,
		struct plate_geom *	plategeom
	)

                                                                           {
  int ltmp;
  double dtmp;
  char varname[MAX_PATH]; // temp var to read dynamic geometry &configuration
  int rchk, mchk;
  sprintf(varname, "%s.l1", geoname); CONFIG_GET_FLOAT(config, varname, plategeom->lx);
  sprintf(varname, "%s.l2", geoname); CONFIG_GET_FLOAT(config, varname, plategeom->ly);
  sprintf(varname, "%s.n1", geoname); CONFIG_GET_INT(config, varname, plategeom->nx);
  sprintf(varname, "%s.n2", geoname); CONFIG_GET_INT(config, varname, plategeom->ny);

  geometries[igeom].type = PLATE;
  geometries[igeom].ptrgeom = (void*) plategeom;
  geometries[igeom].nelmb = (plategeom->nx) * (plategeom->ny);
  geometries[igeom].nnodb = (plategeom->nx + 1) * (plategeom->ny + 1);
  sprintf(varname, "%s.ncntr", geoname);
  ltmp = 0;
  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, " plate\n");
  logger(LOG_DEBUG, "------------------------\n");
  logger(LOG_DEBUG, "  - n1          = %d\n", plategeom->nx);
  logger(LOG_DEBUG, "  - n2          = %d\n", plategeom->ny);
  logger(LOG_DEBUG, "  - l1          = %f\n", plategeom->ly);
  logger(LOG_DEBUG, "  - l2          = %f\n", plategeom->ly);


 
  return 0;


}