! 
! R. Lucas 19-mar-1997 
! Control panel for self calibration of UV data. 
! 
!---------------------------------------------------------------------- 
! 
if exist(self_integ) then 
  del /var self_ampli self_phase self_integ cal_ampli cal_phase cal_gain 
endif 
def integer self_channel /global 
def real self_integ /global 
def logical self_ampli self_phase cal_ampli cal_phase /global 
def cha cal_gain*80 /global 
let self_integ 60. 
let self_phase yes 
let cal_phase yes 
let self_channel 1 
! 
gui\panel "Self calibration control panel" pr:graphic/x_self.hlp 
sic\let name 'name' /prompt "Generic name " 
sic\let type 'type' /prompt "Map file type" 
! 
gui\button "@ p_uvall" UVALL "UV Plots" gildas_help:uv_fit.hlp "UVALL parameters" 
sic\let xtype 'xtype' /prompt "X data" /choice u v angle radius time date scan number amp phase real imag weight 
sic\let ytype 'ytype' /prompt "Y data" /choice u v angle radius time date scan number amp phase real imag weight 
let first 'first' /prompt "First channel" 
let last 'last' /prompt "Last channel" 
let limits 'limits' /prompt "Plot limits" 
let error_bars 'error_bars' /prompt "Plot error bars" 
! 
! 
gui\button "@ p_uv_map" UV_MAP "Mapping from UV data" gildas_help:uv_map.hlp "UV_MAP parameters" 
let uv_taper 'uv_taper[1]' 'uv_taper[2]' 'uv_taper[3]' /prompt "UV Taper" 
let weight_mode 'weight_mode' /prompt "Weighting mode" /choice natural uniform 
let map_size 'map_size[1]' 'map_size[2]' /prompt "Map size (in pixels)" 
let map_cell 'map_cell[1]' 'map_cell[2]' /prompt "Pixel size (in arcsec)" 
let uv_shift 'uv_shift' /prompt "Shift and rotate map on specified center ?" 
let map_ra 'map_ra' /prompt "Right Ascension" 
let map_dec 'map_dec' /prompt "Declination" 
let map_angle 'map_angle' /prompt "Angle from North" 
! 
let uv_cell 'uv_cell[1]' 'uv_cell[2]' /prompt "UV Cell size and min. weight" 
let mcol 'mcol[1]' 'mcol[2]' /prompt "First and last channel" 
let wcol 'wcol' /prompt "Weight channel" 
! 
let convolution 'convolution' /prompt "Convolution function" /index box sinc expo expo_sinc spheroidal 
! 
gui\button "@ p_clean" CLEAN "Cleaning dirty image" gildas_help:clean.hlp "CLEAN input parameters" 
sic\let gain 'gain' /prompt "Loop gain" /range 0 1 
sic\let niter 'niter' /prompt "Number of components" 
sic\let fres 'fres' /prompt "Maximum residual (fraction of peak value)" 
sic\let ares 'ares' /prompt "Maximum residual (absolute value)" 
sic\let blc 'blc[1]' 'blc[2]' /prompt "Bottom left corner of clean box" 
sic\let trc 'trc[1]' 'trc[2]' /prompt "Top right corner of clean box" 
sic\let major 'major' /prompt "Major axis" 
sic\let minor 'minor' /prompt "Minor axis" 
sic\let angle 'angle' /prompt "Axis position angle" 
! 
sic\let method 'method' /choice clark simple mrc /prompt "Cleaning algorithm" 
sic\let keep 'keep' /prompt "Keep cleaning after convergence ?" 
! 
gui\button "vec\run uv_cct pr:graphic/uv_cct.init /now" UV_CCT "Model from clean comp." gildas_help:uv_cct.hlp "No parameters ..." 
! 
gui\button "vec\run uv_gain pr:graphic/uv_gain.init /now" UV_GAIN "Self calibrate" gildas_help:uv_gain.hlp "UV_GAIN input parameters" 
sic\let self_integ 'self_integ' /prompt "Integration time" 
sic\let self_ampli 'self_ampli' /prompt "Do Amplitudes ?" 
sic\let self_phase 'self_phase' /prompt "Do Phases ?" 
sic\let self_channel 'self_channel' /prompt "Channel used ?" 
! 
gui\button "vec\run uv_cal pr:graphic/uv_cal.init /now" UV_CAL "Calibrate" gildas_help:uv_cal.hlp "UV_CAL input parameters" 
sic\let cal_gain 'cal_gain' /prompt "Gain table ?" /file "*-gain.uvt" 
sic\let cal_ampli 'cal_ampli' /prompt "Do Amplitudes ?" 
sic\let cal_phase 'cal_phase' /prompt "Do Phases ?" 
! 
gui\button "@ x_lmv.graphic" DISPLAY "Displaying results" x_whole.greg "DISPLAY Input parameters" 
sic\let first 'first' /prompt "First channel" 
sic\let last 'last' /prompt "Last channel" 
! 
sic\let size 'size' /prompt "Size of plotted area (sec)" 
sic\let spacing 'spacing' /prompt "Contour step" 
sic\let cross 'cross' /prompt "Center cross size (sec)" 
sic\let mark 'mark' /prompt "Box marking type" /choice velocity frequency channel none 
sic\let beam_type 'beam_type' /prompt "Additional beam plotted" /choice dirty clean none 
sic\let extra 'extra' /prompt "Extra data plotted" /choice none coverage 
! 
sic\let do_nice 'do_nice' /prompt "Plot clean beam ?" 
sic\let do_bit 'do_bit' /prompt "Show bitmap images ?" 
sic\let do_contour 'do_contour' /prompt "Draw contours ?" 
sic\let do_grey 'do_grey' /prompt "Fill contours ?" 
sic\let do_tree 'do_tree' /prompt "Structure plot ?" 
! 
on error return 
gui\go