2DARK天文志

1.设置数据路径，启动aips软件：
cd $AIPS_ROOT/DATA/XCOLD97
export MYDATA=`pwd`
cd ~
aips tv=local
900                                           (自己的id)

2.装载数据：
task 'fillm';default;
datain 'MYDATA:AH330_';vlaobs 'AH330'
doweight 1;band ''
nfiles 0;
ncount 10;inp
手动输入go，pcat：先重新查看参数设置，以免出错。
（附注：1988年一月一日以前磁带上的数据在用fillm来读时需要做一个转换，不过如果是下载的AIPS friendly数据的话就不存在这个问题）

3.查看各个表：
task 'prtab';
查看天线阵列:task 'prtan';default;getn 2;prtlev 0;inp;go选择邻近天线作为参考天线refant 23：
imh查看观测频率

4.标记坏点，生成FG表:
task 'quack';default;getn 2
flagv 1
opcode ''
aparm 0
sour '';inp
go

(1)校正antenna based errors:
task 'snplt';default;getn 2;
inext 'ty';timer 0;
optype 'tsys'
sour ''
nplots 16
dotv 1
inp
go                     找到有问题的天线，主要看系统噪声是否有大幅增长
task 'uvflg';default;getn 2;timer 0;antenna 19 0;单独列出异常天线
stokes ''
flagv 1
inp
go
antenna 26 0;timer 0 6 10 0 0 6 20 0;inp;具体时间范围
go

5.查看观测源:
所有源：
task 'listr';default;getn 2
optype 'scan'
docrt 132;inp
go
校正源：calc '*';go
打印源信息：
docrt -1;calc '';OUTPRINT 'MYDATA:iterms/AH330_sources';go
calc '*';OUTPRINT 'MYDATA:calcsources';go;docrt 132;OUTPRINT ''
查看基线长度，计算分辨率，视场，继而计算cellsize和imsize：
task 'uvplt';default;getn 2
sour '3C286' '';dotv 1;初级校正源
bparm 3 2
inp
go
sour '2229+695' '';go     次级校正源
sour 'CEPA' '';go         目标源
查找替换：1.项目名称；2.参考天线；3.初级校正源；4.次级校正源；5.目标源；6.文件名。
spatial resolution=lamda/distance=1/600k~0.34"
cellsize取其三分之一约0.15
imsize=[3E8/(4.8851G*25*cellsize)]*2～2546取2048

(2)校正baseline based errors:
task 'uvplt';default;getn 2;sour '2229+695' ''
bparm 6 7 2 0;docalib -1;dotv 1;go
bparm 0;go
bparm 0 2 0;go
task 'vplot';getn 2;sour '3C286' '';dotv 1;bparm 0;nplots 5;inp;
flag -1;go
flag 1;go                查看quack是否充分
task 'quack';sour '';flagv 1;opcode 'beg';aparm 0 10/60 0剔除前10秒
inp;
go
task 'uvplt';sour '2229+695' ''；bparm 0;dotv 1;flag 1;go
task 'vplot';sour '3C286' '';dotv 1;bparm 0;nplots 5;flag 1;go再次查看
task 'tvflg';stokes 'i';docalib -1;flagv 1;inp;go
AMP V DIFF, LOAD大致看一下(白色是坏点)
——实际，当观测数据足够多时（或者仪器精度足够高，天线数目足够多时），FLAG就不是很必要了，而要把重心放在self-calibration上。

6.UVCOP

7.设置初级校正源的流量值大小:
task 'setjy';default;getn 2
source '3C286',''
optyp 'calc';freqid 1;go
或者手动设置（尤其是当找不到合适的flux矫正源时或者未观测flux矫正源时）：
task 'setjy';default;getn 2;source '3C286','';ZEROSP 10.00,0;optyp '';calcode '';freqid 1;go
查看设置结果：
task 'prtab';default;getn 2
inext 'su';go

8.TASAV

9.run vlaprocs
校正过程:VLACALIB(天线被校正)VLACLCAL(源被校正:把对天线的响应校正外推到对目标源的观测)VLARESET()
初级校正源对天线的校正：
task 'vlacal';
cals '3C286',''
calcode '*'
refant 23
uvrange 0 0
snver 1
docalib -1
doprint -1
minamp 10
minph 10
inp
vlacal              (运行完之后，顺便检查是否有异常的天线EDITA(LISTR-UVFLG, TVFLAG)，对于VLA的观测，并不是很必要。)
(SNPLT看closure error，然后决定是否要重新FLAG并重做CALIB)
task 'edita';default;getn 2
inext 'ty'
inver 0
timer 0
freqid 3
bif 1;eif 0
antennas 0
antuse 1,2,3,4,5,6,7
flagver 1
solint 0
dohist false
dotwo true
crowded true
inp;go

