;					Wayne Pryor,  March 7, 1990
;					(303) 492-1259
;					edited Aug 17,90 WP/KES         <<
;			KES version   Aug 19, 1990			;<*
;			radical edits Aug 20, 1990			;/
;			added WIB and SIB	Aug 21, 90		;//
;			verified operation      Aug 22, 1990  KES
;  		Updated for general encounter mode  Nov 14, 90     	;<c
;	        Add a listing file for saved data file			;c1

;  This program creates EUV spectra by sector in an array called
;  DATA.  These spectra can conveniently be displayed by 
;    first running the program
;    then, using the statement:  SURFACE,DATA

;  Editorial comment: there may be loss of valid sector data due to JPL
;    (original) data record truncation/re-start.  KES
;    This might be recovered by eliminating the statements marked ***

;  Note the first file refers to 2.89 rpm spin, second,third to 3.15 rpm
;  The first file has 37 sectors, the second has 35 sectors

	FILENAM=' '
  	READ,'Supply file name (ex:EUV903171345_3181400)',FILENAM
  	CLOSE,1								;/
	EXT='.DAT' & filloc='torus1:[gll_raw.euv]'
	OPENR,1,filloc+FILENAM+EXT			;Open input data file
  	A=ASSOC(1,BYTARR(1200))			;Input data array access
;;	CLOSE,2 & OPENW,2,'TEST_INFO.LIS'	;Open a printf file
;	CLOSE,3 & OPENW,3,filenam+'.LIS'	;Open a printf file	;c1
;;	PRINTF,2,'Input file is',FILENAM	;record input file name
;;	PRINTF,3,'Input file is',FILENAM	;record input file name

; Four AACS records occur between each EUV data record; 300 F.P.values
  	AA=ASSOC(1,LONARR(300))   		;AACS input record access
  
	NBEG=1					;beginning record 
	READ,'Supply Spin Rate CHOICE (1=2.89 or 2=3.15)',CHOICE
	If choice eq 1 then SPINRATE=2.88295 else SPINRATE=3.1444
        IJK=0
        WHILE NOT EOF(1) DO BEGIN
          T=A(IJK) & IJK=IJK+1
        ENDWHILE
        RECMAX=IJK-1
;	READ,'Begin Rec, Final rec',NBEG,RECMAX
;	printf,3,'beg & end rec=',nbeg,recmax				;c1
	NUMSECT=(RECMAX-NBEG)*4			;guesstimate total # sectors
	READ,' Enter cmded number of sectors (E1=125,Ven=120)',CMD_SECS
;	printf,3,'cmded sectors=',cmd_secs				;c1

; Initilize arrays, params
  	DATA=INTARR(129,128)	    ;Holds valid data (1-127,0-127)(sect,wav)
  	DATA_TEMP=INTARR(129,128)	;Temp until validated	"  "     ;/
        data2=fltarr(129,128)
        oplus=fltarr(126)
        bkg=fltarr(126)
        BOO=flTARR(126)
        he=fltarr(126)
        hb=fltarr(126)
        ha=fltarr(126)
  	NRIMS=0				;Counts number of records (1-i)(@RIM)
  	NSPINS=0			;Cnts # Spin packets=#revs (1-i)
  	ISB=-1				;Cnts Sector boundaries found (0-i)
  	SYNC=0				;=1 when spin packet synced (1st FE)
  	PHOTON_CNT=FLTARR(NUMSECT)-1	;Cnt total # photons in one RIM
	INIT_CMD=-1			;Have seen sector cmd yet?Then hold it
	SECS_SEEN=INTARR(NUMSECT)-1	;Record # sectors seen in each spin
	WIB=FLTARR(129)-1		;IW and IS count 1-128		;//
	SIB=FLTARR(129)-1		;            "     "		;//
        sibprime=sib
	IS=-1 & IW=-1			;Both Addrs unknown;don't use til know
	MOD91=FLTARR(NUMSECT)-1		;MOD91 Time of WIB values
	RTI=FLTARR(NUMSECT)-1		;RTI    "   "   "    "
	HOLD_IS=FLTARR(NUMSECT)-1	;Hold the consecutive IS values
	HOLD_SIB=FLTARR(NUMSECT)-1	;Hold the consecutive SIB values
	HOLD_IW=FLTARR(NUMSECT)-1	;Hold the consecutive IW values
	HOLD_WIB=FLTARR(NUMSECT)-1	;Hold the consecutive WIB values
	M_IS=0 & M_IW=0			;Indicies for above hold arrays
	LAST_IS=-1			;The last IS value
	LAST_IW=-1			;The last IW value

