Construct_r3fp64Stack Subroutine

  subroutine Construct_r3fp64Stack(stack,N,mold)
    class(r3fp64Stack),intent(out) :: stack
    integer,intent(in)            :: N
    real(real64),intent(in)      :: mold(:,:,:)
    ! local
    integer :: l(1:3),u(1:3)

    l = lbound(mold)
    u = ubound(mold)

    allocate(stack%tokens(l(1):u(1), &
                          l(2):u(2), &
                          l(3):u(3), &
                          1:N))
    stack%top_index = 0

  endsubroutine Construct_r3fp64Stack