source: issm/oecreview/Archive/19101-20495/ISSM-20282-20283.diff@ 20498

../trunk-jpl/src/m/coordsystems/gmtmask.py

@@ -38 +38 @@
         
         #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 (:()
         try:
@@ -67 +67 @@
         line=fid.readline()
         oce_vertices=[]
         while line:
-                ind=str2num(line.split()[3])
+                ind=int(float(line.split()[2]))-1;
                 oce_vertices.append(ind)
                 line=fid.readline()
         fid.close()

../trunk-jpl/src/m/coordsystems/gmtmask.m

@@ -30 +30 @@
         
         %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 (:()
         if ismac,