Changeset 15943
- Timestamp:
- 08/26/13 10:33:35 (12 years ago)
- Location:
- issm/trunk-jpl/src/m
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/m/classes/mask.py
r15942 r15943 44 44 md = checkfield(md,'mask.vertexonfloatingice' ,'size',[md.mesh.numberofvertices],'values',[0,1]) 45 45 md = checkfield(md,'mask.vertexongroundedice' ,'size',[md.mesh.numberofvertices],'values',[0,1]) 46 md = checkfield(md,'mask.icelevelset' ,'size',[md.mesh.numberofvertices]) 47 isice=numpy.array(md.mask.icelevelset>0,int) 46 md = checkfield(md,'mask.groundedice_levelset','size',[md.mesh.numberofvertices]) 47 md = checkfield(md,'mask.ice_levelset' ,'size',[md.mesh.numberofvertices]) 48 isice=numpy.array(md.mask.ice_levelset>0,int) 48 49 totallyicefree=(numpy.sum(isice[md.mesh.elements-1],axis=1)==0).astype(int) 49 50 if any(totallyicefree): 50 raise TypeError("elements with no ice not implemented yet, each element should have at least one vertex with md.mask.ice levelset > 0")51 raise TypeError("elements with no ice not implemented yet, each element should have at least one vertex with md.mask.ice_levelset > 0") 51 52 52 53 return md -
issm/trunk-jpl/src/m/classes/model/model.m
r15775 r15943 216 216 md.mask.vertexonrock=project2d(md,md.mask.vertexonrock,1); 217 217 end 218 md.mask.icelevelset=project2d(md,md.mask.icelevelset,1); 218 md.mask.groundedice_levelset=project2d(md,md.mask.groundedice_levelset,1); 219 md.mask.ice_levelset=project2d(md,md.mask.ice_levelset,1); 219 220 220 221 %lat long … … 742 743 md.mask.vertexongroundedice=project3d(md,'vector',md.mask.vertexongroundedice,'type','node'); 743 744 md.mask.vertexonrock=project3d(md,'vector',md.mask.vertexonrock,'type','node'); 744 md.mask.icelevelset=project3d(md,'vector',md.mask.icelevelset,'type','node'); 745 md.mask.groundedice_levelset=project3d(md,'vector',md.mask.ice_levelset,'type','node'); 746 md.mask.ice_levelset=project3d(md,'vector',md.mask.ice_levelset,'type','node'); 745 747 if ~isnan(md.inversion.cost_functions_coefficients),md.inversion.cost_functions_coefficients=project3d(md,'vector',md.inversion.cost_functions_coefficients,'type','node');end; 746 748 if ~isnan(md.inversion.min_parameters),md.inversion.min_parameters=project3d(md,'vector',md.inversion.min_parameters,'type','node');end; -
issm/trunk-jpl/src/m/parameterization/setmask.m
r15749 r15943 46 46 md.mask.elementongroundedice=elementongroundedice; 47 47 md.mask.vertexongroundedice=vertexongroundedice; 48 md.mask.ice levelset=ones(md.mesh.numberofvertices,1);48 md.mask.ice_levelset=ones(md.mesh.numberofvertices,1); -
issm/trunk-jpl/src/m/parameterization/setmask.py
r15749 r15943 46 46 md.mask.elementongroundedice = elementongroundedice 47 47 md.mask.vertexongroundedice = vertexongroundedice 48 md.mask.ice levelset = numpy.ones(md.mesh.numberofvertices,bool)48 md.mask.ice_levelset = numpy.ones(md.mesh.numberofvertices,bool) 49 49 50 50 return md
Note:
See TracChangeset
for help on using the changeset viewer.