; Outer, file, loop: read each rec
	J=NBEG					;1ST rec to use
MORE:   INARR= A(J)				;Get an input rec
	TIM=FIX(INARR)				;convert Hdr info to Int.vals
    
; Use only recs which have RIM start flags at synced location (rest are
; AACS or bad/truncated/offset records			;RIM start
	IF (INARR(80) EQ '7E'XB) AND (INARR(81) EQ '7E'XB) THEN GOTO,SCIREC
	GOTO,FINIS

SCIREC:	NRIMS=NRIMS+1				;Increm RIM cntr (=Hkp cntr)
        PRINT,'record #, tim',J,TIM(0),TIM(2),TIM(4),TIM(6),TIM(8),$
		FIX(INARR,23)+INARR(26)
;;        PRINTF,2,'record #,tim',J,TIM(0),TIM(2),TIM(4),TIM(6),TIM(8),$
;;		FIX(INARR,23)+INARR(26)
;          PRINTF,3,'record #,tim',J,TIM(0),TIM(2),TIM(4),TIM(6),TIM(8),$
;		FIX(INARR,23)+INARR(26)
    	SEC_SEEN=INARR(87) AND '7F'XB 		;Mask Cmded # SECTORS
	IF CMD_SECS NE SEC_SEEN THEN BEGIN
		PRINT,'Sectors seen .ne. cmd',SEC_SEEN
		if sec_seen eq 0 then goto,keepit			;c1
	;	PRINTf,3,'Sectors seen .ne. cmd',SEC_SEEN		;c1
		PRINT,'Choose what to do:'
		PRINT,'		0 = Accept as a data hit'
		Print,'		1 = Configuration changed, finish up'
		Print,'		2 = Just STOP'
		READ,WHAT_NOW
		IF WHAT_NOW EQ 2 THEN STOP
		IF WHAT_NOW EQ 1 THEN GOTO,FINIS
		END
keepit:	CMD_SCAN_SECS=INARR(86) AND '7F'XB	;Mask Cmded # scan per sec
;;;	PRINT,'  Cmded sectors=',CMD_SECS
;;	PRINTF,2,'  Cmded sectors=',CMD_SECS
;;	PRINTF,3,'  Cmded sectors=',CMD_SECS
	IF INIT_CMD LT 0 THEN INIT_CMD=CMD_SECS	;Keep copy orig/last cmd
	IF CMD_SECS NE INIT_CMD THEN BEGIN	;Warn of a changed cmd
;;;	  PRINT,'WARNING: CMDED # SECTORS CHANGED;Orig,new:',INIT_CMD,CMD_SECS
;;	  PRINTF,2,'WARNING:CMDD # SECTORS CHANGD;Orig,new:',INIT_CMD,CMD_SECS
	  ENDIF

; Handle RIM packet info at start of each data rec			;//
	LAST_IS=IS				;Keep tract of prev.addr valu
	IS=INARR(88)				;SIB addr
	SS=INARR(89)*256.+INARR(90)		;SIB value
	IF IS NE LAST_IS+1 THEN $
             PRINT,'                         Last IS, this IS',LAST_IS,IS
	SIB(IS)=SS				;store it
        sibprime(is)=total(data(is,0:125))
;        for i2=0,125 do begin
;           data2(is,i2)=data(is,i2)*(sib(is)/sibprime(is))/(is*60.66)   ;counts per second if spectrum preserved
;        endfor
	HOLD_IS(M_IS)=IS & HOLD_SIB(M_IS)=SS & M_IS=M_IS+1	;Store consec
	LAST_IW=IW			 	;Keep tract of prev.addr valu
	IW=INARR(91)			     ;Addr of WIB in Next spin packet
;;	PRINTF,3,'     IS,SIB,IW=',IS,SS,IW
;;	PRINTf,2,'Comparision of RIM values:Hdr,internal=',INARR(26),INARR(85)
						;Inst.mode is in INARR(82)

; Now handle data part of rec			;Old FOR L=92,1171 DO BEGIN

     	L=92				;data bytes (after FE/FE packet) ;//
LOOP: 	IF L GT 1171 THEN GOTO,FINIS		;test for last byte
      	IVALU=INARR(L)				;This value
	IF IVALU EQ 'FF'XB THEN GOTO,NEXT	;Filler valu, ignore
 SPIN:  IF IVALU EQ '0FE'XB THEN BEGIN		;Spin packet marker?
;;	  PRINTF,2,'  ISB,FE packet:',ISB,INARR(L:L+10)
	  WW=INARR(L+2)*256.+INARR(L+1)		;calc WIB value		  ;//
	  IF IW GT -1 THEN BEGIN		;know valid addr so store WIB
	    IF IW GT 128 THEN IW=IW-128
	    RRTI=INARR(L+4)			;  "	"	"	  ;//*
	    IF IW EQ LAST_IW+1 THEN GOTO,KEEP   ;It's OK                   //*
	    IF IW NE LAST_IW+1 THEN BEGIN       ;See if off by only one    //*
             PRINT,'                         Last IW, this IW,RTI',$
                                                   LAST_IW,IW,RRTI
;;             PRINTF,3,'                         Last IW, this IW,RTI',$
;;                                                   LAST_IW,IW,RRTI
		ENDIF
	    IF (IW EQ LAST_IW+2) AND ((RRTI EQ 90) OR (RRTI EQ 0)) THEN BEGIN  ;//*
;;;	      PRINT,'                        IW ACCEPTED AS IW-1'	  ;//*
;;	      PRINTF,2,'                        IW ACCEPTED AS IW-1'	  ;//*
;;	      PRINTF,3,'                        IW ACCEPTED AS IW-1'	  ;//*
	      IW=IW-1 & ENDIF						  ;//*
KEEP:	    WIB(IW)=WW 
;;	    PRINTF,3,'      IW,WIB=',IW,WW
	    HOLD_IW(M_IW)=IW & HOLD_WIB(M_IW)=WW & M_IW=M_IW+1
	    LAST_IW=IW & IW=IW+1	 	;Keep tract of prev.addr valu
	    MOD91(NSPINS)=INARR(L+3)		;store time of WIB value  ;//
	    RTI(NSPINS)=INARR(L+4)		;  "	"	"	  ;//
	    ENDIF							  ;//
;  When reach new Spin Packet Start marker, store any 'valid' data where
;  valid means if the number of sectors seen agrees with # commanded
;   Recall that the cmded # sectors equels sector zero (=one 21.4 ms inst
;   period used to zero the spectrom registrs) plus [n-1] actual data sectors.
;  Temporarily, take everything after the cmded # is found but don't take
;   less than the cmded # since spin rate caused roll over to next revolution
;   may be in effect

	  IF ISB GE CMD_SECS-1 THEN BEGIN       ;found = expected-1?
