Changeset 2824
- Timestamp:
- 01/13/10 16:20:45 (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/src/m/classes/public/project2d.m
r2687 r2824 1 function projection_value=project2d(md ,value,layer)1 function projection_value=project2d(md3d,value,layer) 2 2 %PROJECT2D - returns the value of a field for a given layer of the mesh 3 3 % … … 9 9 % 10 10 % Usage: 11 % projection_value=project2d(md ,value,layer)11 % projection_value=project2d(md3d,value,layer) 12 12 % 13 13 % Example: 14 % vel2=project2d(md ,md.vel,2);14 % vel2=project2d(md3d,md3d.vel,2); 15 15 16 16 %some checks on list of arguments … … 20 20 end 21 21 22 if md .counter<4,22 if md3d.counter<4, 23 23 error('this model has not been extruded! ... exiting'); 24 24 end 25 25 26 if ~strcmpi(md .type,'3d'),26 if ~strcmpi(md3d.type,'3d'), 27 27 error('wrong model type ... should be ''3d'''); 28 28 end 29 29 30 if ((layer<1) | (layer>md .numlayers)),31 error(['layer must be between 1 and ' num2str(md .numlayers)]);30 if ((layer<1) | (layer>md3d.numlayers)), 31 error(['layer must be between 1 and ' num2str(md3d.numlayers)]); 32 32 end 33 33 34 if size(value,1)==md .numberofgrids,35 projection_value=value((layer-1)*md .numberofgrids2d+1:layer*md.numberofgrids2d,:);34 if size(value,1)==md3d.numberofgrids, 35 projection_value=value((layer-1)*md3d.numberofgrids2d+1:layer*md3d.numberofgrids2d,:); 36 36 else 37 projection_value=value((layer-1)*md .numberofelements2d+1:layer*md.numberofelements2d,:);37 projection_value=value((layer-1)*md3d.numberofelements2d+1:layer*md3d.numberofelements2d,:); 38 38 end
Note:
See TracChangeset
for help on using the changeset viewer.