! 
! DO MAP IN MOSAIC CASE : 
!       FOR I 1 TO NF 
!	    UV_MAP 
!       NEXT 
!       MAKE_MOSAIC 
! 
if (.not.exist(nfield)) then 
  define integer nfield /global 
  define double bmin hppb /global 
endif 
if (.not.exist(show_result)) then 
  define logical show_result /global 
  let show_result .true. 
endif 
let uv_shift .true. 
let nfield 2 
let bmin 0.2 
let hppb 50 
! 
gui\panel "UV_MOSAIC control panel" pr:map/mapping_mosaic.hlp 
! 
let name 'name' /prompt "Generic name " 
let nfield 'nfield' /prompt "Number of fields" 
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 wcol 'wcol' /prompt "Weight channel" 
let map_ra 'map_ra' /prompt "RA center of mosaic" 
let map_dec 'map_dec' /prompt "Dec center of mosaic" 
let hppb 'hppb' /prompt "Half Power Primary Beam (sec)" 
let bmin 'bmin' /prompt "Truncation level" 
let show_result 'show_result' /prompt "Show resulting dirty image?" 
! 
! UV_MAP 
! 
gui\button "@ do_mosaic.map" "DO MOSAIC" "Do mosaic from UV data" gildas_run: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 "RA center of mosaic" 
let map_dec 'map_dec' /prompt "Dec center of mosaic" 
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 
! 
on error return 
gui\go "@ do_mosaic.map"