pro get_stret_Dong2000_5min, stret1_fname, stret1_lwc5, $
							stret1_re5, stret1_error_frac

common five_min_data_input2, hrfrac5,jday5,height5,dbz5,nppts5,npts5,cbase5,ctop5, $
	vel5, width5,temp5, pressure5, n_layers5,bzl15,tzl15,btl15,ttl15,bzl25,tzl25,btl25,ttl25, $
	bzl35,tzl35,btl35,ttl35,bzl45,tzl45,btl45,ttl45,fluxclear5,solrat5,lwp5,vceil5,vap5

; puts the iwc and re onto the time-height grid.

for jjj=0,n_elements(stret1_fname)-1 do begin

cdfid=ncdf_open(strtrim(stret1_fname[jjj],2))	; open the netcdf file

time_did=ncdf_dimid(cdfid,'time')
height_did=ncdf_dimid(cdfid,'height')

ncdf_diminq, cdfid, time_did, char_strng, num_times
ncdf_diminq, cdfid, height_did, char_strng, num_heights

; get the id's of the variables to be read

num_sec_id=ncdf_varid(cdfid,'base_time')
time_offset_id=ncdf_varid(cdfid,'time_offset')
lwc_id=ncdf_varid(cdfid,'lwc')
height_id=ncdf_varid(cdfid,'height')
re_id=ncdf_varid(cdfid,'re')


; get the data



ncdf_varget, cdfid, num_sec_id, num_sec
ncdf_varget, cdfid, time_offset_id, time_offset
ncdf_varget, cdfid, height_id, hgt_stret

ijday_ihrfrac_fm_numsec, num_sec, time_offset, hrfrac, jday, yy, mm, dd, hh, mi, ss


ncdf_varget, cdfid, lwc_id, lwc_in
ncdf_varget, cdfid, re_id, re_in

;if max(radrat) gt 0.15 then lwc_in[where(radrat gt 0.1)]=-9999.

ncdf_close, cdfid


;ijday_ihrfrac_fm_numsec, num_sec, time_offset, hrfrac, jday, yy, mm, dd, hh, mi, ss


; loop through the observations, match times, then go through heights

for j=0,num_times-1 do begin
for k=0,num_heights-1 do begin

	;time_index=where( (abs(hrfrac[j]-hrfrac5)) le 0.17 )
	time_index=where( (abs(jday[j]-jday5)) le 0.0035 )

	z_index=where( (abs(hgt_stret[k]-height5)) le 200. )

	z_index=where( (abs(hgt_stret[k]-height5)) le 100. )

	if max(z_index) eq -1 and max(time_index) gt -1 then z_index=where( (abs(hgt_stret[k]-height5)) le 200. )


if max(time_index) gt -1 and max(z_index) gt -1 then begin

for jjj=0,n_elements(time_index)-1 do begin
for kkk=0,n_elements(z_index)-1 do begin


stret1_lwc5[z_index[kkk],time_index[jjj]]=lwc_in[k,j]
if re_in[k,j] gt 0. then stret1_re5[z_index[kkk],time_index[jjj]]=re_in[k,j]*1.e-6
if re_in[k,j] le 0. then stret1_re5[z_index[kkk],time_index[jjj]]=-9999.
stret1_error_frac[z_index[kkk],time_index[jjj]]=0.15

endfor	; for kkk=0,n_elements(z_index)-1 do begin
endfor	;for jjj=0,n_elements(time_index)-1

endif	;if max(time_index) gt -1 and max(z_index) gt -1 then begin

endfor ; do k=0,num_heights-1 do begin
endfor ;for j=0,n_elements(re_in)-1 do begin




endfor

return
end