C+
C NAME:
C	pr_star_pole_cen_xy_dist_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	pr_star_pole_cen_xy_dist_st_rot(JJ,NSTAR,NOBS,XCENP,YCENP,XCEN,YCEN,alatyy,AROTA,ST)
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 pole on the original ASHI image
C     YCENP (NSTAR,NOBS) real        The Y location of the pole on the original ASHI image
C     XCEN  (NSTAR,NOBS) real        Distance of the pole to the lat, lon center on the image
C     YCEN  (NSTAR,NOBS) real        Distance of the pole to the X,Y center on the ASHI image
C     alatyy(NOBS)       real        The latitude (declination) of ASHI latitude and the location of the star    
C     AROTA (NSTAR,NOBS) real        Rotation of the original image to place the pole at the top of the new grid       
C     ST3   (NOBS)       real        The Sidereal Time of the Image
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        Distance of the pole to the X,Y center on the ASHI image
C     XCEN  (NSTAR,NOBS) real        Distance of the pole to the lat, lon center on the image
C     YCEN  (NSTAR,NOBS) real        The Y location of the lat-lon center on the original ASHI image
C     alatyy(NOBS)       real        The latitude (declination) of ASHI latitude and the location of the star    
C     AROTA (NSTAR,NOBS) real        Rotation of the original image to place the pole at the top of the new grid
C     AROTA (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 pr_star_pole_cen_xy_dist_st_rot(JJ,NSTAR,NOBS,XCENP,YCENP,XCEN,YCEN,alatyy,AROTA,ST)

        real            ST    (NOBS),
     &                  alatyy(NOBS)

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

          print*, ' '
          print*, ' Into pr_star_pole_cen_xy_dist_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) then
               print*,' Rot1 = Rotation to place pole at the top of the image'
               print*,' Rot2 = Rotation (HA - ST) from your center Pol placed at one location'
               write(*,'(A)'),'Star Image PXvalue PYvalue  D(la-Lo)  D(X,Y)  ASHI La   ST     Rot1    Rot2'
             end if
             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)
C            write(*,'(I2,I8,7F9.2)'),JJ,I,XCENP(JJ,I),YCENP(JJ,I),XCEN(JJ,I),YCEN(JJ,I),alatyy(I),-AROTA(JJ,I),ARPST
             write(*,'(I2,I6,F10.2,F8.2,2F9.2,4F8.2)')
     &             JJ,I,XCENP(JJ,I),YCENP(JJ,I),XCEN(JJ,I),YCEN(JJ,I),alatyy(I),ST(I),-AROTA(JJ,I),ARPST
C             if(I.eq.10) stop
            end do

            print*, ' Out of pr_star_pole_cen_xy_dist_st_rot'
	  return
	end