Changeset 27752
- Timestamp:
- 05/15/23 15:56:29 (22 months ago)
- Location:
- issm/trunk-jpl/src/m/coordsystems
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/m/coordsystems/gmtmask.m
r27589 r27752 40 40 41 41 %figure out which vertices are on the ocean, which one on the continent: 42 gmt_select_options='-h0 -Df -R0/360/-90/90 -A0 -JQ180/200 -Nk/s/s/k/s'; 42 % 43 % NOTE: Remove -Ve option to enable warnings if this method is not working 44 % expected 45 % 46 gmt_select_options='-Ve -h0 -Df -R0/360/-90/90 -A0 -JQ180/200 -Nk/s/s/k/s'; 43 47 [status,result]=system(['gmt select ./' filename_all ' ' gmt_select_options ' > ./' filename_oce]); 44 48 if status~=0, -
issm/trunk-jpl/src/m/coordsystems/gmtmask.py
r27589 r27752 44 44 45 45 # Figure out which vertices are on the ocean, which one on the continent: 46 gmt_select_options = '-h0 -Df -R0/360/-90/90 -A0 -JQ180/200 -Nk/s/s/k/s' 46 # 47 # NOTE: Remove -Ve option to enable warnings if this method is not working 48 # expected 49 # 50 gmt_select_options = '-Ve -h0 -Df -R0/360/-90/90 -A0 -JQ180/200 -Nk/s/s/k/s' 47 51 subproc_cmd = 'gmt select ./all_vertices.txt ' + gmt_select_options + ' > ./oce_vertices.txt' 48 52 subproc = subprocess.Popen(subproc_cmd, shell=True, stdout=subprocess.PIPE, stderr=subprocess.PIPE, universal_newlines=True)
Note:
See TracChangeset
for help on using the changeset viewer.