; Use this main program as an example of how to call the get_scat_daily_v4 reading routine
; Reads version-3a data released October 2006
; Reads version-4  data released April 2011

; enter julian day and long year (day and year shown here match the verification file data)
for lyear=2000,2000 do begin		;qscat verify 2000,2000 : seawinds verify 2003,2003
  for jday=11,11 do begin			;qscat verify 11,11		: seawinds verify 115,115

		jul2day, jday,lyear, month,dayofmonth		; jul2day subroutine is from the Earth Space Research Group, UCSB, IDL library available online
		yy=string(lyear,format='(i4.4)')
		mm=string(month,format='(i2.2)')
		dd=string(dayofmonth,format='(i2.2)')

	; create file name of data file to read
		file_path='drive:\provide your own file path\'
		filename=file_path+yy+mm+dd+'.gz'  ;is the file *zipped*  ?
		print,filename

	; read scat data
		Get_scat_daily_v4,filename, mingmt,windspd,winddir,scatflag,radrain
		good=where(mingmt lt 1440, cntgood)
		if (cntgood eq 0) then goto, nofile

	; write out data in verification file

		for iasc=0,0 do begin  ;ascending or descending swath segments
		  for ilat=273,277 do begin		;for qscat and seawinds verify 273,277
			for ilon=169,174 do begin   ;for qscat and seawinds verify 169,174
				print,ilon+1,ilat+1,mingmt[ilon,ilat,iasc],$
				windspd[ilon,ilat,iasc],winddir[ilon,ilat,iasc],$
				scatflag[ilon,ilat,iasc],radrain[ilon,ilat,iasc],$
				format='(2i5,f6.0,2f8.2,i5,f8.2)'
			endfor  ;ilon
		  endfor  ;ilat
		endfor   ;iasc

		nofile: if(cntgood eq 0) then print,'file not found: ', filename

  endfor   ;jday
endfor   ;lyear

end