Changeset 20163


Ignore:
Timestamp:
02/15/16 16:20:30 (9 years ago)
Author:
dlcheng
Message:

CHG: Resolving python syntax errors for test2002

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk-jpl/test/NightlyRun/test2002.py

    r20162 r20163  
    2626love = dlmread('../Data/love_numbers_10k.txt')
    2727nlov=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 
     28print love
     29md.slr.love_h = love[1][1:nlov]  # radial displacement (height)
     30md.slr.love_k = love[2][1:nlov]  # gravitational potential (phi)
    3131#}}}
    3232#mask:  {{{
     
    3636icemask=ones(md.mesh.numberofvertices,1)
    3737pos=find(mask==0)
    38 icemask(pos)=-1
     38icemask[pos]=-1
    3939pos=find(sum(mask(md.mesh.elements),2)<3)
    40 icemask(md.mesh.elements(pos))=-1
     40icemask[md.mesh.elements[pos]]=-1
    4141
    4242md.mask.ice_levelset=icemask
    4343md.mask.ocean_levelset=zeros(md.mesh.numberofvertices,1)
    4444pos=numpy.nonzero(md.mask.ice_levelset==1)
    45 md.mask.ocean_levelset(pos)=1
     45md.mask.ocean_levelset[pos]=1
    4646
    4747#make sure that the ice level set is all inclusive:
Note: See TracChangeset for help on using the changeset viewer.