source: issm/branches/trunk-jpl-damage/src/c/shared/Exp/DomainOutlineWrite.cpp@ 11427

/*!\file:  DomainOutlineWrite.cpp
 * \brief DomainOutlineWrite.c: write the vertex coordinates defined in a domain 
 * outline from Argus (.exp file). The first contour in the file is for 
 * the outside domain outline. The following contours represent holes in
 * the domain.
 */

#include <stdio.h>
#include "../Alloc/alloc.h"
#include "../../include/include.h"
#include "../Exceptions/exceptions.h"

int DomainOutlineWrite(int nprof,int* profnvertices,double** pprofx,double** pprofy,bool* closed,char* domainname,bool whole=true){

        
        /*Error management: */
        int noerr=1;
        int i,counter;

        /*I/O: */
        FILE* fid=NULL;

        /*input: */
//      int nprof; //number of profiles in the domainname file
//      int* profnvertices=NULL; //array holding the number of vertices for the nprof profiles
//      double** pprofx=NULL; //array of profiles x coordinates
//      double** pprofy=NULL; //array of profiles y coordinates
//      bool* closed=NULL; //array holding closed flags for the nprof profiles

        /*open domain outline file for writing: */
        if ((fid=fopen(domainname,"w"))==NULL){
                _error_("%s%s","could not open domain file ",domainname); 
                noerr=0; goto cleanupandreturn;
        }

        /*Start writing profiles: */
        for(counter=0;counter<nprof;counter++){

                /*Write header: */
                fprintf(fid,"%s %s\n","##","Name:");
                fprintf(fid,"%s %s\n","##","Icon:0");
                fprintf(fid,"%s %s %s %s\n","#","Points","Count","Value");
                
                /*Write number of profile vertices: */
                if(closed[counter] && !whole)
                        fprintf(fid,"%u %s\n",profnvertices[counter]+1,"1.");
                else
                        fprintf(fid,"%u %s\n",profnvertices[counter]  ,"1.");
        
                /*Write next line: */
                fprintf(fid,"%s %s %s %s %s\n","#","X","pos","Y","pos");

                /*Write vertices: */
                for (i=0;i<profnvertices[counter];i++){
                        fprintf(fid,"%lf\t%lf\n",pprofx[counter][i],pprofy[counter][i]);
                }

                /*Now check that we are dealing with open contours: */
                if(closed[counter] && !whole)
                        fprintf(fid,"%lf\t%lf\n",pprofx[counter][0],pprofy[counter][0]);

                /*Write blank line: */
                if(counter < nprof-1) fprintf(fid,"\n");
        }

        /*close domain outline file: */
        fclose(fid);

        cleanupandreturn: 
        return noerr;
}