WRITE(IUN_RAD) (HEAD(I),I=1,512)  ! header info            BYTE

	NTL		! NTL: number of fractional series         INTEGER*2
	KSIG = -4       ! new data format, 4 antenna square array  INTEGER*2
	KWIN = -2       ! new data format, subcode -2              INTEGER*2
	NRUN = 1	! first range cell in data set             INTEGER*2
	NROB = NRAN	! last range cell in data set              INTEGER*2
	NRST = 1	! range cell step width                    INTEGER*2

	WRITE(IUN_RAD) NTL,KSIG,KWIN,NRUN,NROB,NRST

	DO IRAN = NRUN,NROB,NRST                      ! loop through range cells
	                                              ! write signal and noise
	 WRITE(IUN_RAD) (SC(IANT),SN(IANT),IANT=1,4)  ! REAL*4 SC(4), SN(4)
	 NUM_ANG = 0                                  ! INTEGER*4 NUM_ANG
	 DO     IANG = 1,360
	  IF(KK(IANG).GT.0) THEN 
	   NUM_ANG = NUM_ANG + 1                      ! REAL*4 V_RADIAL(6,360)
	   V_RADIAL(1,NUM_ANG) = FLOAT(IANG)          ! angle (1...360)
	   V_RADIAL(2,NUM_ANG) = FLOAT(KK(IANG))      ! # of contributions
	   V_RADIAL(3,NUM_ANG) = URM(IANG)            ! mean current velocity
	   V_RADIAL(4,NUM_ANG) = URS(IANG)            ! variance
	   V_RADIAL(5,NUM_ANG) = USM(IANG)            ! sum of weights (S/N)
	   V_RADIAL(6,NUM_ANG) = PWR(IANG)            ! sum of signal power
          END IF
	 END DO 
	 WRITE(IUN_RAD) NUM_ANG                       ! write number of angles
	 IF(NUM_ANG.GT.0) THEN                        ! write existing data
	   WRITE(IUN_RAD) ((V_RADIAL(I,IANG), I=1,6), IANG=1,NUM_ANG)
	 END IF
	END DO