read_acoustic_files()

int read_acoustic_files

(

)

Reads acoustic solution files.

                         {

  FILE*  file_sources  = NULL ;
  FILE*  file_planw    = NULL;
  FILE*  file_dipoles  = NULL ;

  char line[200];
  int i;
  long double d1,d2,d3,d4,d5,d6,d7,d8,d9;

  int nfreq;

  int mpi_ncntr;
  int nphiu;

  nfreq = modsolinfo->nsig * modsolinfo->nome;
  mpi_ncntr  = modgeominfo->ncntr + modgeominfo->nchief;
  nphiu = modgeominfo->ncntr;


  /* open files */
  if(modsolinfo->nsourc > 0){
   logger(LOG_DEBUG," reading file %s\n",fileinfo->sources_file);
   file_sources = fopen(fileinfo->sources_file,"r");
   if(NULL == file_sources ){
     perror("Coudln't find SOURCES file ");
      return 0;
   }
  }

  if(modsolinfo->nplanw > 0){
   logger(LOG_DEBUG," reading file %s\n",fileinfo->planw_file);
   file_planw = fopen(fileinfo->planw_file,"r");
   if(NULL == file_planw){
     perror("Coudln't find PLAN WAVES file ");
      return 0;
   }
  }


  if(modsolinfo->ndipoles > 0){
    logger(LOG_DEBUG," reading file %s\n",fileinfo->dipoles_file);
    file_dipoles = fopen(fileinfo->dipoles_file,"r");
    if(NULL == file_dipoles && modsolinfo->ndipoles > 0){
      perror("Coudln't find DIPOLES file ");
       return 0;
    }
  }
 /* Reading */

 /* SOURCES */
 for(i=0;i<modsolinfo->nsourc;i++){
   fgets(line,100,file_sources);
   sscanf(line,ACPOINT_SCANF,&d1,&d2,&d3,&d4,&d5);
   solution->sources[i].pos.x = d1;
   solution->sources[i].pos.y = d2;
   solution->sources[i].pos.z = d3;
   solution->sources[i].amplitude = d4 + I*d5;
   _DEBUG_ACOUPOINT("SOURCES",rank,solution->sources[i]);
   LOG_ACPOINT(LOG_DEBUG,i,solution->sources[i]);
 }

 /* PLANE WAVES */
 for(i=0;i<modsolinfo->nplanw;i++){
   fgets(line,100,file_planw);
   sscanf(line,ACPOINT_SCANF,&d1,&d2,&d3,&d4,&d5);
   solution->planw[i].pos.x = d1;
   solution->planw[i].pos.y = d2;
   solution->planw[i].pos.z = d3;
   solution->planw[i].amplitude = d4 + I*d5;
   _DEBUG_ACOUPOINT("PLANW",rank,solution->planw[i]);
   LOG_ACPOINT(LOG_DEBUG,i,solution->planw[i]);
 }


 /* DIPOLES */
 for(i=0;i<modsolinfo->ndipoles;i++){
   fgets(line,100,file_dipoles);
   sscanf(line,DIPOLE_SCANF,&d1,&d2,&d3,
                            &d4,&d5,&d6,
                            &d7,&d8,
                            &d9);
   solution->dipoles[i].pos.x = d1;
   solution->dipoles[i].pos.y = d2;
   solution->dipoles[i].pos.z = d3;
   solution->dipoles[i].dir.x = d4;
   solution->dipoles[i].dir.y = d5;
   solution->dipoles[i].dir.z = d6;
   solution->dipoles[i].amplitude = d7 + I*d8;
   solution->dipoles[i].d = d9;
   _DEBUG_DIPOLE("DIPOLES",i,solution->dipoles[i]);
   LOG_DIPOLE(LOG_DEBUG,i,solution->dipoles[i]);
 }

 if(modsolinfo->nsourc  > 0 && NULL != file_sources ) fclose(file_sources);
 if(modsolinfo->ndipoles> 0 && NULL != file_dipoles ) fclose(file_dipoles);
 if(modsolinfo->nplanw  > 0 && NULL != file_planw   ) fclose(file_planw);


 if(modsolinfo->nsourc > 0 )  logger(LOG_DEBUG,"file %s read\n",fileinfo->sources_file);
 if(modsolinfo->ndipoles > 0) logger(LOG_DEBUG,"file %s read\n",fileinfo->dipoles_file);
 if(modsolinfo->nplanw > 0 )  logger(LOG_DEBUG,"file %s read\n",fileinfo->planw_file);
 return 1;
}