Changeset 1237


Ignore:
Timestamp:
07/06/09 11:20:57 (16 years ago)
Author:
Mathieu Morlighem
Message:

minor

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk/src/m/classes/public/BasinConstrain.m

    r1234 r1237  
    11function md=BasinConstrain(md,domain);
    2 %BASINCONSTRAIN constrain basin, using a constraint domain outline, to dirichlet boundary conditions.
     2%BASINCONSTRAIN - constrain basin
    33%
    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.
    78%
    8 % See also: N/A
    9 
    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');
    1314
    1415%now, flag grids and elements outside the domain outline.
     
    2324md.elementonwater(elementnotondomain)=1;
    2425
    25 
    2626%now, make sure all elements on water have grids that are spc'd, otherwise, we'll get a singular problem.
    2727pos=find(md.elementonwater);
     
    3131md.dirichletvalues_diag(grids,1)=md.vx_obs(grids);
    3232md.dirichletvalues_diag(grids,2)=md.vy_obs(grids);
    33 
Note: See TracChangeset for help on using the changeset viewer.