program main

c*********************************************************************72
c
cc MAIN is the main program for ANYPLT_PRB.
c
c  Discussion:
c
c    ANYPLT_PRB tests the ANYPLT library.
c
c  Licensing:
c
c    This code is distributed under the MIT license. 
c
c  Modified:
c
c    03 July 2008
c
c  Author:
c
c    John Burkardt
c
      implicit none

      character*80 carray
      integer icom
      integer iplt1
      integer iplt2
      integer ixplt1
      integer ixplt2
      integer iyplt1
      integer iyplt2
      integer marray
      real xplt1
      real xplt2
      real yplt1
      real yplt2

      common /anycom/ iplt1, iplt2, ixplt1, ixplt2, iyplt1, 
     &                iyplt2, marray, xplt1, xplt2, yplt1, yplt2
      common /anychr/ carray

      call timestamp ( )
      write ( *, '(a)' ) ' '
      write ( *, '(a)' ) 'ANYPLT_PRB'
      write ( *, '(a)' ) '  Test the ANYPLT library.'
c
c  Enable graphics, and select window portion.
c  If we're using the ANYCGM interface, request PS output.
c
      icom = 0

      xplt1 = 0.0E+00
      xplt2 = 1.0E+00
      yplt1 = 0.0E+00
      yplt2 = 1.0E+00
      iplt1 = 1

      carray = 'PS'

      call anyplt ( icom )
c
c  Get the version number
c
      icom = 13
      call anyplt ( icom )

      write ( *, '(a)' ) ' '
      write ( *, * ) carray

      call test01 ( )
      call test02 ( )
      call test03 ( )
c
c  End graphics
c
      icom = 1
      call anyplt ( icom )
c
c  Terminate.
c
      write ( *, '(a)' ) ' '
      write ( *, '(a)' ) 'ANYPLT_PRB'
      write ( *, '(a)' ) '  Normal end of ANYPLT tests.'
      write ( *, '(a)' ) ' '
      call timestamp ( )

      stop
      end
      subroutine test01 ( )

c*********************************************************************72
c
cc TEST01 draws a sine curve.
c
c  Licensing:
c
c    This code is distributed under the MIT license. 
c
c  Modified:
c
c    24 January 2009
c
c  Author:
c
c    John Burkardt
c
      implicit none

      character ( len = 80 ) carray
      integer i
      integer icom
      integer iplt1
      integer iplt2
      integer ixplt1
      integer ixplt2
      integer iyplt1
      integer iyplt2
      integer marray
      real xlen
      real xmax
      real xmin
      real xplt1
      real xplt2
      real ylen
      real ymax
      real ymin
      real yplt1
      real yplt2

      common /anycom/ iplt1, iplt2, ixplt1, ixplt2, iyplt1, 
     &                iyplt2, marray, xplt1, xplt2, yplt1, yplt2
      common /anychr/ carray

      write ( *, * ) ' '
      write ( *, * ) 'TEST01'
      write ( *, * ) '  Draw a sine curve.'
c
c  Begin plot
c
      icom = 2
      call anyplt ( icom )
c
c  Define user coordinate window
c
      xmin = -0.5E+00
      xlen = 11.0E+00
      xmax = xmin + xlen

      ymin = -1.2E+00
      ylen = 2.4E+00
      ymax = ymin + ylen

      icom = 3
      xplt1 = xmin
      xplt2 = xlen
      yplt1 = ymin
      yplt2 = ylen
      call anyplt ( icom )
c
c  Move to first point, and draw
c
      icom = 4
      do i = 1, 101
        xplt1 = real ( i - 1 ) / 10.0E+00
        yplt1 = sin ( xplt1 )
        call anyplt ( icom )
        icom = 5
      end do
c
c  Label every ten-th point
c
      carray = '.'
      do i = 1, 101, 10
        xplt1 = real ( i - 1 ) / 10.0E+00
        yplt1 = sin(xplt1)
        icom = 11
        call anyplt ( icom )
      end do
c
c  Print label
c
      icom = 7
      xplt1 = 0.10E+00
      xplt2 = 0.05E+00
      yplt1 = 0.85E+00
      yplt2 = 0.0
      carray = 'Sine curve'
      marray = 10
      call anyplt ( icom )
c
c  Label Y axis
c
      icom = 7
      xplt1 = 0.05E+00
      xplt2 = 0.025E+00
      yplt1 = 0.5E+00
      yplt2 = 90.0E+00
      carray = 'Y axis'
      marray = 6
      call anyplt ( icom )
c
c  Draw Y axis
c
      icom = 4
      xplt1 = 0.0E+00
      yplt1 = ymin
      call anyplt ( icom )
      icom = 5
      xplt1 = 0.0E+00
      yplt1 = ymax
      call anyplt ( icom )
c
c  Label X axis
c
      icom = 7
      xplt1 = 0.3E+00
      xplt2 = 0.025E+00
      yplt1 = 0.05E+00
      yplt2 = 0.0E+00
      carray = 'X axis'
      marray = 6
      call anyplt ( icom )
