Changeset 27589 for issm/trunk-jpl/src/m/coordsystems/gmtmask.m

Timestamp:

02/17/23 13:41:46 (2 years ago)

Author:

jdquinn

Message:

CHG: Reverting to GMT 6.0.0 because executables under 6.4.0 do not link statically to libgmt; find gmtselect if it is available; removed test with restart from packaging scripts; cleanup

File:

• issm/trunk-jpl/src/m/coordsystems/gmtmask.m (modified) (2 diffs)

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

@@ -5 +5 @@
 %      mask.ocean = gmtmask(md.mesh.lat,md.mesh.long);
 %
+%       TODO: Standardize discovery of GMT bin path and whether or not we have GMT 
+%       modules (i.e. `gmt select`) between this file, gmtmask.py, and 
+%       gmtmaskparallel.m
 
         %are we doing a recursive call? 
@@ -31 +34 @@
         %First, write our lat,long file for gmt:
         nv=length(lat);
-        filename_all = ['all_vertices-' num2str(feature('GetPid')) '.txt']; 
-        filename_oce = ['oce_vertices-' num2str(feature('GetPid')) '.txt']; 
+        filename_suffix=[num2str(feature('GetPid')) '.txt'];
+        filename_all=['all_vertices-' filename_suffix]; 
+        filename_oce=['oce_vertices-' filename_suffix];
         dlmwrite(filename_all,[long lat (1:nv)'],'delimiter','\t','precision',10);
 
         %figure out which vertices are on the ocean, which one on the continent:
-        [status,result] = system(['gmt select ./' filename_all ' -h0 -Df -R0/360/-90/90 -A0 -JQ180/200 -Nk/s/s/k/s > ./' filename_oce]);
+        gmt_select_options='-h0 -Df -R0/360/-90/90 -A0 -JQ180/200 -Nk/s/s/k/s';
+        [status,result]=system(['gmt select ./' filename_all ' ' gmt_select_options ' > ./' filename_oce]);
         if status~=0,
-                error(result);
+                [status,result] = system(['gmtselect ./' filename_all ' ' gmt_select_options ' > ./' filename_oce]);
+                if status~=0,
+                        error(result);
+                end
         end