read_gmsh()

int read_gmsh	(	int	igeom,
		const char *	geoname,
		struct gmsh_geom *	gmshgeom
	)

                                                                        {
  double dtmp;
  char varname[MAX_PATH]; // temp var to read dynamic geometry &configuration

  int j;
  int eltype;
  int netmp;
  int ntmp;
  unsigned long nelcount = 0;
  FILE* fnodes;
  char line[MAX_PATH];

  gmshgeom->filename = NULL;
  sprintf(varname, "%s.filename", geoname); CONFIG_GET_STRING(config, varname, gmshgeom->filename);
  sprintf(varname, "%s.symvec.x", geoname); CONFIG_GET_FLOAT(config, varname, gmshgeom->symvec.x);
  sprintf(varname, "%s.symvec.y", geoname); CONFIG_GET_FLOAT(config, varname, gmshgeom->symvec.y);
  sprintf(varname, "%s.symvec.z", geoname); CONFIG_GET_FLOAT(config, varname, gmshgeom->symvec.z);
  sprintf(varname, "%s.refpoint.x", geoname); CONFIG_GET_FLOAT(config, varname, gmshgeom->refpoint.x);
  sprintf(varname, "%s.refpoint.y", geoname); CONFIG_GET_FLOAT(config, varname, gmshgeom->refpoint.y);
  sprintf(varname, "%s.refpoint.z", geoname); CONFIG_GET_FLOAT(config, varname, gmshgeom->refpoint.z);

  if (runinfo->ksymmi == 2 || runinfo->ksymmi == 3) {
    logger(LOG_ERROR, "Symmetry level ksymmi=%d. Acousto supports only 1 symmetry plane (ksymmi=1) or axialsymmetry (ksymmi>3) with GMSH files \n", runinfo->ksymmi);
    return -1;
  }

  /* Reading .msh file for number of nodes and elements */

  fnodes = fopen(gmshgeom->filename, "r");
  if (NULL == fnodes) {
    logger(LOG_ERROR, MSG_GEOM_FILE_NOTFOUND, gmshgeom->filename);
    return ACOUSTO_ERROR_GEOMETRY;
  }

  nelcount = 0;
  while (fgets(line, 100, fnodes) != NULL) {
    if (0 == strcmp("$Nodes\n", line)) {
      fgets(line, 100, fnodes);
      sscanf(line, "%d\n", &ntmp);
      geometries[igeom].nnodb = runinfo->nsymm * ntmp;
      logger(LOG_INFO, "GMSH FILE -> number of read nodes %d\n", ntmp);
      logger(LOG_INFO, "GMSH FILE -> number of nodes after symmetrization %d\n", geometries[igeom].nnodb);
    }
    if (0 == strcmp("$Elements\n", line)) {
      fgets(line, 100, fnodes);
      sscanf(line, "%d\n", &netmp);
      for (j = 0; j < netmp; j++) {
        fgets(line, 100, fnodes);
        sscanf(line, "%d %d\n", &ntmp, &eltype);
        if (2 == eltype || 3 == eltype) {
          nelcount++;
        }
      }
      geometries[igeom].nelmb = runinfo->nsymm * nelcount;
      logger(LOG_INFO, "GMSH FILE -> number of read triangular and quadrilateral elements %d\n", nelcount);
      logger(LOG_INFO, "GMSH FILE -> number of triangular and quadrilateral elements after symmetrization %d\n", geometries[igeom].nelmb);
    }
  }
  fclose(fnodes);

  geometries[igeom].type = GMSH;
  geometries[igeom].ptrgeom = (void*) gmshgeom;
  geometries[igeom].ncntr = (int) (geometries[igeom].nelmb / (runinfo->nsymm));
  logger(LOG_INFO, "GMSH FILE -> number of collocation points calculated %d\n", geometries[igeom].ncntr);

  logger(LOG_DEBUG, "------------------------\n");
  logger(LOG_DEBUG, " file\n");
  logger(LOG_DEBUG, "------------------------\n");
  logger(LOG_DEBUG, "  - filename    = \"%s\"\n", gmshgeom->filename);


  return 0;


}