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,