; Use this main program as an example of how to call the get_ascat_daily reading routine

version = '_v02.1'

; edit below line to YOUR DATA PATH
file_path = 'G:\ftp\ascat\metopa\bmaps' + version + '\' 

; enter year, month, day (dates shown here match the verification file data)

for year=2007,2007 do begin
  for month=10,10 do begin
    for day=22,22 do begin

		yyyy = string(year,format='(i4.4)')
		mm =   string(month,format='(i2.2)')
		dd =   string(day,format='(i2.2)')


	; create file name of data file to read
		filename = file_path + 'y'+yyyy+'\' + 'm'+mm+'\' + 'ascat_' + yyyy+mm+dd + version + '.gz'
		print,filename

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

	; write out data in verification file

		for iasc=1,1 do begin  ;ascending or descending swath segments
		  for ilat=273,277 do begin
			for ilon=169,174 do begin
				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],$
				sosmap[ilon,ilat,iasc],$
				format='(2i5,f6.0,2f8.2,i5,2f8.2)'
			endfor  ;ilon
		  endfor  ;ilat
		endfor  ;iasc

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

    endfor   ;day
  endfor   ;month
endfor   ;year

end