Changeset 1237
- Timestamp:
- 07/06/09 11:20:57 (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/src/m/classes/public/BasinConstrain.m
r1234 r1237 1 1 function md=BasinConstrain(md,domain); 2 %BASINCONSTRAIN constrain basin, using a constraint domain outline, to dirichlet boundary conditions.2 %BASINCONSTRAIN - constrain basin 3 3 % 4 % usage: md=BasinConstrain(md,constraindomain) 5 % 6 % where constraindomain is an Argus domain outline file enclosing the geographical area of interest. 4 % Constrain basin using a constraint domain outline, 5 % to dirichlet boundary conditions. 6 % constraindomain is an Argus domain outline file enclosing 7 % the geographical area of interest. 7 8 % 8 % See also: N/A9 10 % prepare barycentre positions.11 xelem=md.x(md.elements)*[1;1;1]/3; 12 yelem=md.y(md.elements)*[1;1;1]/3;9 % Usage: 10 % md=BasinConstrain(md,constraindomain) 11 % 12 % Example: 13 % md=BasinConstrain(md,'DomainOutline.exp'); 13 14 14 15 %now, flag grids and elements outside the domain outline. … … 23 24 md.elementonwater(elementnotondomain)=1; 24 25 25 26 26 %now, make sure all elements on water have grids that are spc'd, otherwise, we'll get a singular problem. 27 27 pos=find(md.elementonwater); … … 31 31 md.dirichletvalues_diag(grids,1)=md.vx_obs(grids); 32 32 md.dirichletvalues_diag(grids,2)=md.vy_obs(grids); 33
Note:
See TracChangeset
for help on using the changeset viewer.