build_nodes()

build_nodes

(

struct geometry *

geom

)

Build nodes from external file

Parameters

[in,out]

geom

structure

                                      {
  int inodb,ivbase;
  long double x,y,z;
  int i,j1,j2,j3,j4;
  char line[100];
  FILE* fnodes;
  struct nodes_geom    *nodesgeom;


  nodesgeom  = (struct nodes_geom*) geom->ptrgeom;
  logger(LOG_DEBUG,"loading file %s\n",nodesgeom->filename);
  /* opening file */
  fnodes = fopen(nodesgeom->filename,"r");
  if(NULL == fnodes){
    logger(LOG_ERROR,MSG_GEOM_FILE_NOTFOUND,nodesgeom->filename);
    return ACOUSTO_ERROR_GEOMETRY;
  }


  /* Reading nodes */
  for(inodb=0;inodb<geom->nnodb;inodb++){
    fgets(line,100,fnodes);
    sscanf(line,VECTOR_SCANF,&x,&y,&z);
    geometry->nodes[inodb + _g_inodb].x = x;
    geometry->nodes[inodb + _g_inodb].y = y;
    geometry->nodes[inodb + _g_inodb].z = z;
  }


  /* Reading topology */
  for(i=1;i<=geom->nelmb;i++){
    ivbase = (i-1)*4;
    fgets(line,100,fnodes);
    sscanf(line," %d %d %d %d\n",&j1,&j2,&j3,&j4);
    geometry->jnodb[ivbase   + _g_ielmb] = j1-1;
    geometry->jnodb[ivbase+1 + _g_ielmb] = j2-1;
    geometry->jnodb[ivbase+2 + _g_ielmb] = j3-1;
    geometry->jnodb[ivbase+3 + _g_ielmb] = j4-1;
  }

  geom_fill(geom); // filling geometry

  //adding to global vars
  _g_icntr += geom->ncntr; 
  _g_inodb += geom->nnodb;
  _g_ielmb += geom->nelmb;

  if(nodesgeom->filename) logger(LOG_DEBUG,"file %s read\n",nodesgeom->filename);
  fclose(fnodes);

  return 1;

}