C+
C NAME:
C	read_rot_rm_matc.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_rot_rm_matc(cfilrm,NOBS,doys8,st,rotrmm)
C
C INPUTS:
C	cfilrm      character	      Name of the ASHI input file 
C	NOBS		integer	      Number of image observations
C
C OUTPUTS:
C	st	      real	            Sidereal time
C	rotrmm	real	            Rotation to get the pole to the top of the image
C
C FUNCTIONS/SUBROUTINES:
C
C MODIFICATION HISTORY:
C	November-2022, Bernard Jackson (UCSD)
C-
	subroutine read_rot_rm_matc(cfilrm,NOBS,doys8,st,rotrmm)
        character 	cfilrm*29

        character    cSRC(9)*78

        real         st     (NOBS),             ! Sidereal time for each measurement
     &               rotrmm (NOBS)              ! mean rotation from Matthew's center

        real*8       doys8  (NOBS)              ! Julian date of observation      
        
          print*,' '
          print*, ' Into read_rot_rm_matc.f',iyr,idoyst
	    open (13, file=cfilrm,status='old',recl=120,access='sequential',form='formatted',iostat=iReadashi)
	      if(iReadashi.ne.0) then
		 print *, 'There are no', cfilrm,' files in this directory.' 
             close(13)
		 stop
	      end if
	      NL = 1
            iprint = 0
C            iprint = 1
            do I = 1,9
              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*,'All the above 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,'(I4,F14.5,10X,F8.3,13X,F7.2)',iostat=iReadashi) II,doys8(I),st(I),rotrmm(I)
	        if(iReadashi.eq.0) then
                NL = NL + 1
              else
                close(13)
	          print *, ' '
	          print *, ' Reads not complete, number of ASHI track values read in =',NL
                stop
	        end if
              if(NL.le.30.and.iprint.eq.1)  write (*,'(2I5,F14.5,F8.3,F7.2)') I,II,doys8(I),st(I),rotrmm(I)
	      end do
            close(13)
	      print *, ' '
	      print *,' Reads complete, number of ASHI track values read in =',NL
            print*, ' Out of read_rot_rm_matc.f',iyr,idoyst
	  return
	end