; example script to run coverage calculation code

scan_velocity = 120.
scan_length = 900.
offset = 0.
alternate_direction = 1
bolo_params_file = './bolo_params_feb04.txt'
platescale = 0
beam_fwhm = 30.
array_angle = 0 & scan_offset_angle = 10.7

scan_height = 900.
step_size = 26.
n_scans = ceil(scan_length / step_size)

ps_stub = ''
data_struct = ''

; radec scans, radec source coord, radec plotting system (source epoch)
scan_sys = 2
src_sys = 2
plot_sys = 2
plot_epoch_flag = 0

src_epoch = 2000.
src_fk4_flag = 0

; scan along RA
position_angle = 90.

; specify RA in degrees, then correct

src_x = [ 90, 0, 0]
src_y = [ 60, 0, 0]
src_coord = [src_x, src_y]
if (src_sys eq 2) then src_coord[0:2] = src_coord[0:2]/15.

; day and time are HST
day = [2004, 01, 01]
time = [00, 00, 00]
date = [day, time]

calc_track_cov, $
   scan_velocity, scan_length, n_scans, $
   step_size = step_size, offset = offset, $
   position_angle = position_angle, $
   alternate_direction = alternate_direction, $
   scan_sys = scan_sys, $
   src_coord = src_coord, $
   src_sys = src_sys, $
   src_epoch = src_epoch, $
   src_fk4_flag = src_fk4_flag, $
   plot_sys = plot_sys, $
   plot_epoch_flag = plot_epoch_flag, $
   date = date, $
   bolo_params_file = bolo_params_file, $
   platescale = platescale, beam_fwhm = beam_fwhm, $
   array_angle = array_angle, $
   scan_offset_angle = scan_offset_angle, $
   ps_stub = ps_stub, $
   data_struct = data_struct