c
c  Draw X axis
c
      icom = 4
      xplt1 = xmin
      yplt1 = 0.0E+00
      call anyplt ( icom )
      icom = 5
      xplt1 = xmax
      yplt1 = 0.0E+00
      call anyplt ( icom )
c
c  End plot
c
      icom = 9
      call anyplt ( icom )

      return
      end
      subroutine test02 ( )

c*********************************************************************72
c
cc TEST02 ...
c
c  Licensing:
c
c    This code is distributed under the MIT license. 
c
c  Modified:
c
c    24 January 2009
c
c  Author:
c
c    John Burkardt
c
      implicit none

      real alpha
      character*80 carray
      real del
      real delt
      real dely
      integer i
      integer icom
      integer iplt1
      integer iplt2
      integer ixplt1
      integer ixplt2
      integer iyplt1
      integer iyplt2
      integer j
      integer marray
      integer nx
      integer ny
      real pi
      parameter ( pi = 3.141592653589793E+00 )
      real theta
      real tnorm
      real u1
      real u2
      real uu
      real v1
      real v2
      real vv
      real x1
      real x2
      real xlen
      real xmax
      real xmin
      real xplt1
      real xplt2
      real y1
      real y2
      real ylen
      real ymax
      real ymin
      real yplt1
      real yplt2

      common /anycom/ iplt1, iplt2, ixplt1, ixplt2, iyplt1, 
     &                iyplt2, marray, xplt1, xplt2, yplt1, yplt2
      common /anychr/ carray

      write ( *, * ) ' '
      write ( *, * ) 'TEST02'
      write ( *, * ) '  ???.'
c
c  Begin plot
c
      icom = 2
      call anyplt ( icom )
c
c  Define user coordinate window.
c  For this plot, make sure axes are of quite different lengthsc
c
      xmin = 0.0E+00
      xlen = 10.0E+00
      xmax = xmin + xlen
      ymin = 0.0E+00
      ylen = 1.0E+00
      ymax = ymin + ylen

      icom = 3
      xplt1 = xmin
      xplt2 = xlen
      yplt1 = ymin
      yplt2 = ylen
      call anyplt ( icom )
c
c  Move to first point, and draw
c
      ny = 5
      nx = 5
      delt = 0.5E+00 * ( xmax - xmin ) / real ( nx - 1 )
      dely = 0.5E+00 * ( ymax - ymin ) / real ( ny - 1 )

      carray = '.'

      do i = 1, ny
        do j = 1, nx

          x1 = xmin + real ( j - 1 ) * ( xmax - xmin ) / real ( nx - 1 )
          y1 = ymin + real ( i - 1 ) * ( ymax - ymin ) / real ( ny - 1 )

          xplt1 = x1
          yplt1 = y1
          icom = 4
          call anyplt ( icom )

          x2 = x1 + delt
          y2 = y1 + dely * exp ( x1 * y1 ) / exp ( xmax * ymax )
          xplt1 = x2
          yplt1 = y2
          icom = 5
          call anyplt ( icom )

          uu = x2 - x1
          vv = y2 - y1
          tnorm = sqrt ( uu * uu + vv * vv )

          if ( 0.0E+00 < tnorm )then
            theta = pi - 0.5 * atan ( 2.0E+00 )
            alpha = atan2 ( vv, uu )
            del = sqrt ( 5.0E+00 ) * tnorm / 3.0E+00
            u1 = x1 + del * cos ( alpha - theta )
            v1 = y1 + del * sin ( alpha - theta )
            u2 = x1 + del * cos ( alpha + theta )
            v2 = y1 + del * sin ( alpha + theta )

            xplt1 = u1
            yplt1 = v1
            icom = 4
            call anyplt ( icom )

            xplt1 = x2
            yplt1 = y2
            icom = 5
            call anyplt ( icom )

            xplt1 = u2
            yplt1 = v2
            icom = 5
            call anyplt ( icom )

          end if

          write ( *, * ) i,j

        end do
      end do
c
c  Print label
c
      icom = 7
      xplt1 = 0.10E+00
      xplt2 = 0.05E+00
      yplt1 = 0.85E+00
      yplt2 = 0.0E+00
      carray = 'Comparison of arrows'
      marray = 20
      call anyplt ( icom )
c
c  Label Y axis
c
      icom = 7
      xplt1 = 0.05E+00
      xplt2 = 0.025E+00
      yplt1 = 0.5E+00
      yplt2 = 90.0E+00
      carray = 'y axis'
      marray = 6
      call anyplt ( icom )
c
c  Draw Y axis
c
      icom = 4
      xplt1 = 0.0E+00
      yplt1 = ymin
      call anyplt ( icom )
      icom = 5
      xplt1 = 0.0E+00
      yplt1 = ymax
      call anyplt ( icom )
