vrmlout()

void vrmlout

(

)

Produces the VRML files for visualization and handling of the geometry. Uses VRML 2.0 format (http://www.web3d.org)

              {
  int in1,in2,in3,in4,i;
  FILE *f;
  char fname[MAX_PATH];

  logger(LOG_DEBUG,"About to produce VRML files \n");
  
  if (0 != rank) return; // Does the job only on master process

  get_filename(fname,"%s-surface.wrl",rundetails->title);  // Domain Boundary
  f = fopen(fname,"w");
  fprintf(f,"#VRML V2.0 utf8 \n");     
  fprintf(f,"Transform { \n translation 0.0 0.0 0.0 \n children[ \n Shape { \n");
  fprintf(f,"appearance Appearance{ \n material Material{ \n diffuseColor 1.0 0.0 0.0 \n");
  fprintf(f,"emissiveColor 1.0 0.0 0.0 \n shininess 0.8 }} \n");
  fprintf(f,"geometry IndexedFaceSet { \n coord Coordinate { point [ \n");
  for(i=0;i<modgeominfo->nnodb;i++){
    fprintf(f,"  %f   %f   %f \n",(double)geometry->nodes[i].x,(double)geometry->nodes[i].y,(double)geometry->nodes[i].z);
  }
  fprintf(f,"] } coordIndex [ \n");
  for(i=0;i<modgeominfo->nelmb;i++){
    in1 = geometry->elements[i].idx[0];
    in2 = geometry->elements[i].idx[1];
    in3 = geometry->elements[i].idx[2];
    in4 = geometry->elements[i].idx[3];
        fprintf(f,"%d   %d   %d   %d   %d -1\n",in1,in2,in3,in4,in1);
  }
  fprintf(f," ] } } ] }\n");
  fclose(f);

  get_filename(fname,"%s-controls.wrl",rundetails->title);  // Domain Boundary
  f = fopen(fname,"w");
  fprintf(f,"#VRML V2.0 utf8 \n");     
  fprintf(f,"Transform { \n translation 0.0 0.0 0.0 \n children[ \n Shape { \n");
  fprintf(f,"appearance Appearance{ \n material Material{ \n diffuseColor 1.0 1.0 0.0 \n");
  fprintf(f,"emissiveColor 1.0 1.0 0.0 \n shininess 0.8 }} \n");
  fprintf(f,"geometry PointSet { \n coord Coordinate { point [ \n");
  for(i=0;i<modgeominfo->ncntr;i++){
    fprintf(f,"  %f   %f   %f \n",(double)geometry->cntr[i].x,(double)geometry->cntr[i].y,(double)geometry->cntr[i].z);
  }
  /*fprintf(f,"] } normal Normal { vector [ \n");
  for(i=0;i<modgeominfo->nelmb;i++){
          fprintf(f,"%f   %f   %f \n",geometry->elements[i].n.x,geometry->elements[i].n.y,geometry->elements[i].n.z);
  }*/
  fprintf(f,"] } \n");
  fprintf(f,"} } ] }\n");
  fclose(f);

  if(0 != modgeominfo->nchief){  
          get_filename(fname,"%s-chief.wrl",rundetails->title);  // Domain Boundary
          f = fopen(fname,"w");
          fprintf(f,"#VRML V2.0 utf8 \n");     
          fprintf(f,"Transform { \n translation 0.0 0.0 0.0 \n children[ \n Shape { \n");
          fprintf(f,"appearance Appearance{ \n material Material{ \n diffuseColor 0.7 0.7 0.7 \n");
          fprintf(f,"emissiveColor 0.7 0.7 0.7 \n shininess 0.8 }} \n");
          fprintf(f,"geometry PointSet { \n coord Coordinate { point [ \n");
          for(i=0;i<modgeominfo->nchief;i++){
            fprintf(f,"  %f   %f   %f \n",(double)geometry->chief[i].x,(double)geometry->chief[i].y,(double)geometry->chief[i].z);
          }
          fprintf(f,"] } \n");
          fprintf(f,"} } ] }\n");
          fclose(f);
  }
  if(0 != modgeominfo->nmics){  
          get_filename(fname,"%s-mics.wrl",rundetails->title);  // Domain Boundary
          f = fopen(fname,"w");
          fprintf(f,"#VRML V2.0 utf8 \n");     
          fprintf(f,"Transform { \n translation 0.0 0.0 0.0 \n children[ \n Shape { \n");
          fprintf(f,"appearance Appearance{ \n material Material{ \n diffuseColor 0.0 0.1 0.1 \n");
          fprintf(f,"emissiveColor 0.0 0.1 0.1 \n shininess 0.8 }} \n");
          fprintf(f,"geometry PointSet { \n coord Coordinate { point [ \n");
          for(i=0;i<modgeominfo->nmics;i++){
            fprintf(f,"  %f   %f   %f \n",(double)geometry->mics[i].x,(double)geometry->mics[i].y,(double)geometry->mics[i].z);
          }
          fprintf(f,"] } \n");
          fprintf(f,"} } ] }\n");
          fclose(f);
  }
  

  logger(LOG_DEBUG,"VRML files produced\n");
}