Changeset 15600
- Timestamp:
- 07/24/13 17:09:03 (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/m/classes/mask.py
r15596 r15600 49 49 md = checkfield(md,'mask.vertexonwater' ,'size',[md.mesh.numberofvertices],'values',[0,1]) 50 50 md = checkfield(md,'mask.icelevelset' ,'size',[md.mesh.numberofvertices]) 51 isice=numpy.zeros_like(md.mesh.numberofvertices,int)52 53 54 totalyicefree=(numpy.sum(isice[md.mesh.elements-1],axis=1)==0).astype(int)55 if any(totalyicefree)56 51 isice=numpy.zeros(md.mesh.numberofvertices,int) 52 pos=npy.nonzero(md.mask.icelevelset>0) 53 isice[pos]=1 54 totallyicefree=(numpy.sum(isice[md.mesh.elements-1],axis=1)==0).astype(int) 55 if any(totallyicefree): 56 raise TypeError("elements with no ice not implemented yet, each element should have at least one vertex with md.mask.icelevelset > 0") 57 57 58 58 return md
Note:
See TracChangeset
for help on using the changeset viewer.