Changeset 1298
- Timestamp:
- 07/13/09 08:56:58 (16 years ago)
- Location:
- issm/trunk/src/m
- Files:
-
- 1 deleted
- 5 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/src/m/classes/public/extrude.m
r1241 r1298 140 140 %dead grids 141 141 md.deadgrids=ones(md.numberofgrids,1); 142 md.deadgrids(md.elements(md.elements_type(:,1)~= macayealenum,:))=0;%non macayeal grids are not dead142 md.deadgrids(md.elements(md.elements_type(:,1)~=MacAyealEnum,:))=0;%non macayeal grids are not dead 143 143 md.deadgrids(find(md.gridonbed))=0;%grids from elements on bed are not dead 144 144 end -
issm/trunk/src/m/classes/public/ismodelselfconsistent.m
r1222 r1298 41 41 bool=0;return; 42 42 end 43 if any(ones(md.numberofelements,1)-((md.elements_type(:,1)== hutterenum) + (md.elements_type(:,1)==macayealenum) + (md.elements_type(:,1)==pattynenum)))43 if any(ones(md.numberofelements,1)-((md.elements_type(:,1)==HutterEnum) + (md.elements_type(:,1)==MacAyealEnum) + (md.elements_type(:,1)==PattynEnum))) 44 44 disp(['Types of elements have not been set properly, run setelementstype first']) 45 45 bool=0;return; 46 46 end 47 if any(ones(md.numberofelements,1)-((md.elements_type(:,2)== stokesenum) + (md.elements_type(:,2)==noneenum)))47 if any(ones(md.numberofelements,1)-((md.elements_type(:,2)==StokesEnum) + (md.elements_type(:,2)==NoneEnum))) 48 48 disp(['Types of elements have not been set properly, run setelementstype first']) 49 49 bool=0;return; 50 50 end 51 51 if strcmpi(md.type,'2d'), 52 if (ismember( pattynenum,md.elements_type(:,1)) | ismember(stokesenum,md.elements_type(:,2))),52 if (ismember(PattynEnum,md.elements_type(:,1)) | ismember(StokesEnum,md.elements_type(:,2))), 53 53 disp(['For a 2d model, only MacAyeal''s and Hutter''s elements are allowed']); 54 54 bool=0;return; … … 186 186 187 187 %HUTTER ON ICESHELF WARNING 188 if any(md.elements_type(:,1)== hutterenum & md.elementoniceshelf),188 if any(md.elements_type(:,1)==HutterEnum & md.elementoniceshelf), 189 189 disp(sprintf('\n !!! Warning: Hutter''s model is not consistent on ice shelves !!!\n')); 190 190 end … … 314 314 315 315 %HUTTER 316 if any(md.elements_type(:,1)== hutterenum);316 if any(md.elements_type(:,1)==HutterEnum); 317 317 disp(['The model has Hutter''s elements. Impossible to compute parameters']); 318 318 bool=0;return; -
issm/trunk/src/m/classes/public/mechanicalproperties.m
r1236 r1298 20 20 error('only 2d model supported yet'); 21 21 end 22 if any(md.elements_type(:,1)~= macayealenum),22 if any(md.elements_type(:,1)~=MacAyealEnum), 23 23 disp('Warning: the model has some non macayeal elements. These will be treated like MacAyeal''s elements'); 24 24 end -
issm/trunk/src/m/classes/public/plot/plot_elementstype.m
r27 r1298 25 25 if is2d 26 26 %Hutter elements 27 posH=find(data(:,1)== hutterenum);27 posH=find(data(:,1)==HutterEnum); 28 28 A=elements(posH,1); B=elements(posH,2); C=elements(posH,3); 29 p1=patch( 'Faces', [A B C], 'Vertices', [x y z],'CData', hutterenum,'FaceColor','flat','EdgeColor',edgecolor);29 p1=patch( 'Faces', [A B C], 'Vertices', [x y z],'CData',HutterEnum,'FaceColor','flat','EdgeColor',edgecolor); 30 30 %MacAyeal element 31 posH=find(data(:,1)== macayealenum);31 posH=find(data(:,1)==MacAyealEnum); 32 32 A=elements(posH,1); B=elements(posH,2); C=elements(posH,3); 33 p2=patch( 'Faces', [A B C], 'Vertices', [x y z],'CData', macayealenum,'FaceColor','flat','EdgeColor',edgecolor);33 p2=patch( 'Faces', [A B C], 'Vertices', [x y z],'CData',MacAyealEnum,'FaceColor','flat','EdgeColor',edgecolor); 34 34 %Pattyn element 35 posH=find(data(:,1)== pattynenum);35 posH=find(data(:,1)==PattynEnum); 36 36 A=elements(posH,1); B=elements(posH,2); C=elements(posH,3); 37 p3=patch( 'Faces', [A B C], 'Vertices', [x y z],'CData', pattynenum,'FaceColor','flat','EdgeColor',edgecolor);37 p3=patch( 'Faces', [A B C], 'Vertices', [x y z],'CData',PattynEnum,'FaceColor','flat','EdgeColor',edgecolor); 38 38 %Stokes elements 39 39 alpha=0.35; 40 posS=find(data(:,2)== stokesenum);40 posS=find(data(:,2)==StokesEnum); 41 41 if ~isempty(posS) 42 42 A=elements(posS,1); B=elements(posS,2); C=elements(posS,3); 43 % p4=patch( 'Faces', [A B C], 'Vertices', [x y z],'CData', stokesenum,'FaceColor','flat','EdgeColor',edgecolor,'EdgeAlpha',alpha,'FaceAlpha',alpha);43 % p4=patch( 'Faces', [A B C], 'Vertices', [x y z],'CData',StokesEnum,'FaceColor','flat','EdgeColor',edgecolor,'EdgeAlpha',alpha,'FaceAlpha',alpha); 44 44 % legend([p1 p2 p3 p4],'Hutter''s elements','MacAyeal''s elements','Pattyn''s elements','Stokes''s elements'); 45 45 else … … 48 48 else 49 49 %Hutter elements 50 posH=find(data(:,1)== hutterenum);50 posH=find(data(:,1)==HutterEnum); 51 51 A=elements(posH,1); B=elements(posH,2); C=elements(posH,3); D=elements(posH,4); E=elements(posH,5); F=elements(posH,6); 52 p1=patch( 'Faces', [A B C], 'Vertices', [x y z],'CData', hutterenum,'FaceColor','flat','EdgeColor',edgecolor);53 patch( 'Faces', [D E F], 'Vertices', [x y z],'CData', hutterenum,'FaceColor','flat','EdgeColor',edgecolor);54 patch( 'Faces', [A B E D], 'Vertices', [x y z],'CData', hutterenum,'FaceColor','flat','EdgeColor',edgecolor);55 patch( 'Faces', [B E F C ], 'Vertices', [x y z],'CData', hutterenum,'FaceColor','flat','EdgeColor',edgecolor);56 patch( 'Faces', [C A D F ], 'Vertices', [x y z],'CData', hutterenum,'FaceColor','flat','EdgeColor',edgecolor);52 p1=patch( 'Faces', [A B C], 'Vertices', [x y z],'CData', HutterEnum,'FaceColor','flat','EdgeColor',edgecolor); 53 patch( 'Faces', [D E F], 'Vertices', [x y z],'CData', HutterEnum,'FaceColor','flat','EdgeColor',edgecolor); 54 patch( 'Faces', [A B E D], 'Vertices', [x y z],'CData', HutterEnum,'FaceColor','flat','EdgeColor',edgecolor); 55 patch( 'Faces', [B E F C ], 'Vertices', [x y z],'CData', HutterEnum,'FaceColor','flat','EdgeColor',edgecolor); 56 patch( 'Faces', [C A D F ], 'Vertices', [x y z],'CData', HutterEnum,'FaceColor','flat','EdgeColor',edgecolor); 57 57 %MacAyeal elements 58 posM=find(data(:,1)== macayealenum);58 posM=find(data(:,1)==MacAyealEnum); 59 59 A=elements(posM,1); B=elements(posM,2); C=elements(posM,3); D=elements(posM,4); E=elements(posM,5); F=elements(posM,6); 60 p2=patch( 'Faces', [A B C], 'Vertices', [x y z],'CData', macayealenum,'FaceColor','flat','EdgeColor',edgecolor);61 patch( 'Faces', [D E F], 'Vertices', [x y z],'CData', macayealenum,'FaceColor','flat','EdgeColor',edgecolor);62 patch( 'Faces', [A B E D], 'Vertices', [x y z],'CData', macayealenum,'FaceColor','flat','EdgeColor',edgecolor);63 patch( 'Faces', [B E F C ], 'Vertices', [x y z],'CData', macayealenum,'FaceColor','flat','EdgeColor',edgecolor);64 patch( 'Faces', [C A D F ], 'Vertices', [x y z],'CData', macayealenum,'FaceColor','flat','EdgeColor',edgecolor);60 p2=patch( 'Faces', [A B C], 'Vertices', [x y z],'CData', MacAyealEnum,'FaceColor','flat','EdgeColor',edgecolor); 61 patch( 'Faces', [D E F], 'Vertices', [x y z],'CData', MacAyealEnum,'FaceColor','flat','EdgeColor',edgecolor); 62 patch( 'Faces', [A B E D], 'Vertices', [x y z],'CData', MacAyealEnum,'FaceColor','flat','EdgeColor',edgecolor); 63 patch( 'Faces', [B E F C ], 'Vertices', [x y z],'CData', MacAyealEnum,'FaceColor','flat','EdgeColor',edgecolor); 64 patch( 'Faces', [C A D F ], 'Vertices', [x y z],'CData', MacAyealEnum,'FaceColor','flat','EdgeColor',edgecolor); 65 65 %Pattyn elements 66 posP=find(data(:,1)== pattynenum);66 posP=find(data(:,1)==PattynEnum); 67 67 A=elements(posP,1); B=elements(posP,2); C=elements(posP,3); D=elements(posP,4); E=elements(posP,5); F=elements(posP,6); 68 p3=patch( 'Faces', [A B C], 'Vertices', [x y z],'CData', pattynenum,'FaceColor','flat','EdgeColor',edgecolor);69 patch( 'Faces', [D E F], 'Vertices', [x y z],'CData', pattynenum,'FaceColor','flat','EdgeColor',edgecolor);70 patch( 'Faces', [A B E D], 'Vertices', [x y z],'CData', pattynenum,'FaceColor','flat','EdgeColor',edgecolor);71 patch( 'Faces', [B E F C ], 'Vertices', [x y z],'CData', pattynenum,'FaceColor','flat','EdgeColor',edgecolor);72 patch( 'Faces', [C A D F ], 'Vertices', [x y z],'CData', pattynenum,'FaceColor','flat','EdgeColor',edgecolor);68 p3=patch( 'Faces', [A B C], 'Vertices', [x y z],'CData', PattynEnum,'FaceColor','flat','EdgeColor',edgecolor); 69 patch( 'Faces', [D E F], 'Vertices', [x y z],'CData', PattynEnum,'FaceColor','flat','EdgeColor',edgecolor); 70 patch( 'Faces', [A B E D], 'Vertices', [x y z],'CData', PattynEnum,'FaceColor','flat','EdgeColor',edgecolor); 71 patch( 'Faces', [B E F C ], 'Vertices', [x y z],'CData', PattynEnum,'FaceColor','flat','EdgeColor',edgecolor); 72 patch( 'Faces', [C A D F ], 'Vertices', [x y z],'CData', PattynEnum,'FaceColor','flat','EdgeColor',edgecolor); 73 73 %Stokes elements 74 74 alpha=0.35; 75 posS=find(data(:,2)== stokesenum);75 posS=find(data(:,2)==StokesEnum); 76 76 if ~isempty(posS) 77 77 A=elements(posS,1); B=elements(posS,2); C=elements(posS,3); D=elements(posS,4); E=elements(posS,5); F=elements(posS,6); -
issm/trunk/src/m/classes/public/setelementstype.m
r1111 r1298 77 77 gridonhutter(md.elements(find(hutterflag),:))=1; 78 78 md.gridonhutter=gridonhutter; 79 md.elements_type(find(hutterflag),1)= hutterenum();79 md.elements_type(find(hutterflag),1)=HutterEnum(); 80 80 81 81 %2: MacAyeal elements … … 83 83 gridonmacayeal(md.elements(find(macayealflag),:))=1; 84 84 md.gridonmacayeal=gridonmacayeal; 85 md.elements_type(find(macayealflag),1)= macayealenum();85 md.elements_type(find(macayealflag),1)=MacAyealEnum(); 86 86 87 87 %3: Pattyn elements … … 89 89 gridonpattyn(md.elements(find(pattynflag),:))=1; 90 90 md.gridonpattyn=gridonpattyn; 91 md.elements_type(find(pattynflag),1)= pattynenum();91 md.elements_type(find(pattynflag),1)=PattynEnum(); 92 92 93 93 %4: Stokes elements 94 94 md.gridonstokes=gridonstokes; 95 md.elements_type(find(stokesflag),2)= stokesenum();95 md.elements_type(find(stokesflag),2)=StokesEnum(); 96 96 97 97 %5: None elements (non Stokes) 98 md.elements_type(find(~stokesflag),2)= noneenum();98 md.elements_type(find(~stokesflag),2)=NoneEnum(); 99 99 100 100 %Create the border grids between Pattyn and MacAyeal and extrude them … … 121 121 122 122 %figure out the border stokes grids 123 stokes_elements=find(md.elements_type(:,2)== stokesenum()); %find the elements on the stokes domain123 stokes_elements=find(md.elements_type(:,2)==StokesEnum()); %find the elements on the stokes domain 124 124 borderflags=zeros(md.numberofgrids,1); 125 125 borderflags(md.elements(stokes_elements,:))=1; %find all the grids of the elements on stokes domain, ie stokes grids and borderstokes … … 132 132 133 133 %figure out solution types 134 md.ishutter=double(any(md.elements_type(:,1)== hutterenum()));135 md.ismacayealpattyn=double(any(md.elements_type(:,1)== macayealenum() | md.elements_type(:,1)==pattynenum() ));134 md.ishutter=double(any(md.elements_type(:,1)==HutterEnum())); 135 md.ismacayealpattyn=double(any(md.elements_type(:,1)==MacAyealEnum() | md.elements_type(:,1)==PattynEnum() )); 136 136 md.isstokes=double(any(md.elements_type(:,2))); 137 137
Note:
See TracChangeset
for help on using the changeset viewer.