function obj = Beam(Envelope)
  function obj = Beam(Envelope, fDelta, fTheta, fPhi) 
  fDelta, fTheta, fPhi are with respect to the simulation coods 
  (z - along the propagation and along t, with x, y in the grid)
  default values are zero - corresponding to propagation along z.