C+
C NAME:
C	print_star_pole_xy_rot.f
C PURPOSE:
C	Print out the star and pole x_y location on the original ASHI image and the rotation to get the pole to the top of the new grid.,
C
C CATEGORY:
C	I/O
C CALLING SEQUENCE:
C	print_star_pole_xy_rot(JJ,NSTAR,NOBS,XCENP,YCENP,xOPpos,yOPpos,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     XCENP (NSTAR,NOBS) real        The X location of the star on the original ASHI image
C     YCENP (NSTAR,NOBS) real        The Y location of the star on the original ASHI image
C     xOPpos(NSTAR,NOBS) real        The 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 Y-position of the pole on the original image from the ASHI latitude and the location of the star
C     ROTA  (NSTAR,NOBS) real        Rotation of the original image to place the pole at the top of the new grid       
C
C OUTPUTS:
C     JJ                 integer     Stellar value number 1 for Polaris, 2 for ????, etc
C     I                  integer     The image number measured
C     XCENP (NSTAR,NOBS) real        The X location of the star on the original ASHI image
C     YCENP (NSTAR,NOBS) real        The Y location of the star on the original ASHI image
C     xOPpos(NSTAR,NOBS) real        The 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 Y-position of the pole on the original image from the ASHI latitude and the location of the star
C     ROTA  (NSTAR,NOBS) real        Rotation of the original image to place the pole at the top of the new grid  
C                                   (These outputs are printed so  they can be plotted on a spreadsheet
C
C FUNCTIONS/SUBROUTINES:
C
C MODIFICATION HISTORY:
C	October-2022, Bernard Jackson (UCSD)
C-
	subroutine print_star_pole_xy_rot(JJ,NSTAR,NOBS,XCENP,YCENP,xOPpos,yOPpos,ROTA,AROTA)

        real            xOPpos(NOBS),
     &                  yOPpos(NOBS)

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

          print*, ' '
          print*, ' Into print_star_pole_x-y_rot'
           
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(I.eq.1) print*,'Star Image   Xval    Yval  PXvalue PYvalue Star Rot Pole Rot'
             write(*,'(I3,I8,6F8.2)')JJ,I,xOPpos(I),yOPpos(I),XCENP(JJ,I),YCENP(JJ,I),ROTA(JJ,I),-AROTA(JJ,I)
C             if(I.eq.10) stop

            end do

            print*, ' Out of print_star_pole_x-y_rot'
	  return
	end