gildas数据处理完整汇总(待续)

综合几个教程，做出一个完整的汇总。

---

假设要处理的数据文件名为：mydata.30m

class
step 1:
file in mydata.30m
find
list
会显示大量数字，其中有source name,linename,telescope,offset  in RA and DEC, telescope observation number.
step 2:提取合适的光谱
set source sourcename
set line  linename
set tel telescopename
set offset n1 n2
find
从大量数据出筛选出合适的光谱来。
list
会看到只剩下少部分光谱。
get first
plot
看看第一个光谱的样子
set unit v f
将速度作为x轴坐标值
注意读取头文件信息：RA,DEC,tau(大气不透明度)，tsys(system nosie temperature),time(on-source integration time,in minutes),El(elevation of source海拔)，DV(velocity resolution per pixel)
step 3:set nomatch
average
plot
drop n  删除第n个谱线
average
plot

smooth  该命令每运行一次光谱分辨率减小一半。
line 2
min
vis
做高斯拟合，之前应该先做基线处理
做完高斯拟合会出现一行数据
分别是:area,position,width,intensity
这些数据很有用，area是线的积分强度，position是中心位置，width是线宽
上面的数可以用来计算H2质量，红移（距离），星系中气体动质量(kinematics)
计算公式不在这里给出了，见网页http://loke.as.arizona.edu/~ckulesa/smt/basics.html
step 4 基线剔除
所有的光谱都有自己的基线，移除设备产生的干扰。
file in mydata.30m
find
set win m1 m2
file out reduced.dat m
for i 1 to found
   get next
   base 1
   plot
   write
next

用上面的循环得到基线处理过后的文件reduced.dat

file in reduced.dat
find
set mode x n1 n2
set mode y n3 n4
map /gr
可以看到一个塞满大量光谱的图
下面要做一件很重要的事情，提取出数据中的重要信息，做出moment maps，是展示数据的最好的方法，用到的命令是print moment
print moment n1 n2 /output moments.txt    使用print area也可以做到
n1 n2是选取的速度区间

step 5:
下面在greg中进行
生成的moments.txt里有大量有用的信息，以列来划分
包含：observation number,x and y offsets,integrated intensity,centroid velocity,linewidth.
column x 2 y 3 z 4 /file moments.txt
第二列给x轴，第三列给y轴，第四列给z轴，2和3是坐标RA,DEC的offsets,4是积分强度.
下面的命令取自G53

symbol source "sourcename"
symbol contourlevels "20 to 100 by 10"
symbol range "300 -300 -300 300"

symbol channel "dataname.mean"
symbol spitzer "dataname.gdf"
set plot_page landscape  画出图的尺寸，默认landscape即A4大小
set orientation 0
set character 0.4
set font duplex
set system equa
tick 2 10 24 120
set tick 0.15

set character 0.25
set box 9 15 9 15
image 'channel'
limits 'range' /rg
pencil /colour 7 /weight 1 /dash 1
lev 4 to 16 by 3
rgmap /per 5.5 /pens 0
pencil /colour 0 /weight 1 /dash 1
box /abs

def image c 'spitzer' read
rgdata c /var
pencil /colour 0 /weight 1 /dash 1
pl /sc linear -0.18 1.1
lut rainbow3

set char 0.18
pencil /colour 0 /weight 0 /dash 1
wed top 0.25