C+
C NAME:
C	nsew_max.f
C PURPOSE:
C	print_the North, South, East, and West maxima of the data images and rotation value of each ASHI image. This analysis attempts to find
C     a graohical solution for the x and y scale from the star (Polaris in the first attempt) to from the distance across the image from
C     the low to high values.
C CATEGORY:
C	I/O
C CALLING SEQUENCE:
C	nsew_max(JJ,NSTAR,NOBS,xOPpos,distcx,yOPpos,distcy,XCENP,YCENP,ROTA)
C
C INPUTS:
C     JJ                 integer     Stellar value number 1 for Polaris, 2 for ????, etc
C     NSTAR              integer     Total number of stars possible
C     NOBS               integer     Number of star images used
C     xOPpos(NSTAR,NOBS) real        The maximum X-position of the pole on the original image from the ASHI latitude and the location of the star
C     yOPpos(NSTAR,NOBS) real        The maximum Y-position of the pole on the original image from the ASHI latitude and the location of the star
C     XCENP (NSTAR,NOBS) real        The X distance of the pole from the ASHI latitude in degrees
C     YCENP (NSTAR,NOBS) real        The Y distance of the pole from the ASHI latitude in degrees
C     ROTA  (NSTAR,NOBS) real        Rotation of the original image to the Polaris position       
C
C OUTPUTS:
C     JJ                 integer     Stellar value number 1 for Polaris, 2 for ????, etc
C     I                  integer     The image number measured
C     xOPpos(NSTAR,NOBS) real        The maximum X-position of the pole on the original image from the ASHI latitude and the location of the star
C     yOPpos(NSTAR,NOBS) real        The maximum Y-position of the pole on the original image from the ASHI latitude and the location of the star
C     XCENP (NSTAR,NOBS) real        The X distance of the pole from the ASHI latitude in degrees
C     YCENP (NSTAR,NOBS) real        The Y distance of the pole from the ASHI latitude in degrees
C     ROTA  (NSTAR,NOBS) real        Rotation of the original image to the Polaris position   
C                                   (These outputs are printed so  they can be plotted on a spreadsheet
C
C FUNCTIONS/SUBROUTINES:
C
C MODIFICATION HISTORY:
C	September-2022, Bernard Jackson (UCSD)
C-
	subroutine nsew_max(JJ,NSTAR,NOBS,xOPpos,yOPpos,XCENP,YCENP,ROTA)

        real            xOPpos(NSTAR,NOBS),
     &                  yOPpos(NSTAR,NOBS)

        real            XCENP (NSTAR,NOBS),
     &                  YCENP (NSTAR,NOBS),
     &                  ROTA  (NSTAR,NOBS)

          print*, ' '
          print*, ' Into print_distcdiff'
           
C        Print out max and min and the image rotation to put the north pole at the top in the original image           

            do I=1,NOBS

              if(ROTA(JJ,I).gt.355.0.or.ROTA(JJ,I).lt.5.0) then
                write(*,'(A,I2,I5,F7.1,F8.3,45X,F7.2)')'North',JJ,I,yOPpos(JJ,I),XCENP(JJ,I),ROTA(JJ,I)
              end if
              if(ROTA(JJ,I).gt.175.0.and.ROTA(JJ,I).lt.185.0) then
                write(*,'(A,I2,I5,15X,F7.1,F8.3,30X,F7.2)')'South',JJ,I,yOPpos(JJ,I),-YCENP(JJ,I),ROTA(JJ,I)
              end if
              if(ROTA(JJ,I).gt.265.0.and.ROTA(JJ,I).lt.275.0) then
                write(*,'(A,I2,I5,30X,F7.1,F8.3,15X,F7.2)')'East ',JJ,I,xOPpos(JJ,I),-XCENP(JJ,I),ROTA(JJ,I)
              end if
              if(ROTA(JJ,I).gt.85.0.and.ROTA(JJ,I).lt.95.0) then
                write(*,'(A,I2,I5,45X,F7.1,F8.3,F7.2)')'West ',JJ,I,xOPpos(JJ,I),YCENP(JJ,I),ROTA(JJ,I)
              end if
            end do

            print*, ' Out of print_cent_rot'
	  return
	end