C+
C NAME:
C	read_cen_rm_bvjc.f
C PURPOSE:
C	Read the sidereal time and the mean rotation determined from Polaris and Matthew's center to place the Pole at the map top 
C CATEGORY:
C	I
C CALLING SEQUENCE:
C	call read_cen_rm_bvjc(cfilrb,NOBS,centxp,centyp,centx,centy)
C
C INPUTS:
C	cfilrb      character	      Name of the ASHI input file 
C	NOBS		integer	      Number of image observations
C
C OUTPUTS:
C	centxp	real	            Unrotated pole mean X centroid position
C	centyp	real	            Unrotated pole mean Y centroid position
C	centx	      real	            Unrotated BVJ mean X center location
C	centy	      real	            Unrotated BVJ mean Y center location    
C
C FUNCTIONS/SUBROUTINES:
C
C MODIFICATION HISTORY:
C	November-2022, Bernard Jackson (UCSD)
C-
	subroutine read_cen_rm_bvjc(cfilrb,NOBS,centxp,centyp,centx,centy)
        character 	cfilrb*26

        character    cSRC(1)*78

        real         centxp (NOBS),             ! mean X position of the unrotated Pole
     &               centyp (NOBS),             ! mean Y position of the unrotated Pole
     &               centx  (NOBS),             ! mean X position of Bernie's unrotated center
     &               centy  (NOBS)              ! mean Y position of Bernie's unrotated center

          print*,' '
          print*, ' Into read_cen_rm_bvjc.f'
	    open (13, file=cfilrb,status='old',recl=120,access='sequential',form='formatted',iostat=iReadashi)
	      if(iReadashi.ne.0) then
		 print *, 'There are no', cfilrb,' files in this directory.' 
             close(13)
		 stop
	      end if
	      NL = 1
            iprint = 0
C            iprint = 1
            do I = 1,2
              read (13,'(A78)',iostat=iReadashi) cSRC(I)
 	        if(iReadashi.eq.0) then
                if(iprint.eq.1) then
                  if(I.eq.1) print*,' '
 	            print*, cSRC(I)
                end if
 	        end if
              if(I.eq.9.and.iprint.eq.1) then
                print*,' '
                print*,'The above header information is skipped.' 
                print*,'Only the material below in the file is read: ST and Rotrm'
                print*,' '
              end if
            end do
         
            NL = 0
	      do I=1,NOBS
	        read (13,'(7X,I4,16X,4F8.2)',iostat=iReadashi) II,centxp(I),centyp(I),centx(I),centy(I)
	        if(iReadashi.eq.0) then
                NL = NL + 1
              else
                close(13)
	          print *, ' '
	          print *, ' Reads not complete, number of ASHI observation values read in =',NL
                stop
	        end if
              if(NL.le.10.and.iprint.eq.1)  write (*,'(2I5,4F8.2)') I,II,centxp(I),centyp(I),centx(I),centy(I)
	      end do
            close(13)
	      print *, ' '
	      print *,' Reads complete, number of ASHI observation values read in =',NL
            print*, ' Out of Into read_cen_rm_bvjc.f'
	  return
	end