Changeset 8823 for issm/trunk/src/m/utils/BC/SetMarineIceSheetBC.m

Timestamp:

07/06/11 14:30:00 (14 years ago)

Author:

Eric.Larour

Message:

Split spcvelocity into spcvx, spcvy and spcvz

File:

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

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

@@ -35 +35 @@
         warning('SetMarineIceSheetBC warning: ice front all around the glacier, no dirichlet found. Dirichlet must be added manually')
 end
-md.spcvelocity=zeros(md.numberofnodes,6);
-md.spcvelocity(pos,1:3)=1;
+md.spcvx=NaN*ones(md.numberofnodes,1);
+md.spcvy=NaN*ones(md.numberofnodes,1);
+md.spcvz=NaN*ones(md.numberofnodes,1);
+md.spcvx(pos)=0;
+md.spcvy(pos)=0;
+md.spcvz(pos)=0;
 md.diagnostic_ref=NaN*ones(md.numberofnodes,6);
 
@@ -42 +46 @@
 if (length(md.vx_obs)==md.numberofnodes & length(md.vy_obs)==md.numberofnodes)
         disp('      boundary conditions for diagnostic model: spc set as observed velocities');
-        md.spcvelocity(pos,4:5)=[md.vx_obs(pos) md.vy_obs(pos)]; %zeros for vz
+        md.spcvx(pos)=md.vx_obs(pos);
+        md.spcvy(pos)=md.vy_obs(pos);
 else
         disp('      boundary conditions for diagnostic model: spc set as zero');