仰望星空,脚踏实地分享 http://blog.sciencenet.cn/u/duancg

博文

CLASS使用部分命令

已有 978 次阅读 2022-5-5 17:16 |个人分类:CLASS|系统分类:科研笔记

CLASS手册(http://www.iram.fr/IRAMFR/GILDAS/doc/pdf/class.pdf

blob.png


一、将望远镜观测的.fits文件转换为.bur文件

LAS> clear                           

LAS> set angle sec                   

LAS> set format long                 

LAS> set unit v f                    

LAS> set plot histogram              

                                     

LAS> define integer n1 j             

LAS> define char line1*18            

LAS> define char line1s*18           

LAS> define char source1*18          

LAS> define char bank1*18            

                                     

LAS> let line1 SiS                   

LAS> let bank1 A                     

LAS> let line1s w1                   

                                     

LAS> sic delete alllines.bur         

LAS> file out alllines.bur single    

LAS> sic find *'bank1'_'line1s'*.fits

LAS> say 'dir%nfile'                 

LAS> for i 1 to 'dir%nfile'          

       fits read 'dir%file[i]'       

       modify line 'line1'           

       modify doppler 0              

       write                         

LAS> next                            


二、将总的.bur转换为各源的.bur

LAS> file in alllines.bur               

LAS> find                               

LAS> list                               

LAS> set unit v                         

LAS> sic delete DR21_SiS.bur            

LAS> file out DR21_SiS.bur single       

LAS> find /source DR21 /line SiS        

LAS> average /resample                  

LAS> write                              


三、减基线

LAS> sic delete xxx.base     

LAS> file out xxx.base single

LAS> file in xxx.30m/xxx.bur 

LAS> find                    

LAS> get f                   

LAS> plot                    

LAS> set mode x 1111 1111    

LAS> set win 111 111         

LAS> draw win                

LAS> base 1/2/3 /plot        

LAS> plot                    

LAS> write                   


四、拉平基线

LAS> define real ty /like ry  !定义一个实数数组ty,让它和ry有一样的维数,ry存储的是谱线的强度 

LAS> set window -200 -42.3 -22 0 90 100  !窗口分别设置为基线部分、非基线部分            

LAS> base 2 /pl                                                      

LAS> let ty ry /where rx.gt.XXX  !(横坐标:速度或者频率,可自己根据set unit v f或者f v设置) ,将速度/频率大于和小于XXX的部分减过基线,也就是基线拉平到0线,数组ry存储这部分 

LAS> set win aaa bbb ccc ddd    !两个窗口分别设置为基线部分、非基线部分               

LAS> base 0 /plot                        

LAS> let ry ty /where rx.gt.XXX    !用ty将将大于XXX的部分替换掉                    

LAS> plot                                                                                          

五、谱线证认

LAS> file in xxx.base             

LAS> find                         

LAS> get f                        

LAS> plot                         

LAS> use in cdms/jpl/xxx.db       

LAS> lid                          

LAS> type xxx.mod                 

LAS> modsource xxx.mod 30 /verbose       !30为望远镜口径

LAS> modshow                      


六、创建新的数据库

LAS> type xxx.cat     

LAS> type partfunc.cat

LAS> use in xxx.cat   

LAS> select           

LAS> use out xxx.db   

LAS> insert


七、高斯/SHELL拟合 

LAS> modsource xxx.mod 30 /verbose

LAS> modshow                      

LAS> set unit v          !gauss拟合用v,shell拟合用f         

LAS> set win                      

LAS> base 0 /p                    

LAS> method gauss         !shell拟合用method shell    

LAS> line 1/2/3                   

LAS> mini                         

LAS> vis                          


八、删除坏通道

LAS> draw    !draw坏通道的左右两边,得到2个通道编号channel number

LAS> draw kill 7000    !标记第一个坏通道                                  

LAS> draw kill 7100    !标记第二个坏通道                                  

LAS> draw fill 7000    !填充第一个坏通道                                  

LAS> draw fill 7100    !填充第二个坏通道                                  

LAS> p


九、两个bur合成一个bur                          

LAS> sic delete 1.bur     

LAS> file out 1.bur single

LAS> file in xx1.bur      

LAS> find                 

LAS> copy                 

LAS> file in xx2.bur      

LAS> find                 

LAS> copy                 


十、3sigma线

LAS> greg1\pen /col 0 /dashed 2        !col 0是颜色,0、1、2等;dashed 2是虚线线型,1(实线)、2、3等,一般用col 0 和dashed 2

LAS> greg1\draw relocate -12 0.84 /user   !-12是用greg1\draw命令,使用鼠标右键选择图形左边缘后,命令框中的User Coordinates         

LAS> greg1\draw line 17 0.84 /user      !17是图形右边缘的User Coordinates,0.84为自己计算的3rms值                            


十一、调整class主界面字体

LAS> greg1\set char 0.85   !0.85为字体大小 



https://m.sciencenet.cn/blog-3489179-1337175.html


下一篇:[转载]手把手教你看懂天体的光谱

0

该博文允许注册用户评论 请点击登录 评论 (0 个评论)

数据加载中...
扫一扫,分享此博文

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2022-12-2 21:56

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部