Index: /issm/trunk-jpl/src/m/coordsystems/gmtmask.m
===================================================================
--- /issm/trunk-jpl/src/m/coordsystems/gmtmask.m	(revision 20282)
+++ /issm/trunk-jpl/src/m/coordsystems/gmtmask.m	(revision 20283)
@@ -31,5 +31,5 @@
 	%First, write our lat,long file for gmt:
 	nv=length(lat);
-	dlmwrite('./all_vertices.txt',[long lat (1:nv)'],'delimiter','\t');
+	dlmwrite('./all_vertices.txt',[long lat (1:nv)'],'delimiter','\t','precision',10);
 
 	%Avoid bypassing of the ld library path by Matlab (:()
Index: /issm/trunk-jpl/src/m/coordsystems/gmtmask.py
===================================================================
--- /issm/trunk-jpl/src/m/coordsystems/gmtmask.py	(revision 20282)
+++ /issm/trunk-jpl/src/m/coordsystems/gmtmask.py	(revision 20283)
@@ -39,5 +39,5 @@
 	#First, write our lat,long file for gmt:
 	nv=lenlat
-	savetxt('./all_vertices.txt',[long, lat, arange(1,nv+1)],delimiter='\t')
+	savetxt('./all_vertices.txt',transpose([long, lat, arange(1,nv+1)]),delimiter='\t',fmt='%.10f')
 
 	#Avoid bypassing of the ld library path by Matlab (:()
@@ -68,5 +68,5 @@
 	oce_vertices=[]
 	while line:
-		ind=str2num(line.split()[3])
+		ind=int(float(line.split()[2]))-1;
 		oce_vertices.append(ind)
 		line=fid.readline()