module mixmap

  use types, only: &
       real64, int32

  implicit none

  save

  real(kind=real64), dimension(:), allocatable :: &
       star_y, star_sig
  integer(kind=int32) :: &
       star_n

  real(kind=real64), dimension(:,:), allocatable :: &
       model_data
  integer(kind=int32) :: &
       model_nel, model_nzone

contains

  subroutine set_stardata(y, sig, n)

    use types, only: &
         int32, real64

    implicit none

    real(kind=real64), intent(in), dimension(n) :: &
         y, sig
    integer(kind=int32), intent(in) :: &
         n

    if (allocated(star_y)) then
       deallocate(star_y, star_sig)
    endif

    allocate(star_y(n), star_sig(n))

    star_n = n
    star_y(:) = y(:)
    star_sig(:) = sig(:)

  end subroutine set_stardata

  !----------------------------------------------------------------------

  subroutine set_modeldata(data, nel, nzone)

    use types, only: &
         int32, real64

    implicit none

    real(kind=real64), intent(in), dimension(nel, nzone) :: &
         data
    integer(kind=int32), intent(in) :: &
         nel, nzone

    if (allocated(model_data)) then
       deallocate(model_data)
    endif

    allocate(model_data(nel, nzone))

    model_nel = nel
    model_nzone = nzone
    model_data(:,:) = data(:,:)

    print*, nel, ' elements and ', nzone, 'zones.'

  end subroutine set_modeldata

end module mixmap