Changeset 27201

Timestamp:

08/10/22 13:41:19 (3 years ago)

Author:

jdquinn

Message:

BUG: gmt executable syntax

Location:

issm/trunk-jpl/src/m/coordsystems

Files:

• gmtmask.m (modified) (1 diff)
• gmtmask.py (modified) (1 diff)
• gmtmaskparallel.m (modified) (1 diff)

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

@@ -36 +36 @@
 
         %figure out which vertices are on the ocean, which one on the continent:
-        [status,result] = system(['gmtselect ./' filename_all ' -h0 -Df -R0/360/-90/90  -A0 -JQ180/200 -Nk/s/s/k/s > ./' filename_oce]);
+        [status,result] = system(['gmt select ./' filename_all ' -h0 -Df -R0/360/-90/90  -A0 -JQ180/200 -Nk/s/s/k/s > ./' filename_oce]);
         if status~=0,
                 error(result);

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

@@ -45 +45 @@
 
     #figure out which vertices are on the ocean, which one on the continent:
-    subprocess.call('gmtselect ./ all_vertices.txt -h0 -Df -R0/360/-90/90 -A0 -JQ180/200 -Nk/s/s/k/s > ./oce_vertices.txt', shell=True)
+    subprocess.call('gmt select ./ all_vertices.txt -h0 -Df -R0/360/-90/90 -A0 -JQ180/200 -Nk/s/s/k/s > ./oce_vertices.txt', shell=True)
 
     #read the con_vertices.txt file and flag our mesh vertices on the continent

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

@@ -42 +42 @@
         fid=fopen('xjobs.script','w');
         for i=1:length(nnv)-1,
-                fprintf(fid,'%s gmtselect ./all_vertices%i.txt -h0 -Df -R0/360/-90/90  -A0 -JQ180/200 -Nk/s/s/k/s > ./oce_vertices%i.txt\n',gmtpath,i,i);
+                fprintf(fid,'%s gmt select ./all_vertices%i.txt -h0 -Df -R0/360/-90/90  -A0 -JQ180/200 -Nk/s/s/k/s > ./oce_vertices%i.txt\n',gmtpath,i,i);
         end
         fclose(fid);