Index: /issm/trunk/src/m/utils/BC/SetIceShelfBC.m
===================================================================
--- /issm/trunk/src/m/utils/BC/SetIceShelfBC.m	(revision 4961)
+++ /issm/trunk/src/m/utils/BC/SetIceShelfBC.m	(revision 4962)
@@ -1,3 +1,3 @@
-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 +7,8 @@
 %
 %   Usage:
-%      md=SetIceShelfBC(md,icefrontfile)
+%      md=SetIceShelfBC(md,varargin)
 %
 %   Example:
+%      md=SetIceShelfBC(md);
 %      md=SetIceShelfBC(md,'Front.exp');
 %
@@ -15,9 +16,15 @@
 
 %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);