source: issm/oecreview/Archive/19101-20495/ISSM-20059-20060.diff@ 20498

../trunk-jpl/src/m/coordsystems/gmtmaskparallel.m

@@ +1 @@
+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');