Changeset 20163
- Timestamp:
- 02/15/16 16:20:30 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/test/NightlyRun/test2002.py
r20162 r20163 26 26 love = dlmread('../Data/love_numbers_10k.txt') 27 27 nlov=101 28 md.slr.love_h = love(1:nlov,2) # radial displacement (height) 29 md.slr.love_ k = love(1:nlov,3) # gravitational potential (phi)30 28 print love 29 md.slr.love_h = love[1][1:nlov] # radial displacement (height) 30 md.slr.love_k = love[2][1:nlov] # gravitational potential (phi) 31 31 #}}} 32 32 #mask: {{{ … … 36 36 icemask=ones(md.mesh.numberofvertices,1) 37 37 pos=find(mask==0) 38 icemask (pos)=-138 icemask[pos]=-1 39 39 pos=find(sum(mask(md.mesh.elements),2)<3) 40 icemask (md.mesh.elements(pos))=-140 icemask[md.mesh.elements[pos]]=-1 41 41 42 42 md.mask.ice_levelset=icemask 43 43 md.mask.ocean_levelset=zeros(md.mesh.numberofvertices,1) 44 44 pos=numpy.nonzero(md.mask.ice_levelset==1) 45 md.mask.ocean_levelset (pos)=145 md.mask.ocean_levelset[pos]=1 46 46 47 47 #make sure that the ice level set is all inclusive:
Note:
See TracChangeset
for help on using the changeset viewer.