Index: /issm/trunk/src/m/utils/BC/SetMarineIceSheetBC.m
===================================================================
--- /issm/trunk/src/m/utils/BC/SetMarineIceSheetBC.m	(revision 7966)
+++ /issm/trunk/src/m/utils/BC/SetMarineIceSheetBC.m	(revision 7967)
@@ -38,8 +38,4 @@
 md.spcvelocity(pos,1:3)=1;
 
-md.spcwatercolumn=zeros(md.numberofgrids,2);
-pos=find(md.gridonboundary); 
-md.spcwatercolumn(pos,1)=1;
-
 %Dirichlet Values
 if (length(md.vx_obs)==md.numberofgrids & length(md.vy_obs)==md.numberofgrids)
@@ -49,4 +45,8 @@
 	disp('      boundary conditions for diagnostic model: spc set as zero');
 end
+
+md.spcwatercolumn=zeros(md.numberofgrids,2);
+pos=find(md.gridonboundary); 
+md.spcwatercolumn(pos,1)=1;
 
 %segment on Neumann (Ice Front)