Changeset 7967

Timestamp:

04/19/11 09:32:06 (14 years ago)

Author:

seroussi

Message:

mixed up hydrology and velocity boundary conditions in SetMariceIceSheet

File:

• issm/trunk/src/m/utils/BC/SetMarineIceSheetBC.m (modified) (2 diffs)

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

@@ -38 +38 @@
 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 +45 @@
         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)