C+
C NAME:
C	print_star_pole_cen_xy_st_rot.f
C PURPOSE:
C	Print out the star and pole x_y location on the original ASHI image, the lat lon center on the original image and the 
C     rotation to get the star (usually Polaris) and the pole to the top of the new grid.
C
C CATEGORY:
C	I/O
C CALLING SEQUENCE:
C	print_star_pole_cen_xy_st_rot(JJ,NSTAR,NOBS,XCENP,YCENP,XCEN,YCEN,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     XCEN  (NSTAR,NOBS) real        The X location of the lat-lon center on the original ASHI image
C     YCEN  (NSTAR,NOBS) real        The Y location of the lat-lon center 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     AROTA (NSTAR,NOBS) real        Rotation of the original image to place the pole at the top of the new grid FROM MY CENTER        
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     XCEN  (NSTAR,NOBS) real        The X location of the lat-lon center on the original ASHI image
C     YCEN  (NSTAR,NOBS) real        The Y location of the lat-lon center 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     AROTA  (NSTAR,NOBS) + ST real   Rotation of the original image to place the pole at the top of the new grid and stars in same location 
C     AROTA2 (NSTAR,NOBS) + ST real   Rotation of the original image to place the pole at the top of the new grid FROM MY CENTER and stars in same location
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_cen_xy_st_rot(JJ,NSTAR,NOBS,XCENP,YCENP,XCEN,YCEN,AROTA,ST)

        real            xOPpos(NOBS),
     &                  yOPpos(NOBS)

        real            ST(NOBS)

        real            XCENP (NSTAR,NOBS),
     &                  YCENP (NSTAR,NOBS),
     &                  XCEN  (NSTAR,NOBS),
     &                  YCEN  (NSTAR,NOBS),
     &                  AROTA (NSTAR,NOBS)

          print*, ' '
          print*, ' Into print_star_pole_cen_xy_st_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 Rot = star+ST from your center'
             if(I.eq.1) print*,'Star Image  PXvalue  PYvalue   XCEN     YCEN    Star Rot'
             ARPST = -AROTA(JJ,I)-ST(I)
             if(ARPST.lt.0.0)   ARPST = ARPST + 360.0
             if(ARPST.gt.360.0) ARPST = ARPST - 360.0
C             write(*,'(I3,I8,7F9.2)'),JJ,I,XCENP(JJ,I),YCENP(JJ,I),XCEN(JJ,I),YCEN(JJ,I),ARPST,-AROTA(JJ,I),ST(I)
             write(*,'(I3,I8,7F9.2)'),JJ,I,XCENP(JJ,I),YCENP(JJ,I),XCEN(JJ,I),YCEN(JJ,I),ARPST
C             if(I.eq.10) stop
            end do

            print*, ' Out of print_star_pole_cen_xy_st_rot'
	  return
	end