Changeset 9654
- Timestamp:
- 09/07/11 13:19:25 (14 years ago)
- Location:
- issm/trunk/src/m/model
- Files:
-
- 5 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/src/m/model/bamg.m
r9619 r9654 330 330 md.z=zeros(md.numberofnodes,1); 331 331 md.nodeonbed=ones(md.numberofnodes,1); 332 md. nodeonwater=zeros(md.numberofnodes,1);332 md.mask.vertexonwater=zeros(md.numberofnodes,1); 333 333 md.nodeonsurface=ones(md.numberofnodes,1); 334 334 md.elementonbed=ones(md.numberofelements,1); -
issm/trunk/src/m/model/collapse.m
r9646 r9654 50 50 md.elements_type=project2d(md,md.elements_type,1); 51 51 end 52 md.nodeonhutter=project2d(md,md.nodeonhutter,1);53 md.nodeonmacayeal=project2d(md,md.nodeonmacayeal,1);54 md.nodeonpattyn=project2d(md,md.nodeonpattyn,1);55 md.nodeonstokes=project2d(md,md.nodeonstokes,1);56 52 57 53 %boundary conditions -
issm/trunk/src/m/model/extrude.m
r9646 r9654 158 158 md.elements_type=zeros(number_el3d,1); 159 159 md.elements_type=project3d(md,'vector',oldelements_type,'type','element'); 160 md.nodeonhutter=project3d(md,'vector',md.nodeonhutter,'type','node');161 md.nodeonmacayeal=project3d(md,'vector',md.nodeonmacayeal,'type','node');162 md.nodeonpattyn=project3d(md,'vector',md.nodeonpattyn,'type','node');163 md.nodeonstokes=project3d(md,'vector',md.nodeonstokes,'type','node');164 160 end 165 161 -
issm/trunk/src/m/model/ismodelselfconsistent.m
r9650 r9654 386 386 end 387 387 %CHECK THAT ROTATION IS IN THE (X,Y) PLANE FOR 2D MODELS 388 if any(md.nodeonmacayeal),389 pos=find(sum(isnan(md.diagnostic_ref),2)==0 & md. nodeonmacayeal);388 if md.dim==2, 389 pos=find(sum(isnan(md.diagnostic_ref),2)==0 & md.vertices_type==2); 390 390 if any(md.diagnostic_ref(pos,3:5)~=0); 391 391 message(['model not consistent: model ' md.miscellaneous.name ' has problem with rotated spc. The rotation should be in the (x,y) plane for 2D diagnostic models (nodeonmacayeal)']); -
issm/trunk/src/m/model/setelementstype.m
r9641 r9654 73 73 nodeonhutter=zeros(md.numberofnodes,1); 74 74 nodeonhutter(md.elements(find(hutterflag),:))=1; 75 md.nodeonhutter=nodeonhutter;76 75 md.elements_type(find(hutterflag))=1; 77 76
Note:
See TracChangeset
for help on using the changeset viewer.