read_mics()

int read_mics

(

)

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

  fmics = fopen(fileinfo->mics_file, "r");
  havemicsmesh=0;
  if (NULL == fmics) {
    logger(LOG_ERROR, MSG_GEOM_FILE_NOTFOUND, fileinfo->mics_file);
    return ACOUSTO_ERROR_GEOMETRY;
  }
  fgets(line, 100, fmics);
  if(0 == strcmp("$MeshFormat\n", line)){
    havemicsmesh = 1;
    logger(LOG_INFO,"Mics file appears to be a GMSH FILE\n");
    logger(LOG_INFO,"    Proceed with reading of number of mics and mics elements\n");
    nelcount = 0;
    while (fgets(line, 100, fmics) != NULL) {
      if (0 == strcmp("$Nodes\n", line)) {
        fgets(line, 100, fmics);
        sscanf(line, "%d\n", &ntmp);
        modgeominfo->nmics = ntmp;
        logger(LOG_INFO, "MICS GMSH FILE -> number of microphones %d\n", modgeominfo->nmics);
      }
      if (0 == strcmp("$Elements\n", line)) {
        fgets(line, 100, fmics);
        sscanf(line, "%d\n", &netmp);
        for (j = 0; j < netmp; j++) {
          fgets(line, 100, fmics);
          sscanf(line, "%d %d\n", &ntmp, &eltype);
          if (2 == eltype || 3 == eltype || 4 == eltype) {
            nelcount++;
          }
        }
        modgeominfo->nemics = nelcount;
        logger(LOG_INFO,
              "MICS GMSH FILE -> number of read triangular, quadrilateral and tetrahedra elements %d\n",
              nelcount);
      }
    }
    fclose(fmics);
  }else{
    logger(LOG_INFO,"Mics file appears to be in AcouSTO native format\n");
    logger(LOG_INFO,"    Number of mics is read from config file\n");
    CONFIG_GET_INT(config, "modgeom.nmics", modgeominfo->nmics);
  }


  if (modgeominfo->nmics < size) {
    logger(LOG_WARN, MSG_CFG_NOTENOUGHMICS, modgeominfo->nmics, size);
    modgeominfo->nmics = size;
  }


  return 0;


}