Index: /issm/trunk-jpl/src/m/coordsystems/gmtmask.m
===================================================================
--- /issm/trunk-jpl/src/m/coordsystems/gmtmask.m	(revision 27751)
+++ /issm/trunk-jpl/src/m/coordsystems/gmtmask.m	(revision 27752)
@@ -40,5 +40,9 @@
 
 	%figure out which vertices are on the ocean, which one on the continent:
-	gmt_select_options='-h0 -Df -R0/360/-90/90 -A0 -JQ180/200 -Nk/s/s/k/s';
+	%
+	% NOTE: Remove -Ve option to enable warnings if this method is not working 
+	%		expected
+	% 
+	gmt_select_options='-Ve -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,
Index: /issm/trunk-jpl/src/m/coordsystems/gmtmask.py
===================================================================
--- /issm/trunk-jpl/src/m/coordsystems/gmtmask.py	(revision 27751)
+++ /issm/trunk-jpl/src/m/coordsystems/gmtmask.py	(revision 27752)
@@ -44,5 +44,9 @@
 
     # Figure out which vertices are on the ocean, which one on the continent:
-    gmt_select_options = '-h0 -Df -R0/360/-90/90 -A0 -JQ180/200 -Nk/s/s/k/s'
+    #
+    # NOTE: Remove -Ve option to enable warnings if this method is not working 
+    #       expected
+    #
+    gmt_select_options = '-Ve -h0 -Df -R0/360/-90/90 -A0 -JQ180/200 -Nk/s/s/k/s'
     subproc_cmd = 'gmt select ./all_vertices.txt ' + gmt_select_options + ' > ./oce_vertices.txt'
     subproc = subprocess.Popen(subproc_cmd, shell=True, stdout=subprocess.PIPE, stderr=subprocess.PIPE, universal_newlines=True)