Holmberg I   [userid 1006 on aida32]


GENERAL COMMENTS
----------------

Holmberg I and II have been observed in the same scan in D and C array
configuration - that's why the splitting is done in this file.


SPLIT
=====

task 'SPLIT'; Sources 'UGC5139',''; Qual -1; Calcode ''; 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; Baddisk 0

Sources 'UGC 5139'  for  C, D array
Sources 'UGC5139'  for  B array


UVPLT: inspect _individual_ arrays
=====

tvinit
task 'UVPLT';  sources ''; flagver -1; dotv 1;
docal -1; doband -1; flagver -1; xinc 1; bparm 0;
bchan 60; echan 60;

   -> data look good



UVFLG - flag _all_ arrays
=====

flag data > 4 Jy. 

task 'uvflg'; antennas 0; timerang 0; aparm 0 4 0; stokes ''; bif 0;
eif 0; flagver 0; bchan 0; echan 0; opcode 'FLAG'; sources ''



DBCON
=====

task 'DBCON'; Reweight 0 0; Outname 'Dummy'; Dopos 0; Doarray 0; 
Fqtol 0

Outname 'HO_I_BCD'

UVLIN 
=====

task 'UVLIN'; Shift 0 0; Flux 0; Docont -1; Ichansel 15,65,1,0,100,115,1,0;
Order 1; Prtlev 0

AVSPC
=====

task 'AVSPC'; avoption ''; flagver -1;
bif 0;eif 0; channel 0;
Ichansel 16,65,1,0,100,115,1,0;
outclass 'UVCONT'


IMAGR  - 2 cubes with robust 5 and robust 0.5   + one continuum map
=====

noise = 0.93 mJy/beam  -  robust 5  [14.66" x 12.73"]
noise = 1.11 mJy/beam  -  robust 0  [6.38" x 5.32"]


task 'IMAGR'; cellsi 1.5; imsize 1024; uvwtfn 'un'; flux 0.0023; 
bchan 65; echan 100; robust 5; niter 100000; dotv -1; docalib -1; 
doband -1; flagver -1; outdi 5; outname ''; SUBARRAY 0; bmaj 0;
bmin 0; bpa 0;clr2n;  IMAGRPRM 0
outname 'HO_I_NA'


task 'IMAGR'; cellsi 1.5; imsize 1024; uvwtfn 'un'; flux 0.0023; 
bchan 65; echan 100; robust 0.5; niter 100000; dotv -1; docalib -1; 
doband -1; flagver -1; outdi 5; outname ''; SUBARRAY 0;
outname 'HO_I_RO'; IMAGRPRM 0


also do the same with residual scaling:

IMAGRPRM(5)=1
IMAGRPRM(6)=50
IMAGRPRM(7)=50

outname 'HO_I_NA_RES'
outname 'HO_I_RO_RES'


continuum:

massive sources ouside image -> much larger image

robust 5:

determine noise with imagr (niter 2000, flux 0):
noise = 0.17 mJy/beam

task 'IMAGR'; cellsi 1.5; imsize 2048; uvwtfn ''; flux 0.00042; 
niter 100000; dotv -1; outdi 5; calcode ''; docal -1; 
doband -1; clr2na; robust 5; outname 'HO_I_CONT_NA'

robust 0.5:

determine noise with imagr (niter 2000, flux 0):
noise = 0.18 mJy/beam

task 'IMAGR'; cellsi 1.5; imsize 2048; uvwtfn ''; flux 0.00047; 
niter 100000; dotv -1; outdi 5; calcode ''; docal -1; 
doband -1; clr2na; robust 0.5; outname 'HO_I_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

get rms in convolved cube: 1.7 mJy


### blank cube at 2 sigma

task 'BLANK'; opcode 'SELC'; dparm 1, 0, 10000, 0.0034, 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'

### create 10"x10" image

task 'convl'; bmaj 10; bmin 10; blc 0; trc 0; opcode ''; 
outclass 'CVL10'; doblank 0; factor 0; outse 0; outname 'HO_I_RO_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