C
C FUNCTIONS/SUBROUTINES:
C
C MODIFICATION HISTORY:
C	May-2022, Bernard Jackson (UCSD)
C-
	subroutine read_centxy(cfilcxy,NOBS,XDIST,YDIST)

        character 	cfilcxy*26


        real         XDIST (NOBS),            ! mean X center position of the rotated Pole
     &               YSIST (NOBS)             ! mean Y center position of the rotated Pole

          print*,' '
          print*, ' Into read_centxy.f'
	    open (13, file=cfilcxy,status='old',recl=120,access='sequential',form='formatted',iostat=iReadashi)
	      if(iReadashi.ne.0) then
		 print *, 'There are no', cfilcxy,' files in this directory.' 
             close(13)
		 stop
	      end if
            NL = 0
	      do I=1,NOBS
	        read (13,'(I5,F12.5)',iostat=iReadashi) II,XDIST(I),YDIST(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)') II,XDIST(I),YDIST(I)
	      end do
            close(13)
	      print *, ' '
	      print *,' Reads complete, number of ASHI x y values read in =',NL
            print*, ' Out of read_centxy.f'
	  return
	end