Changeset 20283

Timestamp:

03/01/16 20:18:15 (9 years ago)

Author:

Eric.Larour

Message:

CHG: gmtmask.py: issue with indexing! (matlab bias)
gmtmask.py and .m: increase precision to 10 digits for lat,long! otherwise, gmtmask
is returning biased locations!

Location:

issm/trunk-jpl/src/m/coordsystems

Files:

• gmtmask.m (modified) (1 diff)
• gmtmask.py (modified) (2 diffs)

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

@@ -31 +31 @@
         %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 (:()

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

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