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

Index: ../trunk-jpl/src/m/coordsystems/gmtmask.py
===================================================================
--- ../trunk-jpl/src/m/coordsystems/gmtmask.py  (revision 20282)
+++ ../trunk-jpl/src/m/coordsystems/gmtmask.py  (revision 20283)
@@ -38,7 +38,7 @@
        
        #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,7 +67,7 @@
        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()
Index: ../trunk-jpl/src/m/coordsystems/gmtmask.m
===================================================================
--- ../trunk-jpl/src/m/coordsystems/gmtmask.m   (revision 20282)
+++ ../trunk-jpl/src/m/coordsystems/gmtmask.m   (revision 20283)
@@ -30,7 +30,7 @@
        
        %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,