@gia_hdfread.pro
@gia_findwats.pro
@gia_asciiread.pro
@gia_duct.pro
@gia_idm.pro
@gia_wats.pro
@gia_fileplot.pro
@gia_checkit.pro


;|-----------------------------------------------------------------|
;|                                                                 |
;|   Begin Main Program                                            |
;|                                                                 |
;|-----------------------------------------------------------------|

@gia_var.inc
maxrank=long(2)
ddimsizes=lonarr(maxrank)
idimsizes=lonarr(maxrank)

cyclefile=' '
number=' '

set_plot,'Z'

!p.thick = 2
!x.thick = 2
!y.thick = 2

Print,'Enter Directory Number(01 - 14)'
read, number
cyclefile=strtrim(number,2)+'.txt'
openr,1,cyclefile
loop:error_status=0
;---------> Begins Main Loop
while not eof(1) do begin

	hdfread   ;-------> Call Subroutine HDFREAD
	findwats  ;-------> Call Subroutine FINDWATS

	catch, error_status
        if error_status eq -250 then goto,loop

	filename='/usr33/users/klenzing/data/wats/'+strtrim(year,2)+strtrim(date,2)+'_de2_wats_2s_v01.asc'
	print,filename
	openr, 2, filename

	asciiread ;-------> Call Subroutine ASCIIREAD

	CLOSE,2
	FREE_lun,2

	duct      ;-------> Call Subroutine DUCT
	idm       ;-------> Call Subroutine IDM

	if jcnt eq 0 then goto,loop

	wats      ;-------> Call Subroutine WATS

	checkit   ;-------> Call Subroutine CHECKIT

endwhile
;---------> Ends Main Loop

close,1
free_lun,1

Device, Z_Buffer=1
set_plot,'x'

!p.thick=0
!x.thick=0
!y.thick=0


end;----------------------------------------------------------