;;;	    PRINT,'  at FE, ISB=',ISB
;;	    PRINTF,2,'  at FE, ISB=',ISB
	    DATA=DATA+DATA_TEMP(0:CMD_SECS,0:127)	;Store valid data
	    DATA_TEMP=INTARR(129,128)		;re-zero temp array  	   ;/
;;.............................^ check this again
	    SECS_SEEN(NSPINS)=ISB
;;;	    ISB=0				;sect zero is valid so use -1
	    ISB=-1		;;;		;Zero # Spin Packets read  ;/
	    NSPINS=NSPINS+1			;increment # spin s
	    IF NSPINS GT NUMSECT THEN BEGIN	;too many for array size?
	      PRINT,'TOO MANY SPIN PACKETS for arrays'
	      STOP & ENDIF
	    ENDIF
	  L=L+5		    ;Skip FE/WIB/MOD91/RTI bytes of Spin Packet marker
	  SYNC=1				;Set synced=yes
	  GOTO,LOOP							;<*

          ENDIF ELSE BEGIN
 ;if not a spin packet marker then either data or (sector bit+filler)

	  IF SYNC EQ 0 THEN GOTO,NEXT		;Synced yet?If not,thro away
;;;	  IF ISB LE CMD_SECS THEN BEGIN		;;;have seen < #expected?
                                                ;again, take 'em all
	    IF (IVALU AND '80'XB) EQ 0 THEN BEGIN	;Sector boundary?
	      ISB=ISB+1				;Inc sector cnt
;	      PRINT,'  L,value,ISB',L,IVALU,ISB 
;	      PRINTF,2,'  L,value,ISB',L,IVALU,ISB 
	      ENDIF
	    IVALU=IVALU AND '7F'XB		;Mask data bits
; Pixels 126 and 127 are the "hot" pixels; 127 also represents t/m hits and
; (sector bit+filler) values, which indicates how much T/M rate is not used
;Old	    	                                ; IF IVALU LT 126 THEN BEGIN
;New						;take 'em all
	    IF (ISB LT 0) OR (ISB GT 128) THEN BEGIN	;ISB out of range?
;;;	      PRINT,'  ISB error: out of array range:',ISB
	      GOTO, NEXT 
              ENDIF
	    DATA_TEMP(ISB,IVALU)=DATA_TEMP(ISB,IVALU)+1 ;Inc photon cnt
	    PHOTON_CNT(NRIMS)=PHOTON_CNT(NRIMS)+1	;Cnt photons/rec

          ENDELSE	;end of IVALU testing choices
