Changeset 4962 for issm/trunk/src/m/utils/BC/SetIceShelfBC.m

Timestamp:

08/04/10 08:42:15 (15 years ago)

Author:

Mathieu Morlighem

Message:

Added constrained ice shelf

File:

• issm/trunk/src/m/utils/BC/SetIceShelfBC.m (modified) (3 diffs)

issm/trunk/src/m/utils/BC/SetIceShelfBC.m

@@ -1 @@
-function md=SetIceShelfBC(md,icefrontfile)
+function md=SetIceShelfBC(md,varargin)
 %SETICESHELFBC - Create the boundary conditions for diagnostic and thermal models for a  Ice Shelf with Ice Front
 %
@@ -7 +7 @@
 %
 %   Usage:
-%      md=SetIceShelfBC(md,icefrontfile)
+%      md=SetIceShelfBC(md,varargin)
 %
 %   Example:
+%      md=SetIceShelfBC(md);
 %      md=SetIceShelfBC(md,'Front.exp');
 %
@@ -15 +16 @@
 
 %grid on Dirichlet (boundary and ~icefront)
-if ~exist(icefrontfile)
-        error(['SetIceShelfBC error message: ice front file ' icefrontfile ' not found']);
+if nargin==2,
+        icefrontfile=varargin{1};
+        if ~exist(icefrontfile), error(['SetIceShelfBC error message: ice front file ' icefrontfile ' not found']); end
+        gridinsideicefront=ContourToMesh(md.elements,md.x,md.y,expread(icefrontfile,1),'node',2);
+        gridonicefront=double(md.gridonboundary & gridinsideicefront);
+elseif nargin==1,
+        gridonicefront=zeros(md.numberofgrids,1);
+else
+        help SetIceShelfBC
+        error('bad usage');
 end
-gridinsideicefront=ContourToMesh(md.elements,md.x,md.y,expread(icefrontfile,1),'node',2);
-gridonicefront=double(md.gridonboundary & gridinsideicefront);
 pos=find(md.gridonboundary & ~gridonicefront);
 md.spcvelocity=zeros(md.numberofgrids,6);