Index: /issm/trunk-jpl/src/c/shared/Exp/DomainOutlineRead.cpp.bak
===================================================================
--- /issm/trunk-jpl/src/c/shared/Exp/DomainOutlineRead.cpp.bak	(revision 11886)
+++ /issm/trunk-jpl/src/c/shared/Exp/DomainOutlineRead.cpp.bak	(revision 11886)
@@ -0,0 +1,162 @@
+/*!\file:  DomainOutlineRead.cpp
+ * \brief DomainOutlineRead.c: read 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 "../../objects/objects.h"
+#include "../Exceptions/exceptions.h"
+#include "../../Container/DataSet.h"
+
+int DomainOutlineRead(int* pnprof,int** pprofnvertices,double*** ppprofx,double*** ppprofy,bool** pclosed,char* domainname,bool whole=true){
+
+	
+	/*indexing: */
+	int i,counter;
+
+	/*I/O: */
+	FILE* fid=NULL;
+	char chardummy[256];
+	double  ddummy;
+
+	/*output: */
+	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
+
+	/*For each profile: */
+	int n;
+	double* x=NULL;
+	double* y=NULL;
+	bool cl;
+
+	/*open domain outline file for reading: */
+	if ((fid=fopen(domainname,"r"))==NULL){
+		_error_("%s%s","could not find domain file ",domainname); 
+	}
+
+	/*Do a first pass through the domainname file, to figure out how many profiles 
+	 *we need to read: */
+	nprof=1;
+	for(;;){
+		fscanf(fid,"%256s %256s\n",chardummy,chardummy);
+		fscanf(fid,"%256s %256s\n",chardummy,chardummy);
+		fscanf(fid,"%256s %256s %256s %256s\n",chardummy,chardummy,chardummy,chardummy);
+		fscanf(fid,"%20u %256s\n",&n,chardummy);
+		fscanf(fid,"%256s %256s %256s %256s %256s\n",chardummy,chardummy,chardummy,chardummy,chardummy);
+		for (i=0;i<n;i++){
+			fscanf(fid,"%20lf %20lf\n",&ddummy,&ddummy);
+		}
+		/*Ok, we have faked one profile reading, check whether we are at the end of the file, otherwise, keep fake reading next profile:*/
+		if (feof(fid)){
+			break;
+		}
+		nprof++;
+	}
+	
+	/*Allocate and initialize all the profiles: */
+	profnvertices=(int*)xmalloc(nprof*sizeof(int));
+	pprofx=(double**)xmalloc(nprof*sizeof(double*));
+	pprofy=(double**)xmalloc(nprof*sizeof(double*));
+	for (i=0;i<nprof;i++){
+		pprofx[i]=NULL;
+		pprofy[i]=NULL;
+	}
+	closed=(bool*)xmalloc(nprof*sizeof(bool));
+
+	/*Reaset file pointer to beginning of file: */
+	fseek(fid,0,SEEK_SET);
+
+	/*Start reading profiles: */
+	for(counter=0;counter<nprof;counter++){
+
+		/*Skip header: */
+		fscanf(fid,"%256s %256s\n",chardummy,chardummy);
+		fscanf(fid,"%256s %256s\n",chardummy,chardummy);
+		fscanf(fid,"%256s %256s %256s %256s\n",chardummy,chardummy,chardummy,chardummy);
+		
+		/*Get number of profile vertices: */
+		fscanf(fid,"%20u %256s\n",&n,chardummy);
+	
+		/*Skip next line: */
+		fscanf(fid,"%256s %256s %256s %256s %256s\n",chardummy,chardummy,chardummy,chardummy,chardummy);
+
+		/*Allocate vertices: */
+		printf("number of n: %i\n",n);
+		x=(double*)xmalloc(n*sizeof(double));
+		y=(double*)xmalloc(n*sizeof(double));
+		
+
+		/*Read vertices: */
+		for (i=0;i<n;i++){
+			fscanf(fid,"%20lf %20lf\n",&x[i],&y[i]);
+		}
+
+		/*Now check that we are dealing with open contours: */
+		cl=false;
+		if((x[0]==x[n-1]) && (y[0]==y[n-1])){
+			cl=true;
+			if (!whole) {
+				n=n-1;
+			}
+		}
+
+		/*Assign pointers: */
+		profnvertices[counter]=n;
+		pprofx[counter]=x;
+		pprofy[counter]=y;
+		closed[counter]=cl;
+	}
+
+	/*close domain outline file: */
+	fclose(fid);
+
+	/*Assign output pointers: */
+	*pnprof=nprof;
+	*pprofnvertices=profnvertices;
+	*ppprofx=pprofx;
+	*ppprofy=pprofy;
+	if(pclosed)*pclosed=closed;
+	else       xfree((void**)&closed);
+}
+
+DataSet* DomainOutlineRead(char* domainname,bool whole=true){
+
+	/*indexing: */
+	int i;
+
+	/*intermediary: */
+	int nprof;
+	int* profnvertices=NULL;
+	double** pprofx=NULL;
+	double** pprofy=NULL;
+
+	Contour* contour=NULL;
+
+	/*output: */
+	DataSet* domain=NULL;
+
+	/*get domain outline from intermediary function:*/
+	DomainOutlineRead(&nprof,&profnvertices,&pprofx, &pprofy, NULL,domainname,whole);
+
+	double* a=pprofx[0]; for(i=0;i<4;i++)printf("%g\n",a[i]);
+	a=pprofy[0]; for(i=0;i<4;i++)printf("%g\n",a[i]);
+	printf("%i %i\n",profnvertices[0],nprof);
+
+	/*now create dataset of contours: */
+	domain=new DataSet(0);
+	printf("ok1\n");
+
+	for(i=0;i<nprof;i++){
+		domain->AddObject(new Contour(i,profnvertices[i],pprofx[i],pprofy[i],1));
+	}
+	printf("ok2\n");
+
+	return domain;
+}