Construct_r4fp64Stack Subroutine

  subroutine Construct_r4fp64Stack(stack,N,mold)
    class(r4fp64Stack),intent(out) :: stack
    integer,intent(in)            :: N
    real(real64),intent(in)      :: mold(:,:,:,:)
    ! local
    integer :: l(1:4),u(1:4)

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

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

  endsubroutine Construct_r4fp64Stack