0001 function u_g=VelocityExtrude(md,u_g);
0002
0003
0004
0005
0006
0007
0008
0009 if strcmpi(md.type,'2d'),
0010
0011 else
0012
0013
0014
0015
0016 grids2d=md.gridonmacayeal(1:md.numberofgrids2d) & ~(md.gridonhutter(1:md.numberofgrids2d) | md.gridonpattyn(1:md.numberofgrids2d) | md.gridondirichlet_diag(1:md.numberofgrids2d));
0017
0018 vx2d=u_g(1:6:6*md.numberofgrids2d).*grids2d;
0019 vy2d=u_g(2:6:6*md.numberofgrids2d).*grids2d;
0020
0021
0022 vx3d=project3d(md,vx2d,'node');
0023 vy3d=project3d(md,vy2d,'node');
0024
0025
0026 vx3d(1:md.numberofgrids2d)=0;
0027 vy3d(1:md.numberofgrids2d)=0;
0028
0029
0030 u_g(1:6:end)=u_g(1:6:end)+vx3d;
0031 u_g(2:6:end)=u_g(2:6:end)+vy3d;
0032
0033 end
0034
0035 end