科学网

 找回密码
  注册
找源的算法(工具)
钱磊 2010-11-10 19:44
这段时间一直在用CUPID里的gaussclumps找分子云中的云核。一般来说可以找到大部分椭球状的结构,但是有时候也可以看到一些十分扁的椭球(投影到平面上就是一个长条),也有很多地方无法拟合。也就是说分子云中除了规则的椭球状云核,还有其它一些不规则的结构,比如条状结构和片状结构。这些结构用gaussclumps是无能为力的 ...
个人分类: 总结|3576 次阅读|没有评论
矮星系中的气体
钱磊 2010-11-8 16:19
矮星系是一种质量典型质量 到 太阳质量的小星系,它们中的很多是如银河系这样的大星系的卫星星系。虽然小,但是这些星系对于暗物质的观测研究有重要的作用。 通常像银河系这样的星系,虽然暗物质多于可见物质,但是可见物质对引力也有可观的贡献,尤其是星系的中心部分,而中心部分的密度分布对于区分暗物质模型十分关键。 ...
个人分类: 总结|3906 次阅读|没有评论
运行简单的IDL脚本
钱磊 2010-11-8 14:28
没系统学过IDL,基本上是看着别人的程序学的,所以时常有一些最基本的知识都不知道,直到又看到新的例子。以前写一个IDL脚本总是一板一眼PRO xxx开头,END结尾,不管这个脚本如何简单,运行的时候还每每 IDL.com xxx 然后 IDLxxx 今天看到一个算暗晕质量函数的程序包,主程序是Fortran写的,画图用IDL。那个IDL脚本 ...
个人分类: 总结|16904 次阅读|没有评论
用IDL从Fits文件中读二进制表
热度 2 钱磊 2010-11-4 17:23
一个程序计算的结果会输出到一个.FIT文件中,以前用软件打开看过,就是一个数据表。但是用IDL尝试了很多此,都读不出来,总是提示没有图(image)。于是以前都是手动把数据表搞出来。 最近数据量越来越大,手动已经可以把人烦死了。于是迫不得已再开始尝试用IDL读这种文件。今天的尝试比较顺利,顿悟了一下,终于知道这个 ...
个人分类: 总结|8247 次阅读|6 个评论 热度 2
CUPID FINDCLUMPS 拟合分子云核(三)
钱磊 2010-10-29 21:02
我需要知道FINDCLUMPS输出的云块的线宽是怎么定义的。于是查看findclumps.c中的语句 cupidStoreClumps( OUTCAT, xloc, ndfs, nsig, deconv, backoff, ishape, velax, ...
个人分类: 总结|3516 次阅读|没有评论
编程点滴(二)数据怎么存?
钱磊 2010-10-29 10:54
我的程序要输出一些数据,参数文件也要改,原来都是手动,算完了就建个文件夹,把数据文件和参数文件拷进去。久之就烦了。昨天看到一个别人备份系统的脚本,里面用时间作文件名,这样就一劳永逸了。我觉得这是个好办法。因为每次执行程序的时候时间都不同,文件夹就不同了,把文件拷到这些文件夹里就稳妥了。 用时间做文件 ...
个人分类: 总结|2760 次阅读|没有评论
编程点滴(一)搜索变量名
钱磊 2010-10-29 10:28
1. 为了以后搜索变量名方便, 给变量赋值的时候等号前后最好各加一个空格。 2. 对于不区分大小写的语言也最好统一变量名。 3. 如果一个程序很大,各子程序分布在不同的文件(比如都是. c文件)中,要在所有文件中查找 一个变量(比如找sum变量),可以试试 grep &n ...
个人分类: 总结|3211 次阅读|没有评论
CUPID FINDCLUMPS拟合分子云核(二)
钱磊 2010-10-28 22:36
为了弄清楚输出文件中的宽度和峰值时本征的值还是经过波束平滑的值,需要看一下输出程序中是如何处理的。在cupidgclistclump.c中可以看到 msgSetd( V, par *rms ); msgOutif( MSG__DEBUG, , Peak intensity: ^V, status ); msgSetd( V, par ); msgOutif( MSG_ ...
个人分类: 总结|2874 次阅读|没有评论
CUPID FINDCLUMPS拟合分子云核(一)
钱磊 2010-10-28 21:17
用CUPID里的FINDCLIMPS拟合分子云中的团块的时候发现拟合出了很多虚假的团块,在光谱里看的时候并没有相应的高斯成分。所以决定看看源代码,做点笔记。CUPID乃至整个Starlink的源代码都可以在此获得: http://starlink.jach.hawaii. edu/git/?p=starlink.git;a= tree;f=applic ...
个人分类: 总结|5015 次阅读|没有评论
用IDL生成某个分布的二维数组
热度 1 钱磊 2010-10-27 14:27
假定把一个多维(这里说二维的)数组的指标看作坐标,要让每个元素表示某个函数在相应坐标的值应该怎么办?循环当然是可以的,不过其实没有必要,因为这些坐标都是有规律地排列的,IDL自有处理的方法。 对于一维数组而言,IDL自然就可以生成一个数组,每个元素的值就是其指标值。 IDL a = lindgen(n) 有了这样一个数组 ...
个人分类: 总结|10273 次阅读|1 个评论 热度 1

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

GMT+8, 2024-4-24 15:56

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部