pro get_scat_orbit_v04_test

data_root = 'R:\ftp\qscat\qscat_wind_vectors_v04\' ; edit to YOUR DATA PATH

data_folder = data_root + '04000to04999\'

file_name = 'winvec_04741_v04.gz'

file_name_and_path = data_folder + file_name

print, file_name

get_scat_orbit_v04,file_name_and_path,$
	atime,$
	phi_track,$
	lat,$
	lon,$
	iclass,$
	numamb,$
	select_ambiguity,$
	irain_scat,$
	wind_all,$
	dir_all,$
	sos_all,$
	wind,$
	dir,$
	wind_smooth,$
	dir_smooth,$
	wind_gcm,$
	dir_gcm,$
	rad_rain,$
	min_diff,$
	nudge_ncep,$
	nudge_ecmwf


for iscan=1171,1175 do begin
	for icel=17,22 do begin
		iamb=select_ambiguity[icel,iscan]
		iamb = iamb - 1; ; convert iamb from 1-based(1-4) to 0-based(0-3) index
		if (iclass[icel,iscan] ge 2) then begin   ; use only data for which iclass>=2
			print, iscan+1,icel+1, $
			lat[icel,iscan],lon[icel,iscan], $
			wind_all[iamb,icel,iscan], $
			dir_all[iamb,icel,iscan], $
			sos_all[iamb,icel,iscan], $
			wind_gcm[icel,iscan],dir_gcm[icel,iscan], $
			irain_scat[icel,iscan],rad_rain[icel,iscan], $
			format='(2i5,7f8.2,i4,f8.2,i3)'
		endif
	endfor
endfor

end