vtkout_unc()

void vtkout_unc

(

)

Routine vtkout_unc(icounter,CFREQ). Produces the file 'runname-unormals.vtk', which can bee used to visualize the unit normal vectors on the boundary. Produces the file 'runname-body.raw', which can bee used to import the panels' set into Blender (or other modelers) to cross-check the geometry properties. Arguments: none; Data format is 'DATASET UNSTRUCTURED_GRID'.

                 {

  int i;
  FILE *f;
  char fname[MAX_PATH];
  double xc,yc,zc;

  if (0 != rank) return;
  
  get_filename(fname,"%s-unormals.vtk",rundetails->title);  
  f = fopen(fname,"w");

  fprintf(f,"# vtk DataFile Version 2.0\n");
  fprintf(f,"Unit Normal Vectors\n");
  fprintf(f,"ASCII\n");

  fprintf(f,"DATASET UNSTRUCTURED_GRID\n");
  fprintf(f,"POINTS %d float \n",modgeominfo->nelmb);
  for(i=0;i<modgeominfo->nelmb;i++){
          xc = 0.25*(geometry->elements[i].x[0].x+geometry->elements[i].x[1].x+geometry->elements[i].x[2].x+geometry->elements[i].x[3].x);
          yc = 0.25*(geometry->elements[i].x[0].y+geometry->elements[i].x[1].y+geometry->elements[i].x[2].y+geometry->elements[i].x[3].y);
          zc = 0.25*(geometry->elements[i].x[0].z+geometry->elements[i].x[1].z+geometry->elements[i].x[2].z+geometry->elements[i].x[3].z);
          fprintf(f,"  %f   %f   %f \n",xc,yc,zc);
  }
  fprintf(f,"POINT_DATA %d\n", modgeominfo->nelmb);
  fprintf(f,"VECTORS unit_normal float\n");
  for(i=0;i<modgeominfo->nelmb;i++){
            fprintf(f,"%f   %f   %f \n",(double)geometry->elements[i].n.x,(double)geometry->elements[i].n.y,(double)geometry->elements[i].n.z);
  }

  fclose(f);

  // Raw data output for Blender import
  get_filename(fname,"%s-body.raw",rundetails->title);  
  f = fopen(fname,"w");
  for(i=0;i<modgeominfo->nelmb;i++){
          fprintf(f,"  %f   %f   %f  %f   %f   %f  %f   %f   %f  %f   %f   %f\n",(double)geometry->elements[i].x[0].x,
                                                                           (double)geometry->elements[i].x[0].y,
                                                                           (double)geometry->elements[i].x[0].z,
                                                                           (double)geometry->elements[i].x[1].x,
                                                                           (double)geometry->elements[i].x[1].y,
                                                                           (double)geometry->elements[i].x[1].z,
                                                                           (double)geometry->elements[i].x[2].x,
                                                                           (double)geometry->elements[i].x[2].y,
                                                                           (double)geometry->elements[i].x[2].z,
                                                                           (double)geometry->elements[i].x[3].x,
                                                                           (double)geometry->elements[i].x[3].y,
                                                                           (double)geometry->elements[i].x[3].z);
  }

  fclose(f);

}