0001 function md=Loadresults(md,fem,u_g);
0002
0003
0004
0005
0006
0007
0008
0009
0010 if strcmpi(md.type,'2d')
0011 if fem.ismacayealpattyn,
0012 m_dh=fem.m_dh;
0013 indx=1:6:m_dh.gridset.gsize;
0014 indy=2:6:m_dh.gridset.gsize;
0015 indp=4:6:m_dh.gridset.gsize;
0016 indx=indx(m_dh.tpart);
0017 indy=indy(m_dh.tpart);
0018 indp=indp(m_dh.tpart);
0019 else
0020 m_dhu=fem.m_dhu;
0021 indx=1:6:m_dhu.gridset.gsize;
0022 indy=2:6:m_dhu.gridset.gsize;
0023 indp=4:6:m_dhu.gridset.gsize;
0024 indx=indx(m_dhu.tpart);
0025 indy=indy(m_dhu.tpart);
0026 indp=indp(m_dhu.tpart);
0027 end
0028
0029 md.vx=u_g(indx)*md.yts;
0030 md.vy=u_g(indy)*md.yts;
0031 md.vel=sqrt(md.vx.^2+md.vy.^2);
0032 md.pressure=u_g(indp);
0033
0034 else
0035 m_dv=fem.m_dv;
0036 indx=1:6:m_dv.gridset.gsize;
0037 indy=2:6:m_dv.gridset.gsize;
0038 indz=3:6:m_dv.gridset.gsize;
0039 indp=4:6:m_dv.gridset.gsize;
0040 indx=indx(m_dv.tpart);
0041 indy=indy(m_dv.tpart);
0042 indz=indz(m_dv.tpart);
0043 indp=indp(m_dv.tpart);
0044
0045 md.vx=u_g(indx)*md.yts;
0046 md.vy=u_g(indy)*md.yts;
0047 md.vz=u_g(indz)*md.yts;
0048 md.vel=sqrt(md.vx.^2+md.vy.^2+md.vz.^2);
0049 md.pressure=u_g(indp);
0050 end