NEXT:	L=L+1 
  	GOTO,LOOP

FINIS:  J=J+1					;finished reading records?
	IF J LE RECMAX THEN GOTO,MORE

	PRINT,'final record #',J
;	PRINTf,3,'final record #',J
;	PRINTf,3,' NRIMS=',NRIMS
;	PRINTf,3,' Cmded # Sectors',CMD_SECS
;	PRINTf,3,' Total # SPIN Boundaries',NSPINS
;	PRINTf,3,' # Sectors seen in each spin packet',SECS_SEEN(0:NSPINS)
;	PRINTf,3,' Photons/RIM',PHOTON_CNT(0:NRIMS)
					;note that Sec=0 should be empty:
	goto,calc				;;;
	PRINT,' Matrix of good data:',DATA(0:CMD_SECS,0:128)  
;	PRINT,' Matrix of remaining data:',DATA_TEMP(0:CMD_SECS,0:128)
	PRINT,' SIB (valid=1-128):',SIB
	PRINT,' WIB (valid=1-128):',WIB
	PRINT,' Number of IS & SIB values stored=',M_IS
	PRINT,' Stored SIB ADDR values, (in order):',HOLD_IS(0:M_IS)
	PRINT,' Stored SIB, (in order):',HOLD_SIB(0:M_IS)
	PRINT,' Number of IW & WIB values stored=',M_IW
	PRINT,' Stored WIB ADDR values, (in order):',HOLD_IW(0:M_IS)
	PRINT,' Stored WIBs, (in order):',HOLD_WIB(0:M_IW)
	

;Calc total time and sector size
;Data=Spin rate=3.15 
;time in secs=revs*(secs/minute)*(minutes/rev)

calc:	NUM_SEC_PER_SCAN=21.3E-3
            for i=0,125 do begin
               oplus(I)=total(data2(i,44:54))
               he(I)=total(data2(i,1:11))
               hb(I)=total(data2(i,75:85))
               BOO(I)=TOTAL(DATA2(I,87:97))
               ha(I)=total(data2(i,108:118))
               bkg(I)=total(data2(i,20:30))
            endfor
          he2=fltarr(nrims+10) & bkg2=he2 & oplus2=he2 & ha2=he2 & hb2=he2
          BOO2=HE2
;          for i=0,nrims+9 do begin
;            he2(i)=total(he(*,i))
;            bkg2(i)=total(bkg(*,i))
;            oplus2(i)=total(oplus(*,i))
;            ha2(i)=total(ha(*,i))
;            hb2(i)=total(hb(*,i))
;            BOO2(I)=TOTAL(BOO(*,I))
;          endfor
          SPECTRA=FLTARR(128)
          FOR I=0,127 DO SPECTRA(I)=TOTAL(DATA(*,I))
    	  T=NRIMS*(60./SPINRATE)
    	  PRINT,T, 'TOTAL SECONDS INCLUDING GAP ON EACH SPIN'
;    	  PRINTF,3,T, 'TOTAL SECONDS INCLUDING GAP ON EACH SPIN'
;sectorsize=deg/DATA=(rev/min)*(deg/rev)*(min/sec)*(sec/sector)
    	  SECTORSIZE=SPINRATE*(360./60.)*(CMD_SCAN_SECS*NUM_SEC_PER_SCAN) 
					;IN DEGREES
    	  PRINT,T*(CMD_SECS*SECTORSIZE)/360.,'TOTAL SECONDS OF ACTUAL DATA'
;    	  PRINTF,3,T*(CMD_SECS*SECTORSIZE)/360.,'TOTAL SECONDS OF ACTUAL DATA'
  
	CLOSE,2 & CLOSE,3
	PRINT,' For Surface plot of DATA, use .CONT'
	STOP
;	SURFACE,DATA(0:cmd_secs-1,0:125),55,30,.4   ;Plot only expected area
	print,filloc+FILENAM+EXT		
	read,'copy screen; ready',iry
	plot,sib
	print,'SIB for ',filloc+FILENAM+EXT		
	read,'copy screen; ready',iry
	plot,wib
	print,'WIB for ',filloc+FILENAM+EXT		
	read,'copy screen; ready',iry
;	close,3 & openw,3,filenam+'.wib',129*4
;	new=assoc(3,fltarr(129))
;	new(0)=wib
;	new(1)=sib
	close,3 & openw,3,filenam+'.surf',128*4
	new=assoc(3,fltarr(CMD_SECS,128))
	new(0)=data
	close,3
	END