c
c  Label X axis
c
      icom = 7
      xplt1 = 0.3E+00
      xplt2 = 0.025E+00
      yplt1 = 0.05E+00
      yplt2 = 0.0E+00
      carray = 'X axis'
      marray = 6
      call anyplt ( icom )
c
c  Draw X axis
c
      icom = 4
      xplt1 = xmin
      yplt1 = 0.0E+00
      call anyplt ( icom )
      icom = 5
      xplt1 = xmax
      yplt1 = 0.0E+00
      call anyplt ( icom )
c
c  End plot
c
      icom = 9
      call anyplt ( icom )

      return
      end
      subroutine test03 ( )

c*********************************************************************72
c
cc TEST03 ...
c
c  Licensing:
c
c    This code is distributed under the MIT license. 
c
c  Modified:
c
c    24 January 2009
c
c  Author:
c
c    John Burkardt
c
      implicit none

      character*80 carray
      integer i
      integer icom
      integer iplt1
      integer iplt2
      integer ixplt1
      integer ixplt2
      integer iyplt1
      integer iyplt2
      integer marray
      real pi
      parameter ( pi = 3.14159265E+00 )
      real xlen
      real xmax
      real xmin
      real xplt1
      real xplt2
      real ylen
      real ymax
      real ymin
      real yplt1
      real yplt2

      common /anycom/ iplt1, iplt2, ixplt1, ixplt2, iyplt1, 
     &                iyplt2, marray, xplt1, xplt2, yplt1, yplt2
      common /anychr/ carray

      write ( *, * ) ' '
      write ( *, * ) 'TEST03'
      write ( *, * ) '  ???.'
c
c  Begin plot
c
      icom = 2
      call anyplt ( icom )
c
c  Define user coordinate window.
c  For this plot, make sure the axes are of quite different lengthsc
c
      xmin = 0.0E+00
      xlen = 1.0E+00
      xmax = xmin + xlen

      ymin = 0.0E+00
      ylen = 1.0E+00
      ymax = ymin + ylen

      icom = 3
      xplt1 = xmin
      xplt2 = xlen
      yplt1 = ymin
      yplt2 = ylen
      call anyplt ( icom )
c
c  Draw the arrows.
c
      icom = 14
      do i = 1, 12
        xplt1 = 0.5E+00
        yplt1 = 0.5E+00
        xplt2 = real ( i - 1 ) * 2.0E+00 * pi / 12.0E+00
        yplt2 = 0.5E+00 * real ( i ) / 12.0E+00
        call anyplt ( icom )
      end do
c
c  Print label
c
      write ( *, * ) 'DEBUG: TEST03, ICOM = 7'
      icom = 7
      xplt1 = 0.10E+00
      xplt2 = 0.03E+00
      yplt1 = 0.85E+00
      yplt2 = 0.0E+00
      carray = 'Equal angles despite scale?'
      marray = 27
      call anyplt ( icom )
c
c  Label Y axis
c
      icom = 7
      xplt1 = 0.05E+00
      xplt2 = 0.025E+00
      yplt1 = 0.5E+00
      yplt2 = 90.0E+00
      carray = 'Y axis'
      marray = 6
      call anyplt ( icom )
c
c  Draw Y axis
c
      icom = 4
      xplt1 = xmin
      yplt1 = ymin
      call anyplt ( icom )
      icom = 5
      xplt1 = xmin
      yplt1 = ymax
      call anyplt ( icom )
c
c  Label X axis
c
      icom = 7
      xplt1 = 0.3E+00
      xplt2 = 0.025E+00
      yplt1 = 0.05E+00
      yplt2 = 0.0E+00
      carray = 'X axis'
      marray = 6
      call anyplt ( icom )
c
c  Draw X axis
c
      icom = 4
      xplt1 = xmin
      yplt1 = ymin
      call anyplt ( icom )
      icom = 5
      xplt1 = xmax
      yplt1 = ymin
      call anyplt ( icom )
c
c  End the plot.
c
      icom = 9
      call anyplt ( icom )

      return
      end
      subroutine timestamp ( )

c*********************************************************************72
c
cc TIMESTAMP prints out the current YMDHMS date as a timestamp.
c
c  Licensing:
c
c    This code is distributed under the MIT license.
c
c  Modified:
c
c    12 January 2007
c
c  Author:
c
c    John Burkardt
c
c  Parameters:
c
c    None
c
      implicit none

      character * ( 8 ) ampm
      integer d
      character * ( 8 ) date
      integer h
      integer m
      integer mm
      character * ( 9 ) month(12)
      integer n
      integer s
      character * ( 10 ) time
      integer y

      save month

      data month /
     &  'January  ', 'February ', 'March    ', 'April    ', 
     &  'May      ', 'June     ', 'July     ', 'August   ', 
     &  'September', 'October  ', 'November ', 'December ' /

      call date_and_time ( date, time )

      read ( date, '(i4,i2,i2)' ) y, m, d
      read ( time, '(i2,i2,i2,1x,i3)' ) h, n, s, mm

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

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

      return
      end