Changeset 19981

Timestamp:

01/22/16 09:47:38 (9 years ago)

Author:

Eric.Larour

Message:

CHG: address issue of gmtmask not accepting large input files (i.e. larger than ~50-100,000 vertices). Implement
recursive call.

File:

• issm/trunk-jpl/src/m/coordsystems/gmtmask.m (modified) (1 diff)

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

@@ -6 +6 @@
 %
 
-
+        %Check lat and long size is not more than 50,000; If so, recursively call gmtmask: 
+        if length(lat)>50000,
+                for i=1:50000:length(lat),
+                        j=i+50000-1;
+                        if j>length(lat),
+                                j=length(lat);
+                        end
+                        mask(i:j)=gmtmask(lat(i:j),long(i:j));
+                end
+                return
+        end
+        
         %First, write our lat,long file for gmt:
         nv=length(lat);