Index: sm/trunk-jpl/src/m/mesh/planet/planetmesher.m
===================================================================
--- /issm/trunk-jpl/src/m/mesh/planet/planetmesher.m	(revision 17475)
+++ 	(revision )
@@ -1,35 +1,0 @@
-function md=planetmesher(md,varargin)
-%PLANETMESHER - create planet mesh using several packages
-%
-%   This routine creates a planet mesh using several packages: a custom made 
-%   ISSM mesher (called planetmixedmesh.m, method 'mixed') and the planettrimesh
-%   (method 'tria')
-%   where md is a @planet object, varargin is a list of options
-%
-% Usage: md=planetmesher(md,'method','mixed','radius',6378000,'angleresol',1);
-%        md=planetmesher(md,'method','tria','shape','iso','radius',6378000,'refinement',5);
-%
-%
-
-	%recover options
-	options=pairoptions(varargin{:});
-
-	method=getfieldvalue(options,'method','mixed');
-
-	if strcmpi(method,'mixed'),
-		%recover radius and angleresol: 
-		radius=getfieldvalue(options,'radius',6378000); %earth radius as default
-		angleresol=getfieldvalue(options,'angleresol',10); %10 degree resolution
-
-		%call mixed mesh 
-		md=planetmixedmesh(md,radius,angleresol);
-	else
-		%recover radius, shape and level of refinmenet
-		radius=getfieldvalue(options,'radius',6378000); %earth radius as default
-		refinement=getfieldvalue(options,'refinement',5); %refine 5 times
-		shape=getfieldvalue(options,'shape','ico'); 
-
-		%call triangular mesh
-		md=planettrimesh(md,shape,radius,refinement);
-	end
-end
Index: sm/trunk-jpl/src/m/mesh/planet/planetmixedmesh.m
===================================================================
--- /issm/trunk-jpl/src/m/mesh/planet/planetmixedmesh.m	(revision 17475)
+++ 	(revision )
@@ -1,97 +1,0 @@
-function md=planetmixedmesh(md,radius,angleresol)
-%PLANETMIXEDMESH: build 2d shell mesh
-%
-% Usage: md=planetmixedshell(md,radius,angleresol)
-%
-
-conv=pi/180;
-
-r=radius;
-theta=(0:angleresol:360)';
-theta=theta*conv;
-phi=(0:angleresol:180)';
-phi=phi(2:end-1)*conv;
-
-nphi=length(phi);
-nthetha=length(theta);
-
-x=zeros(nphi*nthetha,1);
-y=zeros(nphi*nthetha,1);
-z=zeros(nphi*nthetha,1);
-
-for i=1:nphi,
-	phii=phi(i);
-	x((i-1)*nthetha+1:i*nthetha)=r.*cos(theta).*sin(phii);
-	y((i-1)*nthetha+1:i*nthetha)=r.*sin(theta).*sin(phii);
-	z((i-1)*nthetha+1:i*nthetha)=r.*cos(phii);
-end
-
-quads=zeros((nthetha-1)*(nphi-1),4);
-
-%build elements
-for i=1:nphi-1,
-	for j=1:nthetha-1,
-		count=(i-1)*(nthetha-1)+j;
-
-		A=(i-1)*nthetha+j;
-		B=(i-1)*nthetha+j+1;
-		C=(i)*nthetha+j+1;
-		D=(i)*nthetha+j;
-
-		quads(count,:)=[A B C D];
-	end
-end
-
-%now, add polar elements:
-%first north: phi = 0;
-x=[r.*cos(0).*sin(0);x];
-y=[r.*sin(0).*sin(0);y];
-z=[r.*cos(0);z];
-
-%add 1 to quads: 
-quads=quads+1;
-
-%add tria elements:
-trias=zeros(nthetha-1,4);
-
-for i=1:nthetha-1,
-	A=1;
-	B=i+1;
-	C=i+2;
-	trias(i,:)=[A B C NaN];
-end
-
-quads=[trias;quads];
-
-% now add south pole: 
-phii=180*conv;
-x=[x;r.*cos(phii).*sin(phii)];
-y=[y;r.*sin(phii).*sin(phii)];
-z=[z;r.*cos(phii)];
-nods=length(x);
-
-%add tria elements:
-trias=zeros(nthetha-1,4);
-
-start=nods-nthetha;
-for i=1:nthetha-1,
-	A=start+i-1;
-	B=start+i;
-	C=nods;
-	trias(i,:)=[A B C NaN];
-end
-
-quads=[quads;trias];
-
-md.mesh.elements=quads;
-md.mesh.x=x;
-md.mesh.y=y;
-md.mesh.z=z;
-md.mesh.r=sqrt(x.^2+y.^2+z.^2);
-md.mesh.theta=acos(z./md.mesh.r);
-md.mesh.phi=atan2(y,x);
-
-md.mesh.numberofvertices=length(md.mesh.x);
-md.mesh.numberofelements=size(md.mesh.elements,1);
-
-md.mesh.dimension=3;
Index: /issm/trunk-jpl/src/m/mesh/planet/planettrimesh.m
===================================================================
--- /issm/trunk-jpl/src/m/mesh/planet/planettrimesh.m	(revision 17475)
+++ /issm/trunk-jpl/src/m/mesh/planet/planettrimesh.m	(revision 17476)
@@ -6,4 +6,5 @@
 
 results = sphere_tri(shape,refinement,radius);
+md.mesh=meshplanet(); %???
 md.mesh.x=results.vertices(:,1);
 md.mesh.y=results.vertices(:,2);
@@ -17,4 +18,2 @@
 md.mesh.numberofvertices=length(md.mesh.x);
 md.mesh.numberofelements=size(md.mesh.elements,1);
-
-md.mesh.dimension=3;