real(4), allocatable :: x(:)
  real(wp) , allocatable::y(:,:)