program main

!*****************************************************************************80
!
!! tsp_display_test() tests tsp_display().
!
!  Licensing:
!
!    This code is distributed under the MIT license.
!
!  Modified:
!
!    21 July 2026
!
!  Author:
!
!    John Burkardt
!
  implicit none

  integer, parameter :: rk8 = kind ( 1.0D+00 )

  integer, parameter :: n = 42 

  integer i
  integer i2
  integer, save, dimension ( n ) :: order = (/ &
    42, 41, 40, 39, 38, 37, 36, 35, 34, 33, &
    32, 31, 30, 29, 28, 27, 26, 25, 24, 23, &
    22, 21, 20, 19, 18, 17, 16, 15, 14, 13, &
    12, 11, 10,  9,  8,  7,  6,  5,  4,  3, &
     2,  1 /)
!
!  Stupid Fortran requires me to COUNT the size of a string.
!
  character ( len = 9 ) :: prefix = 'dantzig42'
!
!  Awkward Fortran convention means we list matrix data as though it were
!  (nx2) row-wise, but it will be stored as (2xn), so we have to transpose it.
!
  real ( kind = rk8 ), save, dimension ( n, 2 ) :: x = transpose ( reshape ( (/ &
      170.0,  85.0, &
      166.0,  88.0, &
      133.0,  73.0, &
      140.0,  70.0, &
      142.0,  55.0, &
      126.0,  53.0, &
      125.0,  60.0, &
      119.0,  68.0, &
      117.0,  74.0, &
       99.0,  83.0, &
       73.0,  79.0, &
       72.0,  91.0, &
       37.0,  94.0, &
        6.0, 106.0, &
        3.0,  97.0, &
       21.0,  82.0, &
       33.0,  67.0, &
        4.0,  66.0, &
        3.0,  42.0, &
       27.0,  33.0, &
       52.0,  41.0, &
       57.0,  59.0, &
       58.0,  66.0, &
       88.0,  65.0, &
       99.0,  67.0, &
       95.0,  55.0, &
       89.0,  55.0, &
       83.0,  38.0, &
       85.0,  25.0, &
      104.0,  35.0, &
      112.0,  37.0, &
      112.0,  24.0, &
      113.0,  13.0, &
      125.0,  30.0, &
      135.0,  32.0, &
      147.0,  18.0, &
      147.0,  36.0, &
      154.0,  45.0, &
      157.0,  54.0, &
      158.0,  61.0, &
      172.0,  82.0, &
      174.0,  87.0 /), (/ 2, n /) ) )
  real ( kind = rk8 ), allocatable :: x2(:,:)

  call timestamp ( )

  write ( *, '(a)' ) ''
  write ( *, '(a)' ) 'tsp_display_test():'
  write ( *, '(a)' ) '  Fortran90 version.'
  write ( *, '(a)' ) '  Display the cities and route solving a traveling'
  write ( *, '(a)' ) '  salesman problem (TSP).'
  write ( *, '(a)' ) '  Dataset is "' // prefix // '"'
!
!  Create the circuit data.
!
  allocate ( x2 ( 1 : n, 1 : 2 ) )

  do i = 1, n
    i2 = order(i)
    x2(i,1:2) = x(i2,1:2)
  end do
!
!  Send the data to the plotter.
!
  call tsp_display ( prefix, n, x2 )

  deallocate ( x2 )
!
!  Terminate.
!
  write ( *, '(a)' ) ''
  write ( *, '(a)' ) 'tsp_display_test():'
  write ( *, '(a)' ) '  Normal end of execution.'

  call timestamp ( )
 
  stop
end
subroutine timestamp ( )

!*****************************************************************************80
!
!! timestamp() prints the current YMDHMS date as a time stamp.
!
!  Example:
!
!    31 May 2001   9:45:54.872 AM
!
!  Licensing:
!
!    This code is distributed under the MIT license.
!
!  Modified:
!
!    15 August 2021
!
!  Author:
!
!    John Burkardt
!
  implicit none

  character ( len = 8 ) ampm
  integer d
  integer h
  integer m
  integer mm
  character ( len = 9 ), parameter, dimension(12) :: month = (/ &
    'January  ', 'February ', 'March    ', 'April    ', &
    'May      ', 'June     ', 'July     ', 'August   ', &
    'September', 'October  ', 'November ', 'December ' /)
  integer n
  integer s
  integer values(8)
  integer y

  call date_and_time ( values = values )

  y = values(1)
  m = values(2)
  d = values(3)
  h = values(5)
  n = values(6)
  s = values(7)
  mm = values(8)

  if ( h < 12 ) then
    ampm = 'AM'
  else if ( h == 12 ) then
    if ( n == 0 .and. s == 0 ) then
      ampm = 'Noon'
    else
      ampm = 'PM'
    end if
  else
    h = h - 12
    if ( h < 12 ) then
      ampm = 'PM'
    else if ( h == 12 ) then
      if ( n == 0 .and. s == 0 ) then
        ampm = 'Midnight'
      else
        ampm = 'AM'
      end if
    end if
  end if

  write ( *, '(i2.2,1x,a,1x,i4,2x,i2,a1,i2.2,a1,i2.2,a1,i3.3,1x,a)' ) &
    d, trim ( month(m) ), y, h, ':', n, ':', s, '.', mm, trim ( ampm )

  return
end