Changeset 11315
- Timestamp:
- 02/03/12 12:08:36 (13 years ago)
- Location:
- issm/trunk-jpl/src/m/model
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
TabularUnified issm/trunk-jpl/src/m/model/extrude.m ¶
r11234 r11315 194 194 for i=2:numlayers-1, 195 195 md.mesh.elementconnectivity((i-1)*md.mesh.numberofelements2d+1:(i)*md.mesh.numberofelements2d,:)... 196 =md.mesh.elementconnectivity((i-1)*md.mesh.numberofelements2d+1:(i)*md.mesh.numberofelements2d,:)+md.mesh.numberofelements2d;196 =md.mesh.elementconnectivity((i-1)*md.mesh.numberofelements2d+1:(i)*md.mesh.numberofelements2d,:)+md.mesh.numberofelements2d; 197 197 end 198 198 md.mesh.elementconnectivity(find(isnan(md.mesh.elementconnectivity)))=0; … … 216 216 md.mask.vertexonwater=project3d(md,'vector',md.mask.vertexonwater,'type','node'); 217 217 if ~isnan(md.inversion.cost_functions_coefficients),md.inversion.cost_functions_coefficients=project3d(md,'vector',md.inversion.cost_functions_coefficients,'type','node');end; 218 219 %Put lithostatic pressure is there is an existing pressure 218 if ~isnan(md.qmu.partition),md.qmu.partition=project3d(md,'vector',md.qmu.partition','type','node');end 219 220 %Put lithostatic pressure if there is an existing pressure 220 221 if ~isnan(md.initialization.pressure), 221 222 md.initialization.pressure=md.constants.g*md.materials.rho_ice*(md.geometry.surface-md.mesh.z); -
TabularUnified issm/trunk-jpl/src/m/model/partition/partitioner.m ¶
r10232 r11315 99 99 %extrude if we are in 3D: 100 100 if md.mesh.dimension==3, 101 md3d.qmu.vertex_weight=md.qmu.vertex_weight; 102 md3d.qmu.adjacency=md.qmu.adjacency; 101 103 md=md3d; 102 104 part=project3d(md,'vector',part','type','node');
Note:
See TracChangeset
for help on using the changeset viewer.