;--------------------------------------------------------
;Title:		fy_scs_t.pro
;Purpose:	Main level script for running the SHEBA
;			cloud temperature statistics program for the
;			full year.
;
;Author:	Matthew Shupe
;Date:		2/10/00
;-------------------------------------------------------

;Set up common blocks
;-----------------------
common com_in,tim,jtim,lbase,ltop,ldepol,rbase,rtop
common com_raddn,raddn
common com_outt,rbaset,rtopt,rmidt,rmid,lbaset,ltopt,lmidt,lmid

;Create a vector of all Julian days for the full SHEBA year
; extending from 10/20/97 (2450742) to 10/2/98 (2451089).
;-----------------------------------------------------------
juliandays=lindgen(2451089L-2450742+1)+2450742
;juliandays=[2450875L,2450876L]
numdays=n_elements(juliandays)
restore,'radardown.dat'

;Now loop over the full year and .....
;-------------------------------------------
for i=0,numdays-1 do begin
	lai_scs,juliandays[i]
	stat_scs_t
	if i eq 0 then begin
		;The input data
		time=tim
		jtime=jtim
		baser=rbase
		topr=rtop
		basel=lbase
		topl=ltop
		depol=ldepol
		;The output data
		basert=rbaset
		toprt=rtopt
		midrt=rmidt
		midr=rmid
		baselt=lbaset
		toplt=ltopt
		midlt=lmidt
		midl=lmid
	endif else begin
		;The input data
		time=[tim,tim]
		jtime=[jtime,jtim]
		baser=[baser,rbase]
		topr=[topr,rtop]
		basel=[basel,lbase]
		topl=[topl,ltop]
		depol=[depol,ldepol]
		;The output data
		basert=[basert,rbaset]
		toprt=[toprt,rtopt]
		midrt=[midrt,rmidt]
		midr=[midr,rmid]
		baselt=[baselt,lbaset]
		toplt=[toplt,ltopt]
		midlt=[midlt,lmidt]
		midl=[midl,lmid]
	endelse
endfor

;Save the final years worth of data
;----------------------------------
save,basert,toprt,midrt,midr,baselt,toplt,midlt,midl,jtime,$
	baser,topr,basel,topl,depol,filename='scs_t.fy.dat'

end ;prog