pdf_push_sigma_2d Subroutine

subroutine pdf_push_sigma_2d(self,sigma)
  type(pdf_kernel_2d) :: self    
  real(8)             :: sigma(2,2)
  if(.not.self%status)stop "PDF_PUSH_SIGMA: PDF not allocated"
  self%sigma = sigma
  self%variance=.true.
end subroutine pdf_push_sigma_2d