IBWriteBodySTL.c

Go to the documentation of this file.

#include <stdio.h>
#include <stdlib.h>
#include <basic.h>
#include <mpivars.h>
#include <immersedboundaries.h>


int IBWriteBodySTL(
                    Body3D *body,     
                    char   *filename, 
                    void   *m,        
                    int    rank,      
                    int    *stat      
                  )
{
  MPIVariables  *mpi = (MPIVariables*) m;

  if (mpi->rank == rank) {
    FILE *out;
    out = fopen(filename,"w");
    if (!out) {
      fprintf(stderr,"Error in IBWriteBodySTL(): Unable to open file ");
      fprintf(stderr,"%s for writing on rank %d.\n",filename,mpi->rank);
      *stat = 1;
    } else {
      fprintf(out,"solid TEST\n");
      int n;
      for (n = 0; n < body->nfacets; n++) {
        fprintf(out,"  facet normal %1.16e %1.16e %1.16e\n",
                body->surface[n].nx,body->surface[n].ny,body->surface[n].nz);
        fprintf(out,"    outer loop\n");
        fprintf(out,"      vertex %1.16e %1.16e %1.16e\n",
                body->surface[n].x1,body->surface[n].y1,body->surface[n].z1);
        fprintf(out,"      vertex %1.16e %1.16e %1.16e\n",
                body->surface[n].x2,body->surface[n].y2,body->surface[n].z2);
        fprintf(out,"      vertex %1.16e %1.16e %1.16e\n",
                body->surface[n].x3,body->surface[n].y3,body->surface[n].z3);
        fprintf(out,"    endloop\n");
        fprintf(out,"  endfacet\n");
      }
      fprintf(out,"endsolid TEST\n");
      *stat = 0;
    }
  }
  return(0);
}