procedure dunpack1(d1: integer; var u: ltuple; var ve,v1,v2: vector);
var
  st,te,t1: scalar;
  vt1,vt2: vector;
begin { dunpack1 }
  sabs(u[1],st);
  ismult(4,st);
  st.l:=st.u;
  squot(sone,st,te);
  squot(rho,u[2],t1);

  vtrans(d1,rho,u[3],v1,vt2);
  vscale(u[2],vt2,vt1);
  vprod(te,t1,ve,vt1,v2);
end { dunpack1 };

