Index: /issm/trunk-jpl/src/m/coordsystems/gmtmaskparallel.m
===================================================================
--- /issm/trunk-jpl/src/m/coordsystems/gmtmaskparallel.m	(revision 20060)
+++ /issm/trunk-jpl/src/m/coordsystems/gmtmaskparallel.m	(revision 20060)
@@ -0,0 +1,49 @@
+function mask = gmtmaskparallel(lat,long,ncores)
+%GMTMASKPARALLEL- parallel driver for the gmtmask utility
+%
+%   Usage:
+%      mask.ocean = gmtmaskparallel(md.mesh.lat,md.mesh.long,8);
+%
+
+	%First, write our lat,long file for gmt:
+	nv=length(lat);
+
+	%Split: 
+	nnv=1:floor(nv/ncores):nv;
+	nnv(end)=nv+1;
+
+	%For each segment, write all vertices file: 
+	for i=1:length(nnv)-1,
+		dlmwrite(['./all_vertices' num2str(i) '.txt'],[long(nnv(i):nnv(i+1)-1) lat(nnv(i):nnv(i+1)-1) (nnv(i):nnv(i+1)-1)'],'delimiter','\t');
+	end
+
+	if ismac,
+		dyld_library_path_old=getenv('DYLD_LIBRARY_PATH');
+		setenv('DYLD_LIBRARY_PATH',[ issmdir '/externalpackages/curl/install/lib:' issmdir '/externalpackages/hdf5/install/lib:' issmdir '/externalpackages/netcdf/install/lib' ]);
+	end
+
+	%Build xjobs script:
+	fid=fopen('xjobs.script','w');
+	for i=1:length(nnv)-1,
+		fprintf(fid,'gmt gmtselect ./all_vertices%i.txt -h0 -Df -R0/360/-90/90  -A0 -JQ180/200 -Nk/s/s/k/s > ./oce_vertices%i.txt\n',i,i);
+	end
+	fclose(fid);
+
+	%Call xjobs: 
+	system(sprintf('xjobs -j %i -s ./xjobs.script',ncores));
+
+	%reset DYLD_LIBRARY_PATH to what it was: 
+	if ismac,
+		setenv('DYLD_LIBRARY_PATH',dyld_library_path_old);
+	end
+
+	%concatenate: 
+	system('cat oce_vertices*.txt | grep -v Command | awk ''{printf("%s\n",$3);}''> vertices.txt');
+
+	%read the vertices.txt file and flag our mesh vertices on the continent
+	flags=dlmread('./vertices.txt');
+
+	mask=zeros(nv,1);
+	mask(flags)=1;
+	
+	system('rm -rf ./all_vertices*.txt ./oce_vertices*.txt vertices.txt ./gmt.history');