CI  Read raw real data and convert to complex data for calculation.
	PROGRAM SKYREF
      COMPLEX TARGET(6,10,4,21)
      REAL RAW(42),kgateleft,kgateright
      INTEGER NPTS,NTRACE
	CHARACTER STEM*11,FILENAME(6)*12,dscrp*4,outfile(6)*12,TYPE*1
	character ang*1
CI               
      NA=6                            
	print*,'Enter stem (first 11 letters) of skycal files:'
	read(*,'(a11)')stem    
C
	print*,'Enter 4-letter descriptor for output files:'
	read(*,'(a4)')dscrp
c
c

      OPEN(1,FILE=stem//'2')
      OPEN(2,FILE=stem//'3')
      OPEN(3,FILE=stem//'4')
      OPEN(4,FILE=stem//'5')
      OPEN(5,FILE=stem//'6')
      OPEN(6,FILE=stem//'7')
      NPTS = 21
	DO 5 I=1,NA
      READ(I,*) NTRACE,NFREQ,KGATELEFT,KGATERIGHT
c2000    format(1x,i1,13x,i1,13x,f5.1,9x,f5.1)
C
5       CONTINUE
      PRINT*,NTRACE,NFREQ,KGATELEFT,KGATERIGHT
         DO 12 IA=1,NA
	  do 14, ip = 1,4
            READ(IA,*) (RAW(K),K=1,NPTS*2)
            DO 14 J=1,NPTS
               J2=J*2-1
   14          TARGET(IA,1,IP,J)=CMPLX(RAW(J2),RAW(J2+1))
   12    CONTINUE

      DO 11 IA=1,NA  
	CLOSE(IA)
   11 CONTINUE
C
C
	IF(NFREQ.EQ.1)THEN
		TYPE = 'l'
	ELSEIF(NFREQ.EQ.2)THEN
		TYPE = 'c'
	ELSE
		TYPE = 'x'
	ENDIF
C
      OPEN(13,FILE=dscrp//'_20s.'//TYPE)
      OPEN(14,FILE=dscrp//'_30s.'//TYPE)
      OPEN(15,FILE=dscrp//'_40s.'//TYPE)
      OPEN(16,FILE=dscrp//'_50s.'//TYPE)
      OPEN(17,FILE=dscrp//'_60s.'//TYPE)
      OPEN(18,FILE=dscrp//'_70s.'//TYPE)
C
      DO 20 IA=1,NA
	II=12+IA
C
	DO 20 IP=1,4                       
	   WRITE(II,*) ' P= ',IP
	 DO 20 ID=1,NPTS
	   WRITE(II,*) TARGET(IA,1,IP,ID)
   20 CONTINUE
C
      DO 35 I=13,17
	CLOSE(I)
   35 CONTINUE
      STOP
      END
C
C After this job, you may need to calibrate the raw data.
C Those files are too long to mail, I'd better transfer by 'ftp' later.
C Good luck !
C (If you have questions about these files, don't hesitate to ask me.)