Changeset 10802
- Timestamp:
- 11/18/11 16:56:39 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/src/m/model/mesh/meshyams.m
r9734 r10802 79 79 field=sqrt(vx_obs.^2+vy_obs.^2); 80 80 81 %set nodeonwaterfield81 %set mask.vertexonwater field 82 82 if ~strcmp(groundeddomain,'N/A'), 83 83 nodeground=ContourToMesh(md.mesh.elements,md.mesh.x,md.mesh.y,groundeddomain,'node',2); 84 md. nodeonwater=ones(md.mesh.numberofvertices,1);85 md. nodeonwater(find(nodeground))=0;84 md.mask.vertexonwater=ones(md.mesh.numberofvertices,1); 85 md.mask.vertexonwater(find(nodeground))=0; 86 86 else 87 md. nodeonwater=zeros(md.mesh.numberofvertices,1);87 md.mask.vertexonwater=zeros(md.mesh.numberofvertices,1); 88 88 end 89 89 90 90 %adapt according to velocities 91 disp(' adapting...');91 disp(' Yo, Im adapting...'); 92 92 md=YamsCall(md,field,hmin,hmax,gradation(i),epsilon); 93 93 … … 121 121 if ~strcmp(groundeddomain,'N/A'), 122 122 nodeground=ContourToMesh(md.mesh.elements,md.mesh.x,md.mesh.y,groundeddomain,'node',2); 123 md. nodeonwater=ones(md.mesh.numberofvertices,1);124 md. nodeonwater(find(nodeground))=0;123 md.mask.vertexonwater=ones(md.mesh.numberofvertices,1); 124 md.mask.vertexonwater(find(nodeground))=0; 125 125 else 126 md. nodeonwater=zeros(md.mesh.numberofvertices,1);126 md.mask.vertexonwater=zeros(md.mesh.numberofvertices,1); 127 127 end 128 128 if strcmpi(Names.interp,'node'),
Note:
See TracChangeset
for help on using the changeset viewer.