Index: /issm/trunk/src/m/classes/public/BasinConstrain.m
===================================================================
--- /issm/trunk/src/m/classes/public/BasinConstrain.m	(revision 1236)
+++ /issm/trunk/src/m/classes/public/BasinConstrain.m	(revision 1237)
@@ -1,14 +1,15 @@
 function md=BasinConstrain(md,domain);
-%BASINCONSTRAIN constrain basin, using a constraint domain outline, to dirichlet boundary conditions.
+%BASINCONSTRAIN - constrain basin
 %
-% usage: md=BasinConstrain(md,constraindomain)
-%          
-%            where constraindomain is an Argus domain outline file enclosing the geographical area of interest.
+%   Constrain basin using a constraint domain outline, 
+%   to dirichlet boundary conditions.
+%   constraindomain is an Argus domain outline file enclosing 
+%   the geographical area of interest.
 %
-% See also: N/A
-
-%prepare barycentre positions.
-xelem=md.x(md.elements)*[1;1;1]/3;
-yelem=md.y(md.elements)*[1;1;1]/3;
+%   Usage: 
+%      md=BasinConstrain(md,constraindomain)
+%
+%   Example:
+%      md=BasinConstrain(md,'DomainOutline.exp');
 
 %now, flag grids and elements outside the domain outline.
@@ -23,5 +24,4 @@
 md.elementonwater(elementnotondomain)=1;
 
-
 %now, make sure all elements on water have grids that are spc'd, otherwise, we'll get a singular problem.
 pos=find(md.elementonwater);
@@ -31,3 +31,2 @@
 md.dirichletvalues_diag(grids,1)=md.vx_obs(grids);
 md.dirichletvalues_diag(grids,2)=md.vy_obs(grids);
-