用VERAーー源観測(cè)ー

1、AIPS automatic pipeline 用VERA源観測(cè)像合成2007年9月18日今井裕/原産地（相関処理局）工場(chǎng)（解析者）原材料（DVD中）輸送経路/行/観測(cè)生及較正用補(bǔ)助集約/中身各種図生成/基本的較正及結(jié)果表示/処理較正（今日実習(xí)）解析分析（研究）実習(xí)手順/作業(yè)環(huán)境/原材料入手 (FITS files, calibration files)/実行/AIPS手作業(yè)/VERA実行/Velocity tracking (SETJY, CVEL).手作業(yè)/Reference channel選択flagging.手作業(yè)/MCAL実行/Image cube (X,Y,V)作成表示.手作業(yè)結(jié)局

2、、VLBI像合成原理理解、正image cube作成、自分判斷：像描。NML TauriSiO maserVERA観測(cè)諸元詳 r05017a.vex 見(jiàn)/観測(cè)提案者研究進(jìn)行責(zé)任者(principal investigator) ：清水理絵松本尚子（論文執(zhí)筆研究遂行責(zé)任持）/観測(cè): r05017a (2005年通日17日行最初観測(cè)）/研究目標(biāo)天體： 型変光星 NML Tauri(中NML-Tau公式名稱)/観測(cè)電波輝線： 一酸化珪素(SiO) v=1, v=2 (J=1-0)(実習(xí)、振動(dòng)勵(lì)起狀態(tài) v=1 ）/使用電波望遠(yuǎn)鏡：VERA局望遠(yuǎn)鏡/観測(cè)周波數(shù)幅： 16MHz （視線速度幅110 km/

3、s）/分光點(diǎn)數(shù) 512點(diǎn) (P.I.相関処理時(shí)指定）見(jiàn)方/“ ”中身任意。codename書(shū)言、実際”codename”馬鹿正直入力。観測(cè)名入。/環(huán)境変數(shù)（directory名）/$FITS FITS名前環(huán)境変數(shù) (in LINUX)/FITS: 同上 (in AIPS)/入出力： “adverb”AIPS呼indisk, inname, inseq, docal, 作業(yè)環(huán)境/FITS file置場(chǎng)所： $FITS(default: $FITS=$AIPS_ROOT/FITS）/観測(cè)情報(bào)較正用保管場(chǎng)所： $INexport IN=/home/user name/calibs/ （最後”/”忘）

4、/解析結(jié)果表示(TEXT/PostScript files)出力先: $OUTExport OUT=/home/user name/plots/.bashrc (or .cshrc) 書(shū)込便利DVD指定良。 (e.g. $DVD=/Volumes/CDROM/)/$IN, $OUT指定directory作成観測(cè)情報(bào)入手（、？何起？)/観測(cè)指示書(shū)(VEX file)/sftp vlbiveraserver.mtk.nao.ac.jpuser ID: vlbi passwd 誰(shuí)聞！/cd ./mmmyyyy (e.g. ./mar2005 ./apr2005 ./may2005)/get ryyd

5、dd.xxx yy=year 20yy, ddd=Day of the year (DOY)xxx: =vex: VEX file(本） =CHECK: （天體仰角、等）e.g. r05084b.vex r05116b.vex r05151a.CHECK/観測(cè)実施狀況：: VERA文章/./.全局狀況：保持者問(wèn)合相関生及較正入手/普通DVD送 所定directory（後述）/天文臺(tái)(133.40. )接続/ftp mtksp1 (51)user ID: vlbi passwd: 誰(shuí)聞！/home/fits1002 DIR1K, fits2000 DIR2K/*.FITS

6、.1 : FITS data/*.TBIN.1111 : phase calibration/home/work1/analyfiles/Tsys/obscode/station/obscodestation.tsys or obscode_Tsys_station.dat /home/work1/analyfiles/calib2B /obscode/station/ obscode_2BCAL_station.dat 相関生及較正加工（）/FITS files ($FITS 內(nèi)) file名変更 (実行) fitsname: xxxxx.fitsN (N=1,2,3)/AIPS読込用Tsy

7、s files作成: veratsys 作成 file $IN /2-beam calibration (SN) table noiseana 作成 file $IN 相関生及較正加工（2）/Flag file (in $IN)：削除 cp ch_spectral channel number.chflag obscode.uvflg/ obscode.uvflg中身 (band edge ）! The format is appropriate for UVFLG in AIPS.! The outer channels of each IF are flagged.opcode = FLA

8、G dtimrang = 0 timeoff = 0antennas=1 timerang= 0,0,0,0,400,0,0,0 bif=1 eif=16 bchan=1 echan=8 reason=Channels at band edges/antennas=1 timerang= 0,0,0,0,400,0,0,0 bif=1 eif=16 bchan=1017 echan=1024 reason=Channels at band edges/（以降局?jǐn)?shù)分）Pipeline script 使方/RUN scripts: JNET.001, VERA.001, MPLOT.001, MC

9、AL.001, MFIT.001, MMAP.001, SQASH2, MOMNT2, IMGSAV.001, ATM.001/RUN scripts in /usr/local/aips/AIPS version/RUNAIPS version: e.g. 31DEC06/HELP files (*.HLP) in /usr/local/aips/AIPS version/HELP/Read a RUN script: RUN VERA/Input adverbs: INP VERA保存版書(shū)込 Copy&Paste代入/adverbs意味確認(rèn)： HELP VERA/Execute the R

10、UN script: VERADemonstration of the pipeline “VERA”/tmask 0/nsrc 6; srclst3C454.30854+2010352+1234C13.210423-013NML-Tau/calsrclst3C454.30854+2010352+1234C13.210423-013/nfits 1; fitsfile FITS:r05017a_v1.fits/doantab 1; indi 2; in2d 2;/refant 2;refan2 1; plotref 2; oint 1/solcal 5; plotavg 0.25;sciter

11、 2;pols ll;/dir IN:r05017a/dir2 OUT:r05017a/obsband Q; doprint=-2; expname r05017a実際./文字列前後 single quotation mark 區(qū)切 小文字大文字変換/同行幾 adverbs指定： ; 代入文/実習(xí)使. FITLD終了 全 adverbs 再入力 + tmask=-2 MSORT再開(kāi)/最後終 必要較正完了済。心配。r05017a .MSORT . 1 必要。 殘消： getn catalog#; zap生成確認(rèn) (VERA)/R05017a.SCAN/R05017a.DTSUM/R05017a.

12、PRTAN/R05017a.AUTOCORR/R05017a_POSSM(or VPLOT)_UNCAL(or CAL).PS/R05017a_BDTPOWER.PS/R05017a_ACCOR.PS/R05017a_APCAL-GAIN(or GAIN).PS/R05017a_FRING_DELAY(or RATE, PHAS, SNR).PSPostScript (PS) filesGhostview (gv) 図示図示Velocity tracking for a maser source/SETJY: define a source velocitytask setjy;indisk

13、2;getn catalog#; source NML-Tausysvel=32.0; restfreq 4.312e10, 2027000; veltyp lsr;veldef radio;aparm 255 0/CVEL: velocity trackingtask cvel; doband 1; bpver 1; gainuse 1; flagver 1;outdisk 1;getn catalog#; geto catalog#; outclass cvel; SETJY/CVEL付近付近OK=左右大左右大源較正Reference (velocity) channel選択/選択條件/短

14、積分時(shí)間検出 cross-power spectrum 確認(rèn)/輝度分布構(gòu)造単純 相関振幅時(shí)間変化単調(diào) 閉口位相(closure phase)近/両隣spectral channel 輝度分布単純 visibility 位相一定値保 cross-power spectrum確認(rèn)表示tget possmgetn 26; dotv 1doband=-1bchan 200;echan 300Interactive flagging（余裕人）/SPLIT: splitting only the reference channeltask split; getn CVEL cat#; source NML

15、-Taubchan 263; echan 264; docal 1; gainuse 6; doband=-1outclass split; aparm 2 1 0/MULTI + INDXR/IBLED: interactive flaggingtask ibled; getn MULTI cat#; flagver 1; bchan 1; echan 0; antennas 0; baseline 0; docal=-1; doband=-1; dparm 0 0 0 10 10source ; stokes iTV Server畫(huà)面橫畫(huà)面橫/TACOP: copy a FG table

16、to original (u,v) datatask tacop; getn new catalog#; inext fg;inver 1;geto CVEL file cat#; outver 1;ncount 1Demonstration of the pipeline MCALtmask 0; nsrc 1; srclst NML-Taucalsrc NML-Tauindi 2; in2d 2; refant 2; refan2 1; plotref 2; refif 1; refbchan 263; refechan 264; bcl 8; bsn 6doband=-1; pols L

17、L; plotavg 0.1; solcal 2; solsub 2; oint 1 snrcut 5; sciter 4; dir2 OUT:r05017a/expname r05017aobsband Q; doprint=-2生成確認(rèn) (MCAL)/R05017a_FRING_MASER_RATE (or SNR).PS/R05017a_VPLOT_MCAL.PS/r05017a_NML-Tau_ICLN (or IMAPU, IMAPN).PS/r05017a_NML-Tau_EL (or UVCOV, UVPLT).PS/r05017a_NML-Tau_VPLOT_MODEL.PS/

18、r05017a_MCAL_PHAS1(3, 5, or 7).PS/r05017a_MCAL_AMP2(4, 6, or 8).PS/r05017a_MCAL_TXTn.TXT他！TXINC 2; TYINC 2TVALLTVLABEL場(chǎng)所確認(rèn)TVBOXMCAL直/全 adverbs 再入力/次 adverbs 指定 (CLEAN box掛再補(bǔ)正）tmask=-5nboxes 2clbox 501 505 521 521 789 605 807 617 （左下及右下座標(biāo)、）Search maser emission in spectra/IMHEAD使用CL table確認(rèn)/存在velocit

19、y channel 確認(rèn)tget possmdotv 1; solint 400; gainuse 11; bchan 150; echan 350SPLIT: Splitting necessary velocity channelstget splitgetn CVEL cat#; source NML-Taubchan 236; echan 275; docal 1; gainuse 11; doband=-1; flagver 1; outclass split; aparm 0IMAGR: making an image cube I(x,y,v)task imagr; getn S

20、PLIT cat#;docal=-1; doband=-1; flagver 0; source NML-Tausmooth 0; stokes I; nchav 1; chinc 1; channel 0bif 1; eif 0; bchan 0; echan 0; 番號(hào)、開(kāi)始1,2,3,改割振outname r05017a_v1; outdisk 1; outver 1; clr2n; outver 0cellsize 0.0001 0.0001; imsize 1024 1024; fldsize 0; rashift 0; decshift 0uvwtfn na; gain 0.05;

21、 dotv=-1; niter=1000; flux=-1.0; bmaj 0; bmin 0; bpa 0; clbox 0; nboxes 0SQASH2: Making a velocity-integrated imagerun sqash2getn image cat#; xg 0; yg 0;nedge 10; bvel 0; evel 0; snrcut 7dir OUT:r05017a/tget kntrblc 350 400;trc 850 900;dotv 1VLBI処理講習(xí)會(huì)學(xué)/複數(shù)電波望遠(yuǎn)鏡同時(shí)使天體観測(cè)、後日処理詳細(xì)電波像得/可視光赤外線観測(cè)方、畫(huà)像出來(lái)分易 （単発想）/合成像畫(huà)質(zhì)：処理者力量決（、主観決）/80%情報(bào)処理（雑音入時(shí)系列処理）、19%自然科學(xué)（地球自転公転、空気）1%天文學(xué)（源関知識(shí)）/AIPS task 途中止/“abort” 右側(cè)status: WRIT or READ（書(shū)込禁止?fàn)顟B(tài)）/“clrstat”: clear of status (READ, WRIT)/Dowait=1場(chǎng)合 : /LINUX上process id探(ps)/kill process id AIPS上“clrstat” /解析方法間違/display見(jiàn)下言質(zhì)問(wèn)法、時(shí)間消耗得殆/何起 (visi