Loadresults

PURPOSE ^

LOADRESULTS - load results onto model

SYNOPSIS ^

function md=Loadresults(md,fem,u_g);

DESCRIPTION ^

LOADRESULTS - load results onto model

   The solution (vel,pressure,...) are on the g-set,
   use index to get the value on each node

   Usage:
      md=Loadresults(md,fem,u_g)

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

0001 function md=Loadresults(md,fem,u_g);
0002 %LOADRESULTS - load results onto model
0003 %
0004 %   The solution (vel,pressure,...) are on the g-set,
0005 %   use index to get the value on each node
0006 %
0007 %   Usage:
0008 %      md=Loadresults(md,fem,u_g)
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

Generated on Sun 29-Mar-2009 20:22:55 by m2html © 2003