- Timestamp:
- 06/07/17 10:50:54 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/branches/trunk-larour-NatGeoScience2016/src/m/coordsystems/gmtmask.py
r21069 r21759 1 1 from MatlabFuncs import * 2 2 from model import * 3 from numpy import * 3 import numpy as np 4 4 from os import getenv, putenv 5 5 import subprocess … … 12 12 # 13 13 lenlat=len(lat) 14 mask= empty(lenlat)14 mask=np.empty(lenlat) 15 15 16 16 #are we doing a recursive call? … … 38 38 #First, write our lat,long file for gmt: 39 39 nv=lenlat 40 savetxt('./all_vertices.txt',transpose([long, lat,arange(1,nv+1)]),delimiter='\t',fmt='%.10f')40 np.savetxt('./all_vertices.txt',np.transpose([long, lat, np.arange(1,nv+1)]),delimiter='\t',fmt='%.10f') 41 41 42 42 #Avoid bypassing of the ld library path by Matlab (:() … … 72 72 fid.close() 73 73 74 mask= zeros([nv,1])74 mask=np.zeros([nv,1]) 75 75 mask[oce_vertices]=1 76 76
Note:
See TracChangeset
for help on using the changeset viewer.