次级校正源对天线的校正(添加在SN中):
tget 'vlacalib'
cals '2229+695' ''
calcode '*'
refant 23
uvrange 0 0
snver 1
docalib -1
minamp 10
minph 10
doprint -1
inp
vlacal
运行aips算出次级校正源的流量值:
task 'getjy'
source '2229+695',''
cals '3C286', ''
go
查看流量信息:task 'prtab';getn 2;inext 'su';go
对目标源进行校正:
task 'vlaclcal'
getn 2
source 'CEPA'  '2229+695'  ''
cals '2229+695' ''
opcode 'cali';timer 0;interp '2pt'
(SMOTYPE ''
snver 1
gainuse 2
)
doprint -1;vlaclcal
初级校正源的SN表中信息写入Cl表中:?SN table    1 has already been applied
wait;tget 'vlaclcal';sour '3C286' '';cals '3C286' '';inp;
vlaclcal
查看校正后校正源的振幅和rms:
task 'listr';optype 'matx';getn 2
sour '2229+695' ''
docal 2
docrt 132
dparm 3 1 0
uvra 0
anten 0
baseli 0
bif 1;inp
go
tget listr;bif 2;inp;go
［重做:inp VLARESET;VLARESET
删除表:task 'extdest';getn 2;(SN表);inext 'sn';inver -1;extdest;(CL表);inext 'cl';inver 2;extdest］

10.查看校正效果:
task 'uvplt';inp
sour 'CEPA'  ''
dotv 1
calc ''
xinc 10
bparm(1) 11
docal 2
bif 1;inp
go
bparm 6 7 2 0;go

11.生成只含单个校正后的源的uv数据文件:
task 'split'
getn 2
source 'CEPA'  '2229+695' ''
calcode ''
uvrang 0;timer 0;
docalib 2;flagver 1
gainuse 2
dopol -1
doband -1
bif 0
eif 0
stokes ''
blver -1
aparm 0
douvcom 1
ichansel 0;inp
go split;pcat

12.图像处理:
CLEAN过程:
pcat;getn 6
cellsize和imsize的推荐值:task 'setfc';inp;go
task 'imagr';
docalib -1;NBOXES 0;
先用低分辨率查看一下:niter 0;cellsize 1 1;imsize 512 512;robust -5;dotv 1;go
正式CLEAN:
niter 500;cellsize 0.15 0.15
imsize 512 512
robust -5
dotv 1;go
框出源所在的位置:TVBOX
检验CLEAN效果:IMEAN,IMSTAT
task 'imean';pcat;getn 11;
dohist 1;blc 0;trc 0;nboxes 50;go
删除图像文件重做？
imh查看综合波束大小：conv size=
pcat;for I=6 to 11;getn I;zap;end
作等高图:
pcat;getn 11
task 'kntr';
Docont 1
Dogrey -1;Dovect -1
xinc 1;yinc 1;plev 0
clev   5.6169E-05 （使用imean找到rms值）
levs -3 3 6 9 12 18 24 30 36 45;inp
task 'kntr';dotv 1;go
tvwin，tvini，tvall联合使用来选择适当大小的窗口
blc 0;trc 0选择整个窗口
go
task 'vplot';getn 6  UV见度文件
get2n 9 clean之后的图像文件
nmaps 1;nplots 8;aparm 0;bparm 0;go
bparm 0 2 0;go
输出图像结果:
task 'kntr';dotv -1;go
imh查看是否新增了一个PL表:
task 'lwpla';outfile 'PRTFIL:iras22543+6145.ps';inp         文件保存在/home/sxj/AIPS/PRINT
go             lwpla对此文件的读写方式是追加。
tvall;tvwin;task 'imean';
数据拟合:
tvall;tvwin;imstat
task 'imfit';inp;go(记录peak位置，大小，流量等信息)
生成fits文件，很多工作可以以此为基础:
task 'fittp';dataout 'FITS:iras22543+6145.fits';inp         文件保存在/home/sxj/AIPS/FITS
go             fittp对此文件的读写方式是生成新文件，故文件已存在的话会报错。

附录：
删除文件:for I=6 to 11;getn I;zap;end;pcat;
删除表:task 'extdest';getn 2;inext 'sn';inver -1;extdest;inext 'cl';inver 2; extdest;imh
删除项目:
删除信息:
删除用户:

............................................................................................................................................
自校正过程:imagr, calib, imagr, calib...
1.
UVSRT(自校正要求表的序列为"TB")

............................................................................................................................................
查看磁盘使用量:
cd $AIPS_ROOT/DATA
export MYDATA=`pwd`
cd ~
aips tv=local

user 32000;indisk 0;go disku

实时帮助:
help, about, explain, apropos