pro read_CPL_CRS_boundaries

gbin = 22 ;ground return is in bin 22.
numrecs = 2000
xaxt = strarr(numrecs)
nlay = lonarr(numrecs)
gnd_ht = fltarr(numrecs)
lay_bot = fltarr(numrecs,5)  ;CPL finds up to 5 layers
lay_top = fltarr(numrecs,5)
radnum_lay = intarr(numrecs)
radlay_top = fltarr(numrecs,10)  ;CRS finds up to 10 layers.
radlay_bot = fltarr(numrecs,10)

x1 = ' '
x2 = 0L
x3 = 0.0
x4 = 0.0
x5 = 0.0
x6 = 0
openr,1,'Macintosh HD:CPL NRB programs:CRS_data_jul29:layer_boundaries_jul29.dat'

r = 0
while (not eof(1)) do begin
	readf,1,x1
	xaxt(r) = x1    ;time
	readf,1,x2,x3
	nlay(r) = x2    ;number of layers found in CPL data.
	gnd_ht(r) = x3  ;in km.  if = -999, then no ground return was seen for the lidar.
	for i = 0,nlay(r)-1 do begin
		readf,1,x4,x5
		lay_top(r,i) = x4  ;for lidar, in km.
		lay_bot(r,i) = x5  ;for lidar, in km.
		endfor
	readf,1,x6
	radnum_lay(r) = x6  ;number of layers found in CRS data.
	if (radnum_lay(r) NE 0) then begin
		for i = 0,radnum_lay(r)-1 do begin
			readf,1,x4,x5
			radlay_top(r,i) = x4  ;for radar, in km.
			radlay_bot(r,i) = x5  ;for radar, in km.
			endfor
		endif		
			
	r = r + 1
	endwhile
				
			
stop

return
end