0001 function u_g=CieloVelocityExtrude(md,u_g);
0002
0003
0004 if strcmpi(md.type,'2d'),
0005
0006 else
0007
0008
0009
0010 grids2d=md.proj2d_coll_grids(1:md.number_grids2d_coll);
0011 vx2d=u_g(1:6:6*md.number_grids2d_coll);
0012 vy2d=u_g(2:6:6*md.number_grids2d_coll);
0013
0014
0015 vx=zeros(md.numberofgrids2d,1);
0016 vx(grids2d)=vx2d;
0017
0018 vy=zeros(md.numberofgrids2d,1);
0019 vy(grids2d)=vy2d;
0020
0021
0022 vx3d=project3d(md,vx,'node');
0023 vx3d(1:md.number_grids2d_coll)=0;
0024 vy3d=project3d(md,vy,'node');
0025 vy3d(1:md.number_grids2d_coll)=0;
0026
0027
0028 for i=1:length(md.penalties),
0029 vx3d(md.penalties(i).grids3d_noncollapsed)=0;
0030 vy3d(md.penalties(i).grids3d_noncollapsed)=0;
0031 end
0032
0033
0034 u_g(1:6:end)=u_g(1:6:end)+vx3d;
0035 u_g(2:6:end)=u_g(2:6:end)+vy3d;
0036
0037 end
0038
0039 end