NGC 1569 [userid 1039 on aida33] GENERAL COMMENTS ---------------- data are to be found in AIPS 1: 1 1276 NGC1569-B .CH 0 . 1 UV 22-JAN-2004 13:09:44 AIPS 1: 2 1276 NGC1569-B .LINE . 1 UV 22-JAN-2004 18:43:35 AIPS 1: 3 1276 NGC1569-C .CH 0 . 1 UV 22-JAN-2004 18:50:04 AIPS 1: 4 1276 NGC1569-C .LINE . 1 UV 22-JAN-2004 20:17:41 AIPS 1: 5 1276 NGC1569-D .CH 0 . 1 UV 22-JAN-2004 14:45:12 AIPS 1: 6 1276 NGC1569-D .LINE . 1 UV 22-JAN-2004 18:42:38 COMBINING THE DATA ------------------ SPLIT split data first (all 3 LINE files) NGC1569-B .LINE task 'SPLIT'; Sources 'NGC1569',''; Qual -1; Calcode '-CAL'; Timerang 0; Stokes ' '; Selband -1; Selfreq -1; Freqid 1; Bif 0; Eif 0; Bchan 1; Echan 0; Subarray 0; Docalib 2; Gainuse 2; Dopol -1; Blver -1; Flagver 1; Doband 1; Bpver 1; Smooth 0; Douvcomp 1; Aparm 0; Nchav 1; Chinc 1; Ichansel 0; NGC1569-C .LINE task 'SPLIT'; Sources 'NGC1569',''; Qual -1; Calcode '-CAL'; Timerang 0; Stokes ' '; Selband -1; Selfreq -1; Freqid 1; Bif 0; Eif 0; Bchan 1; Echan 0; Subarray 0; Docalib 2; Gainuse 2; Dopol -1; Blver -1; Flagver 1; Doband 1; Bpver 1; Smooth 0; Douvcomp 1; Aparm 0; Nchav 1; Chinc 1; Ichansel 0; NGC1569-D .LINE task 'SPLIT'; Sources 'NGC1569',''; Qual -1; Calcode '-CAL'; Timerang 0; Stokes ' '; Selband -1; Selfreq -1; Freqid 1; Bif 0; Eif 0; Bchan 1; Echan 0; Subarray 0; Docalib 2; Gainuse 2; Dopol -1; Blver -1; Flagver 1; Doband 1; Bpver 1; Smooth 0; Douvcomp 1; Aparm 0; Nchav 1; Chinc 1; Ichansel 0; ------------------------------------------------------- DBCON then combine 1st 2 split files task 'DBCON'; Reweight 0 0; Outname 'Dummy'; Dopos 0; Doarray 0; Fqtol -1; doarray 1 then dbcon the output with the remaining file and delete dummy file outna 'ngc1569-bcd' ---------------------------------------------- UVPLT: inspect DBCON file: ===== tvinit task 'UVPLT'; sources ''; flagver -1; dotv 1; docal -1; doband -1; flagver -1; xinc 1; bparm 0; bchan 50; echan 50; UVFLG - flag dbcon file > 5 Jy ===== task 'uvflg'; antennas 0; timerang 0; aparm 0 5 0; stokes ''; bif 0; eif 0; flagver 0; bchan 0; echan 0; opcode 'FLAG'; sources ''; infile '' ------------------------------------------------------------------------------- UVLIN task 'UVLIN'; Shift 0 0; Flux 0; Docont -1; Ichansel 12,26,1,0,100,112,1,0; Order 1; Prtlev 0 ------------------------------------------------------------------------------- AVSPC task 'AVSPC'; avoption ''; flagver -1; bif 0;eif 0; channel 0; Ichansel 12,26,1,0,100,112,1,0; outclass 'UVCONT' ------------------------------------------------------------------------------- IMAGR noise = 0.56 mJy/beam - robust 5 noise = 0.70 mJy/beam - robust 0 task 'IMAGR'; sources ''; flagver -1; cellsi 1.5; imsize 1024; uvwtfn ''; robust 5; niter 100000; flux =0.0014; dotv -1; bchan 26; echan 100; outname 'N1569_NA'; outdi 2; calcode ''; docal -1; doband -1;clr2na task 'IMAGR'; sources ''; flagver -1; cellsi 1.5; imsize 1024; uvwtfn ''; robust 0.5; niter 100000; flux =0.0017; dotv -1; bchan 26; echan 100; outname 'N1569_RO'; outdi 2; calcode ''; docal -1; doband -1;clr2na also do the same with residual scaling: IMAGRPRM(5)=1 IMAGRPRM(6)=50 IMAGRPRM(7)=50 outname 'N1569_NA_RES' outname 'N1569_RO_RES' continuum: robust 5: determine noise with imagr (niter 2000, flux 0): noise = 0.14 mJy/beam NOTE: strong sources require huge images! getn task 'IMAGR'; cellsi 1.5; imsize 2048; uvwtfn ''; flux 0.00035; niter 100000; dotv -1; outdi 2; calcode ''; docal -1; doband -1; clr2na; robust 5; outname '1569_CONT_NA' robust 0.5: noise = 0.15 mJy/beam getn task 'IMAGR'; cellsi 1.5; imsize 2048; uvwtfn ''; flux 0.00037; niter 100000; dotv -1; outdi 2; calcode ''; docal -1; doband -1; clr2na; robust 0.5; outname '1569_CONT_RO' ------------------------------------------------------- check that - all coordinates are 2000.0 - velocities are set to km/s subim first plane from beam cube - delete beam cube and move extracted plane to correct slot task 'subim'; blc 0 0 1; trc 0 0 1; opcode ''; xinc 1; yinc 1 outclass '' - delete beam cube - use renumber to move subim cube to old position - rename with outclass 'IBM001'; go renam --------------------------------------------------------- ### convolve NA cube: task 'convl'; blc 0; trc 0; outname ''; opcode ''; bmaj 30; bmin 30 outclass 'CVL30' get rms in convolved cube: 2.15e-3 ### blank cube at 2 sigma task 'BLANK'; opcode 'SELC'; dparm 1, 0, 10000, 0.0043, 0 outclass 'CVL_BL'; trc 0; blc 0; bchan 0; echan 0 ### blank cube by hand - output is the master cube task 'BLANK'; opcode 'TVCU'; doinvers -1; outclass 'master' txinc 2; tyinc 2 ### blank residuum-rescaled cubes with master cube: get2n master-cube task 'BLANK'; opco 'IN2C'; outclass 'MA_BLA'; dparm 0; blc 0; trc 0 ### transpose cubes task 'TRANS'; outclass ''; blc 0; trc 0; transcod '312' ### moment maps (take all values): task 'xmom'; flux -10000; icut -10000; blc 0; trc 0 ### pbcor both momnt 0 maps [first need to do 'altsw'] altsw task 'pbcor'; doinvers -1; coord 0; bparm 0; outclass 'X0_PBC' # replace blanks with 0 task 'remag'; pixval 0; blc 0; trc 0; outclass 'X0_P_R' # create 10"x10" image task 'convl'; bmaj 10; bmin 10; blc 0; trc 0; opcode ''; outclass 'CVL10'; doblank 0; factor 0; outse 0; outname 'N1569_NA_RES' ----------------------------------------------------------- ### write files out: RUN WRTPROCS task 'FITDISK'; doeot -1; DOSTOKES -1; DONEWTAB 1; format 0; blocking 10; catno 0 0; OUTFILE 'AIPS' for i=1 to 5; getn i; go fitdisk; end