#!/bin/csh -f # an example of script to produce MFS continuum image # from calibrated SMA data. # invert, clean, restor and cgdisp # jhz 2005-02-02 goto continue continue: #make dirty map \rm -r sgra-star.map sgra-star.beam sgra-star.icmp invert vis=sgra.selfcal map=sgra-star.map \ beam=sgra-star.beam \ imsize=512,512 cell=.2 sup=0 options=mfs,sdb,systemp \ line=channel,6,1,4 #clean map mfclean map=sgra-star.map beam=sgra-star.beam \ out=sgra-star.icmp gain=0.1 \ cutoff=0 niters=1500 region='boxes(200,200,312,312)' \rm -r sgra-star.icln #restore the image with a gaussian beam restor model=sgra-star.icmp beam=sgra-star.beam \ map=sgra-star.map out=sgra-star.icln #display the image cgdisp in=sgra-star.icln type=pixel region=quarter xybin=2 \ device=/xs nxy=1 labtyp=hms,dms \ options=full,beambr,wedge,fiddle,trlab \ range=0,0,lin,2 csize=1,1,0.8 #do statistics in an area off the source imstat in=sgra-star.icln region='boxes(10,10,30,30)' \ axes=ra,dec #do statistics in an area on the source imstat in=sgra-star.icln region=quarter \ axes=ra,dec #make a contore map and display on an x window device cgdisp in=sgra-star.icln type=c \ region=arcsec,'boxes(-15,-15,15,15)' \ xybin=1 device=/xs nxy=1 labtyp=arcsec,arcsec \ range=0,0,lin,2 options=full,beambr,wedge,trlab \ csize=1,1,0.8 slev=a,0.01 \ levs1=-5,3,4,6,9,13,18,24,31,41,51,71,91.100,150,300 #make a contore map and save in a postscript file. cgdisp in=sgra-star.icln type=c \ region=arcsec,'boxes(-15,-15,15,15)' \ xybin=1 device=sgra.ps/vps nxy=1 labtyp=hms,dms \ range=0,0,line,2 options=beambr \ csize=1,1,0.8 slev=a,0.01 \ levs1=-5,3,4,6,9,13,18,24,31,41,51,71,91.100,150,300 exit