Index: /issm/trunk/src/m/classes/public/ThicknessCorrection.m
===================================================================
--- /issm/trunk/src/m/classes/public/ThicknessCorrection.m	(revision 40)
+++ /issm/trunk/src/m/classes/public/ThicknessCorrection.m	(revision 41)
@@ -19,5 +19,5 @@
 end
 
-in=ArgusContourToMesh(md.elements,md.x,md.y,expread(filename,1),'node',1);
+in=ContourToMesh(md.elements,md.x,md.y,expread(filename,1),'node',1);
 pos_shelf=find(in & ~md.gridonicesheet);
 pos_sheet=find(in & md.gridonicesheet);
Index: /issm/trunk/src/m/classes/public/domaintomesh.m
===================================================================
--- /issm/trunk/src/m/classes/public/domaintomesh.m	(revision 40)
+++ /issm/trunk/src/m/classes/public/domaintomesh.m	(revision 41)
@@ -29,5 +29,5 @@
 end
 
-[rinterp]=ArgusContourToMesh(md.elements,md.x,md.y,expread(domainname,1),interptype,1);
+[rinterp]=ContourToMesh(md.elements,md.x,md.y,expread(domainname,1),interptype,1);
 interp=IMdb('select matrix from rinterp');
 IMdb('drop rinterp');
Index: /issm/trunk/src/m/classes/public/geography.m
===================================================================
--- /issm/trunk/src/m/classes/public/geography.m	(revision 40)
+++ /issm/trunk/src/m/classes/public/geography.m	(revision 41)
@@ -50,5 +50,5 @@
 	elementoniceshelf=ones(md.numberofelements,1);
 else
-	[gridoniceshelf,elementoniceshelf]=ArgusContourToMesh(elements,x,y,expread(iceshelfname,1),'element and node',2);
+	[gridoniceshelf,elementoniceshelf]=ContourToMesh(elements,x,y,expread(iceshelfname,1),'element and node',2);
 end
 
@@ -57,5 +57,5 @@
 	elementonicesheet=zeros(md.numberofelements,1);
 else
-	[gridonicesheet,elementonicesheet]=ArgusContourToMesh(elements,x,y,expread(icesheetname,1),'element and node',2);
+	[gridonicesheet,elementonicesheet]=ContourToMesh(elements,x,y,expread(icesheetname,1),'element and node',2);
 end
 
Index: /issm/trunk/src/m/classes/public/meshexprefine.m
===================================================================
--- /issm/trunk/src/m/classes/public/meshexprefine.m	(revision 40)
+++ /issm/trunk/src/m/classes/public/meshexprefine.m	(revision 41)
@@ -39,5 +39,5 @@
 
 %Read domainname file into a matlab array (x,y):
-refinearea=ArgusContourToMesh(md.elements,md.x,md.y,expread(domainname,1),'element',1);
+refinearea=ContourToMesh(md.elements,md.x,md.y,expread(domainname,1),'element',1);
 aires=area(md);
 
Index: /issm/trunk/src/m/classes/public/modelextract.m
===================================================================
--- /issm/trunk/src/m/classes/public/modelextract.m	(revision 40)
+++ /issm/trunk/src/m/classes/public/modelextract.m	(revision 41)
@@ -43,5 +43,5 @@
 		%ok, flag_elem elements
 		A=expread(area,1);
-		flag_elem=ArgusContourToMesh(md1.elements(:,1:3),md1.x,md1.y,A,'element',1);
+		flag_elem=ContourToMesh(md1.elements(:,1:3),md1.x,md1.y,A,'element',1);
 		%check that the elements do not cross the contour -> flag as 0
 		pos=find(flag_elem);
Index: /issm/trunk/src/m/classes/public/recover_areas.m
===================================================================
--- /issm/trunk/src/m/classes/public/recover_areas.m	(revision 40)
+++ /issm/trunk/src/m/classes/public/recover_areas.m	(revision 41)
@@ -60,5 +60,5 @@
 			end
 			%ok, flag elements
-			flag=ArgusContourToMesh(md.elements(:,1:3),md.x,md.y,expread(area,1),'element',1);
+			flag=ContourToMesh(md.elements(:,1:3),md.x,md.y,expread(area,1),'element',1);
 		end
 		if invert,