matavs(eletypes,conn,xz,ELE_VAR,NODAL_VAR) matavs(eletypes,conn,xz,ELE_VAR) % No NODAL_VAR ELE_VAR = all element variables in a simple matrix NODAL_VAR = all nodal variables in a simple matrix matavs(eletypes,conn,xz,disarz(:,totalstep)) matavs(eletypes,conn,xz,disarz(:,150:158)) matavs(eletypes,conn,xz,[disarz(:,150:152) disarx(:,150:152)]) eletypes = 1; for 'quadrilateral', or 2 for triangle matavs(2,nodes,gcoord,-velocity,phi) matavs(1,nodes,gcoord,-velocity) matavs(2,nodes,gcoord,0,phi) matavs(2,nodes,gcoord,0,0) matavs(2,nodes,gcoord,velocity,phi) matavs(2,nodes,gcoord,velocity,0) matavs(2,nodes,gcoord,velocity,1) matavs(1,conn,xz,[porewater esarelez(:,150)],[disarx(:,100) disarz(:,200)]) matavs(1,nodes,gcoord,velocity,phi) -------------------------- Reza, 2004 conchange = 0; % need to change conn order? Put 1 for changing (later)